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FOREWORD 

This  report  has  been  prepared  by  the  Elkton  Division  of  Thlokol  Chemical 
Corporation  for  the  Air  Force  Systems  Command,  Rocket  Propulsion  Laboratory 
(DGSMD),  at  Edwards,  California.  Lt.  Harold  W.  Gale  is  the  Air  Force  Project 
Officer. 

The  report  number  assigned  to  this  document  by  Thlokol  Chemical  Corporation 
is  E92-63.  The  principal  contributors  to  this  report  are  Messrs.  W.  G.  Andrews, 
D.  H.  Frederick,  F.  E.  Mcore,  D.  Saylak,  A.  Stornelli,  D.  D.  Thomas,  and  R.  H. 
Thompson.  Mr.  James  F.  Hoebel  is  the  Program  Manager. 


ABSTRACT 


The  work  conducted  under  Contract  AF  33(616)-6530  during  the  period  from 
July  1,  1982,  to  July  31,  1863,  is  reported.  This  program  wa»  a  research  study  to 
advance  the  state-of-the-art  of  solid  propellant  grain  design  and  consisted  primarily 
of  theoretical  and  applied  analyses  of  solid  propellant  grain  designs,  Heat  transfar 
factors,  effect  of  strain  on  burning  rate,  and  analyses  of  propellant  stresses  were 
Specifically  considered. 

An  analytical  study  was  undertaken  to  develop  a  two-dimensional  heat  transfer 
analysis  for  solid  propellant  motors  subjected  to  a  nonlinear  circumferentially  variant 
and  time  dependent,  external  thermal  environment,  including  convection,  radiation, 
or  known  heat  flux,  or  any  combination  of  these. 

Finite  difference  equations  representing  the  irregular  external  boundary  condition, 
the  temperature  of  the  case,  the  case-propellant  interface,  the  propellant  interior 
points,  and  the  heat  flow  across  the  Irregular  internal  boundary  were  developed.  These 
equations  were  programmed  and  Incorporated  into  a  numerical  solution.  The  program 
further  considers  temperature  dependent  thermal  properties  and  the  internal  heat 
generation  due  to  the  exothermic  polymerization  reaction  of  propellant. 

A  study  was  performed  to  evaluate  the  apparent  approaches  to  solid  propellant  grain 
design  through  digital  techniques  and  then  to  implement  the  best  approach  into  a 
FORTRAN  language  computer  program.  The  Elkton-Moore  Method,  a  Vectorial 

i  » 

approach  to  the  geometrical  portion  of  the  problem,  was  conceived  and  developed  and 


ill 


has  served  as  the  theoretical  base  for  the  Advanced  Grain  Design  Analysis  program. 
Program  AGDA  evaluates  motors  with  arbitrary  center  perforations,  with  or  without 
head-  and  aft-end  webs,  and,  in  addition,  considers  two-dimensional  temperature 
variations  existing  in  the  grain  and  accounts  for  their  effect  on  propellant  burning  rate. 
Motors  are  described  to  the  program  in  terms  of  readily  available  parametric  values, 
which  are  the  standard  configuration  dimensions  stated  on  all  fully  loaded  case  drawings. 

Four  propellants,  including  two  hydrocarbons,  one  polyurethane,  and  a  polyvinyl 
plastisol  type,  were  subjected  to  tensile  strains  ranging  from  0  to  12  percent  and 
burned  in  the  strained  condition  in  a  strand  burner  at  pressures  ranging  from  500  to 
1500  psi.  Results  established  that  not  all  propellants  are  affected  to  the  same  degree 
by  strain.  A  correlation  between  percent  change  in  burning  rate  and  volumetric  change 
index  is  given.  The  latter  term  is  a  representation  of  the  departure  from  Incompressibility 
that  the  material  experiences  due  to  st>  _n,  't'o  achieve  these  data,  an  optical  technique 
for  measuring  Poisson's  ratio  was  developed  and  is  described  in  this  report.  In 
addition  to  linear  burning  rate,  the  maOi  «._;rning  rate  of  the  four  propellants  was 
investigated. 

A  study  was  undertaken  to  find  a  method  by  which  the  effects  of  internal  pressurization 
thermal  shrinkage,  and  axial  acceleration  loads  upon  propellant  structural  stability 
could  be  determined.  Nominal  effort  was  expended  in  specifying  a  means  for  analytically 
characterizing  the  viscoelastic  properties  of  a  solid  propellant  b>  means  of  linear 
differential  operators  and  the  time  dependent  elastic  modulus,  E(t).  The  method  and  its 


iv 


accuracy  u^e  illustrated.  Major  8ffort  was  directed  toward  the  development  of  a 
technique  for  structurally  analyzing  axisymmetric  grains  having  arbitrary  end 
geometries  and  straight-through  ports.  The  analysis  is  formulated  in  terms  of  two 
stress  functions,  $  and  Y  »  and  a  coupled  pair  of  elliptio  partial  differential 
equations.  The  pertinent  boundary  conditions  are  specified,  the  solution  by  means  of 
finite  difference  equations  is  outlined,  and  a  corresponding  computer  program  is 
discussed, 
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The  objective  of  this  program  has  been  to  advance  the  state-of-the-art  of  solid 


rocket  propellant  grain  design  from  a  thermal-structural  as  well  as  a  ballistic  stand¬ 
point,  Previous  studies  under  this  contract  have  included: 

1)  The  development  of  an  IBM-650  computer  program  to  compute 
theoretical  burning  surface  area  as  a  function  of  time. 

2)  The  adaptation  of  photoelastic  principles  to  determine  experi-  . 
mentally  strains  in  scale  rocket  grains  under  thermal  loading. 

3)  The  development  of  an  analytical  method  of  calculating  the 
temperature  at  any  point  in  a  rocket  grain  subjected  to  motor 
curing,  environmental  conditioning,  or  aerodynamic  heating. 

4)  A  mathematical  study  of  the  combined  internal  and  external 
loads  imposed  on  a  system  of  concentric  elastic  cylinders  to 
determine  the  contribution  of  the  propellant  to  the  reduction 
of  the  effect  of  chamber  pressure  on  the  motor  case  wall. 

5)  The  development  of  a  generalized  three-dimensional  computer 
program  to  evaluate  the  ballistic  properties  of  internal  burning, 
singly  connected  solid  propellant  grains,  including  consideration 
of  head-  and  aft-end  web  effects  on  the  burning  surface. 

6)  The  development  and  application  of  the  PhotoStress  technique 
as  a  two-dimensional  experimental  method  for  propellant  strain 
analysis  through  the  extension  of  the  previous  photoelastic  -  . 
studies.  This  effort  at  Thiokol- Elkton  represented  the  first 
application  of  the  PhotoStress  technique  to  solid  propellant 
grains  in  the  industry. 
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1)  The  analysis  of  the  effect  of  propellant  defects  on  ballistic 
performance  through  window  bomb  and  motor  tests,  High¬ 
speed  motor*  pictures  were  utilized  to  analyse  the  effects, 

$)  A  heat  transfer  analysis  in  which  numerical  and  electrical 

analog  technique!  were  Used  to  simulate  the  thermal  Conditioning 
history  of  a  motor  grain, 

Th«  current  program  has  been  divided  into  three  interdependent  aro&a  of  investigation, 
Each  area  Is  described  in  separate  sections  of  this  reports 

•  Ar&a  A  -  Heat  Transfer  Analysis 

a  Area  B  -  Effect  of  Strain  on  Burning  Bate 

•  Area  C  -  Propellant  Slump  Analysis 


II.  HEAT  TRANSFER  ANALYSIS 


A.  TWO-DIMENSIONAL  HEAT  TRANSFER  ANALYSIS 

The  application  of  solid  propellant  rocket  motors  to  system  requirements 
where  the  motor  is  subjected  to  nonuniform  environmental  conditions,  either  natural 
or  induced,  leads  to  many  complex  problems  not  easily  solved  by  routine  techniques. 
For  example,  a  mission  may  be  undertaken  where  the  vehicle  is  subjected  to  heating 
on  a  portion  of  the  circumference  and  cooling  on  the  remaining  portion  depending  upon 
the  attitude  of  the  system  with  respect  to  the  sun,  or  where  variations  in  the  angle 
of  attack  may  cause  circumferential  variations  in  aerodynamic  convective  heating. 
Alternatively,  a  propulsion  system  nestled  on  the  side  of  a  larger  vehicle  could  well 
be  subjected  to  aerodynamic  or  radiation  heating  on  a  portion  of  the  body  and  conductive 
cooling  on  the  alternate  side.  These  and  other  possible  missions  naturally  lead  to 
grain  environments  and  associated  effects  upon  performance,  both  structural  and 
ballistic,  which  previously  were  not  traceable  because  of  the  complexity  of  the  problem. 

Previous  efforts  under  the  original  Air  Force  Contract  AF  33(616)-6530 
resulted  in  the  development  of  a  two-dimensional  (r,  0)  thermal  analysis  for  solid 
propellant  grains  subjected  to  a  uniform  external  environment.  This  analysis  has 
been  extended  under  continued  Air  Force  support  to  consider  nonlinear  circum¬ 
ferentially  variant  and  time  dependent  boundary  conditions  having  environmental 
temperatures,  film  coefficients,  and  other  pertinent  thermal  parameters  of  arbitrary 
magnitude. 
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Further,  the  analysts  will  consider  the  internal  heat  generation  within  the 
propellant  due  to  polymerization.  Thus,  the  analysis  will  provide  the  capability  to 
predict  adverse  temperature  gradients  due  not  only  to  environment  but  also  to  curing, 

Th8  solution  within  the  analysis  to  the  two-dimensional  temperature  field 
existing  in  a  oircular  grain  having  a  star  or  similar  Irregular  Internal  perforate  and 
subjected  to  a  nonlinear  external  boundary  condition  is  obtained  through  the  use  of 
finite  difference  techniques.  Having  represented  the  heat  conduction  equation  and 
boundary  functions  in  finite  difference  notation,  the  solution  is  obtained  by  "marching 
oi'*"  in  finite  time  steps  from  a  specified  initial  thermal  configuration  to  the  time  at 
which  the  temperature  profile  through  the  grain  its  desired.  With  each  time  step 
thermal  properties  and  boundary  parameters  are  adjusted  for  the  change  In  time  and 
temperature. 

This  solution  was  programmed  for  the  IBM-7070  digital  computer  in  the 
FORTRAN  H  language. 

1,  The  Heat  Conduction  Equation 

The  transient  temperature  field  in  a  solid  having  a  circular  oross 
section,  variable  thermal  properties,  and  internal  heat  generation  may  be  expressed 
by  the  following  partial  differential  equation: 


er  _i_  fill  £X  +1  sfc  £1' 

60  “  pCp  Lor  8r  r2  60  60  j 


K 


pCr 


**&  +  !  /9T\+  i_  82Tl 
,Br2  +  r\BrJ  r2  8 0 2 J 


l  +  c- 
cp 


where 


T  -  Temperature,  *R 

0  -  Time,  hr 


(1.1) 
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K  “  K(T)  »  Thermal  Conductivity,  BTU/hr-ft-‘R 

P  “  Density,  lb/ft3 

Cp  *  Cp(T)  *  Specifio  Heat  at  Constant  Pressure,  BTU/lb-*R 
r  *  Radial  Coordinate,  in. 

0  *  Angular  Coordinate,  radians 

The  Internally  generated  heat  q  is  considered  for  this  analysis  to  be  the  exothermic  heat 
of  polymerization  experienced  within  a  solid  propellant  grain  during  curing.  This  may 
be  expressed  as 

q  -  Q[B(1  -P)n  exp  (-A/RT)]  BTU/lb-hr  (1.2) 

where 

Q  3  Heat  of  polymerization  per  unit  mass  BTU/lb-hr. 

P  a  Polymer  Fraction 

K  3  Gas  Constant 

B,A,n  =  Empirical  Constants 

To  account  for  the  extremes  of  thermal  environment  experienced  by 
solid  propellant  rocket  motors  and  to  permit  more  accurate  representation  of  thermal 
properties,  both  the  thermal  conductivity,  K,  and  specific  heat,  Cp,  are  considered  to 
be  functions  of  temperature  of  the  following  form 

f(t)  =  A  +  BT  +  CT2  +  DT3  +  ET4  (1. 3) 

where  each  parameter  has  a  separate  set  of  coefficients.  Tbit'  will  provide  for  an 
extreme  variation  in  thermal  properties  with  temperature. 
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In  order  to  compensate  for  the  variation  in  thermal  properties  during 
cure,  a  relationship  between  cure  temperature  and  the  necessary  parameters  must  be 
obtained  experimentally. 

Having  defined  at  this  point  the  analytical  expression  for  the  heat 
conduction  within  the  propellant,  equation  1. 1  can  now  be  expressed  in  finite  difference 
notation  for  the  single  nodal  point  illustrated  in  Figure  1.  This  produces  the  following 
relationship  in  accdx dance  with  technique  outlined  in  the  literature.  ^ 
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FIGURE  I.  PROPELLANT 


where 


TI,k  " 

T  (r,  0,0  +A  6) 

TJ.k  3 

T  *r,  0,  0) 

TJ+l.k  - 

T  (r+Ar,  0,  6) 

TJ,  k+1  * 

T  <r,  0  +  A0,  0) 

TJ.W  ' 

T  (r,  0  -  A0,  0) 

TJ-l,k 

T  (r-Ar,  0,  9) 

Ar?  and  A0j 

«  Forward  differences  in  r  and  0  respectively 

Aij  and  A0g 

a  Backward  differences  In  r  and  0  respectively 

While  this  equation  primarily  represents  an  energy  balance  for  an  internal  node  it  is, 
at  the  same  time,  the  basic  relationship  for  the  analysis.  By  making  appropriate 
substitutions  within  this  equation  and  by  performing  certain  necessary  adjustments, 
this  relationship  car.  be  adapted,  to  satisfy  the  following  nodal  points  within  the  grain  as 
well: 

1)  External  surface 

2)  Case-propellant  interface 

3)  Internal  surface 

To  illustrate  the  manner  in  which  this  ip.  performed,  each  of  the  above 
will  be  discussed  separately. 

2.  External  Surface  Nodal  Point 

The  boundary  condition  at  the  outer  surface  of  the  grain  resulting  from 
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the  external  environment  to  which  the  motor  is  exposed  is  considered  to  be  any 
combination  of  convoction,  radiation,  or  known  flux.  Since  the  outer  boundary  is 
circular,  the  heat  flow  into  or  from  the  motor  may  be  expressed  in  terms  of  the  radial 
gradient  at  the  surface  alone. 

Thus 

22.  a  Jl£L  (X  -  T  )  +  ^2-  (T  .  T  +  ~  it  tv 

8r  K  '  conv  os'  k  'xrad  Aos  >  k  l1*6) 

where 

*0 

* 

E  3 

a  = 

Tconv 

Toa  = 

^rad  = 


ho  (0,0)  3  Outer  Boundary  Film  Coefficient, 

BTU/hr-ft2-*R 

^  (0,8)  ■  3  Gray-Body  Shape  Factor  (Dimension¬ 

less) 

E  (0,  0)  *  Heat  Flux  per  Unit  Time,  BTU/hr-ft2 

Stephan  Boltzmann  Constant,  BTU/hr-ft2-*!!4 

TCOnv$*®)  =  Environmental  Convective 

Temperature,  °R 

Tos(0,0)  =  Outer  Surface  Temperature,  °R 

Trad(*>  9)  -  Environmental  Radiant  Temperature, 


To  simulate  a  nonlinear  circumferentially  variant,  external  environment, 
about  all  or  a  portion  (thermally  symmetric)  of  the  motor,  each  of  the  principal 
boundary  parameters,  h0,  Tconv  Tracp  and  E,  are  considered  to  be  functions  of 
angular  position  0  of  the  following  form 

f(0)  =  a+  b0+c02  +  d  03  +  e  04  +  f  sin  0  +  g  cos  0  (1, 6) 
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Each  parameter  naturally  has  a  separate  set  of  coefficients.  Furthermore,  to  give 
added  versatility,  the  outer  surface  of  the  motor  may  be  divided  into  two  separate 
"zones"  for  which  each  boundary  parameter  may  have  a  different  functional  relation¬ 
ship.  Thus,  for  example,  a  particular  motor  may  have  the  following  circumferential 
variation  in  external  film  coefficient. 

Zone  I  0  =>  0O  to  0-^ 

hoj  =  +  bj  0  +  Cj  02 


Zone  II  0 


01  to  02 
a  ■  +  sr  cos  0 

n  *n 


To  permit  the  external  environment  to  be  time  dependent,  the  coefficients 
of  each  boundary  parameter  function,  f(0),  are  considered  to  be  functions  of  time. 

This  is  accomplished  within  the  analysis  and  resulting  computer  program 
by  establishing  a  table  for  each  boundary  parameter  with  ten  time  entries  for  each 
coefficient.  This  is  illustrated  below. 

Time  f(0) 

®  [ho>  *?*>  E,  Tconvl  Trad] 

min,  a  b  c  d  e  f  g 

el  al  bl  C1  dl  el  fl  &1 

&2  b2  c2  d2  e2  %  &2 


01O  a10  b10  c10  d10  e10  f10  2io 
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Intermediate  values  of  the  coefficient.and,  consequently,  the  boundary 
parameters  would  be  obtained  by  linear  interpolation.  Likewise,  separate  tables  are 


permitted  for  each  "zone. " 

Having  defined  the  manner  in  which  the  external  boundary  conditions  are 


simulated,  the  finite  difference  equation  for  the  gradient  at  the  outer  surface  can  be 


written  as  follows: 


SliUllLJs.  =  J!s_  .  T  ,  +?h_  ,T  j4  T  4,  +  5 _ 

Ari  Kj,|t(^onv  T°8'  Kj1k(Tred  T°s)  KJik 

Substituting  this  relationship  into  the  basic  equation  1. 1,  with  Tjt  ^  “ 
Tog*  gives  the  following  expression  for  T'os  in  accordance  with  the  nomenclature 
illustrated  in  Figure  2. 


3.  CASE-PBOPEULANT  INTERFACE  POINT 


3,  Case-Propellant  Interface 


As  can  be  seen  from  Figure  3,  the  case-propellant  interface  nodal  point 
is  influenced  in  its  thermal  response  by  the  characteristics  of  both  the  case  material 
and  the  propellant.  Consequently,  the  finite-difference  expression  for  the  interface 
temperature  can  be.  obtained  by  substituting  into  the  basic  equation  1. 1  for  pCP) 

Kj  ^+1  and  Kj  ^  the  following  relationships: 


.  P2  <cd>2  2rl  Ar2  -  Ar2  *  Pi  <Cp)l  2rl  Arl  +  Arl 

^ *  22 
2rj  Ar2  -  Ar2  +  2rj  Ar^  +  Ar^ 


Kj+i.k+1"  V2  Kl-l.k+l 
Arx  +  Ar2 

Arl  KJ+l,k-l  *  Ar2K]-l,k-l 


(1.9) 

(1.10) 

(1. 11) 


Ar^  + 


Ar. 


These  relationships  permit  the  proper  proportioning  of  each  material 
and  corresponding  thermal  properties  adjacent  to  the  interface.  Upon  making  these 
substitutions  Into  the  general  expression,  the  finite  difference  notation  for  the  interface 
becomes  the  following: 

+  Aafcrj  Ar2  -  Ar22  +  2rj  Arx  +  Ari2] 


Tif  =  TIf 


P2  (Cp)2  [2rJ  “r2  -  ir22]  +  Pi  (Cp)lC2rJarl  +  Arl2] 

’Kl.l.k-KM.fcYTm.k-Tl.k  t  Tj.k-Tj-i.k’ 


Ai*i  +  Ar2 


Ar, 


Arl  (Kj+i.kn-  KJ+l.k-l)  +  Av2  <Kj+l.k-l  '  Ki-1,  k-1 ) 
rj,k2  <Ari  +  Ar2)  ^1  +  A02> 
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In  this  expression,  all  forward  differences  in  r  and  terms  subscripted  with 
a  1  refer  to  the  case  material.  Likewise,  all  backward  differences  in  r  and  terms  sub¬ 
scripted  with  a  2  refer  to  the  propellant.  This  nomenclature  is  in  accordance  with  Figure  3. 

4.  Internal  Surface 

The  internal  surface  of  a  rocket  motor,  being  in  most  cases  irregular 
or  noncircular,  experiences  heat  flow  at  the  boundary  which  is  not  precisely  radial. 
Therefore,  the  boundary  condition  at  the  Internal  surface  of  the  grain  is  expressed  in 
terms  of  the  gradient  normal  to  the  surface,  as  follows: 

(Tc  -  Tjs)  (113) 

where 

fy  =  Internal  Boundary  Film  Coefficient,  BTU/hr-ft^-’R 

Tc  =  Internal  Environmental  Temperature,  BR 

Tis  =  Internal  Surface  Temperature,  BR 
n  =  Direction  Normal  to  Internal  Boundary 


•!<}- 


Since  all  finite  differences  within  the  basic  equation  1, 1  are  expressed  with  respect  to 
the  radial  and  circumferential  directions,  the  normal  derivative  (1, 13)  at  the  internal 


boundary  must  be  resolved  as  shown  in  Figure  4  to  obtain  the  following  radial  and 
circumferential  boundary  gradients, 


3T  6  T 
30  ' 


iLr2  M21 

n  [  \dr/ J 

IT  3T|  1  /drf] 

I  r  3  3  n  [;  2  \d  0/  J 


-1/2 

-1/2 


3T 
3 


.  JU¬ 
KI  J 

r  /, 

a2T1/2 

(Tc  “  Tjt  k) 

1  +  r2  |  df 

L  vdi 

')] 

*'1/2 

(1.14) 

=  Ju¬ 
ki,] 

(Tc  -TJ>k) 

i  +i  ($i 

.  r2W 

III 

(1.15) 

These  derivatives  expressed  in  either  forward  or  backward  finite 
difference  notations,  depending  upon  the  grain  geometry,  must  be  substituted  into  the 
general  heat  conduction  equation  (1. 1)  together  with  other  modifications  necessary  to 
simulate  an  internal  boundary  nodal  point.  In  summary,  the  various  types  of  points 
that  may  exist  and  the  necessary  nodal  substitutions  that  are  required  for  their  simulation 
are  presented  below.  Iterations  are  assumed  to  begin  at  0  =  0  and  r  =>  rQ  progressing 
in  a  counterclockwise  direction  and  with  increasing  r. 


8T  =  jL 

1)  Internal  surface  Implies  8n 


a)  Entering  propellant  on  constant  0: 
Replace  backward  differences  in  r 

•1/2 


Ti.k-Tj-i,k 


■i 


Ar2 


withKj7k(Tj’k'Tc) 


+ rp  y 


Ar2  with  Arx,and  Kj_ljkwith  Kj^1>k 


b)  Leaving  propellant  on  constant  0  implies: 
Replace  forward  differences  in  r 

-rl/2 


Ti  +i  k  -  Ti  k  hi  1  /M  . 

^±k^wllh  yj 


ArxwithAr2,  andKj  +1;kwithKj_i>k 


J 
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c)  Entering  propellant  cn  constant  r  implies: 
Replace  backward  differences  in  0 

2,-l/3 


Ti  I,  -  T<  i.  i  hj 
-iJL  X.wUh-3-  (Tj  k  -  Tc) 

Lh  Kj,k 


1+r 


- 


j  \Ar  J 


A02  withA01>  and  Kj  k.^wlth  Kj  k+i 


d)  Leaving  propellant  on  constant  r  implies: 

Replace  forward  differences  in  0 

r  /a  \2i  "l/2 

Tl’'~WltbK~  <T°  'Tj>k)  1  +rH5r)  ’  i0iwlthA02'  andKJ,knwlthKj,k- 
^1  j}»^  •  \  /  * 

1  w 

2)  $  =  0  implies  «  0  and  consequently  the  following  substitutions: 
Repiace  Tj>  k-1  with  Tjf  k+1 

Kj  k-i  with  Kj>k+1 
and  A02  with  A0j 

X  rn 

3)  0  =  0max  l3nP^es  =  0  and  consequently  the  following  substitutions: 
Replace  TJiknw!th  Tj  ^ 

KJ,k+lwlthKJ1k-l 
and  A0^  with  A0g 

4)  Propellant-case  interface  implies  the  following  substitutions: 


Replace  puwith 


KJ,k+l  wlth 


P2  =2  f.2rl  Ar2  *  Ar22J+  "l  °1  [2rJ  Arl  +  Arl2] 

2  2 
2rj  Ar2  -  +  2r j  Arj  +  Ar^ 

Ari  Ki+i,kn  *Ar2  Ki-i,k+i 


Ari 


Arc 


and 


Kj,  k-1  with 


Arl  Kj+1,  k-1  +  Ar2  Kj-1  k-1 


Ar^  + 


Arr. 

A 


-19- 


5)  External  surface  implies  the  following  substitutions: 

ReplaCe  ^  Art~  (Tconv  '  TJ-k>  * 

<TracTTj,k4)  +K^"k*  Arl  w*th  and  Kj+l,kwith  Kj.1|k 

6)  Internal  point  implies  no  change  in  basic  equation, 

The  above  substitutions  may  be  performed  separately  or  in  combinations 
as  required  to  describe  completely  the  nodal  point  condition.  For  the  internal  surface, 
the  values  of  A0  and  Ar  used  in  the  directional  derivative  are  those  which  best  describe 
the  boundary  geometry.  How  well  the  geometry  is  simulated  depends  upon  the  fineness 
of  the  grid  and  the  proper  selection  of  A0  and  Ar. 

Typical  nodal  points  are  illustrated  on  a  hypothetical  grain 
configuration  given  in  Figure  5  and  identified  as  follows: 

Point  No.  Identification 

1  0  =  0,  Internal  surface,  and,  entering  propellant 

on  constant  0,  adiabatic  along  constant  r 

2  0  =  0,  Adiabatic  along  constant  r 

3  0  =  0,  Propellant-case  Interface,  adiabatic  along 

constant  r 

4  0  =  0,  External  surface,  pdiabatic  along  constant  r 

5  Internal  surface,  entering  propellant  on 

constant  0 

6  Internal  point 

7  Propellant-case  interface 

8  External  surface 
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A-6556 


FIGURE  5.  SEGMENT  OF  GRAIN  WITH  PARTIAL  GRID  ILLUSTRATING  TYPICAL 

NODAL  POINTS 


Point  No. 


Identification 


9 

Internal  surface,  entering  propellant  on 
constant  0,  leaving  propellant  on  constant  r 

10 

Internal  surface,  leaving  propellant  on 
constant  r 

11 

Internal  surface,  entering  propellant  on 
constant  ^entering  propellant  on  constant  0 

12 

Internal  surface,  entering  propellant  on 
constant  r 

13 

0  -  0max,  internal  surface,  entering  propellant 
on  constant  0 

14 

0-0 

max 

15 

0  -  0max,  propellant  case -interface 

16 

0  -  0  ,  external  surface 

SiaJv 

Each  nodal  point  is  identified  within  the  program  by  a  code  which  directs  >. 
the  substitutions  and  modifications  of  the  basic  equation. 

5,  Description  of  Computer  Program 

The  method  outlined  in  the  previous  sections  for  the  modification  of  the 
basic  heat  conduction  equation  to  obtain  the  finite  difference  equations  representing 
the  temperatures  at  various  nodal  points  throughout  a  propellant  grain  has  been 
programmed  for  the  IBM-7070  digital  computer  in  the  FORTRAN  II  language. 

Basically,  there  are  two  separate  programs.  The  main  program  is 
capable  of  considering  a  grid  configuration  consisting  of  30  nodal  points  on  any  radial 
lino  and  as  many  circumferential  points  (radial  lines)  as  are  required  to  describe  the 
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grain  completely.  Each  nodal  point  must  be  formed  by  the  intersection  of  a  radial  and 
a  circumferential  line.  Consequently,  whenever  a  radial  lin^  passes  through  a  boundary 
line,  a  circumferential  line  must  also  be  constructed  through  the  same  point  (see 
Figure  5), and  conversely. 

In  the  iterative  process,  temperature  computation  begins  at  the  internal 
surface  on  the  first  radial  line  and  progresses  in  the  direction  of  increasing  r  until  the 
case  surface  is  reached.  The  value  of  0  is  then  incremented,  and  computation  begins 
again  at  the  inner  surface.  This  process  is  repeated  until  the  entire  configuration  is 
scanned.  At  this  point,  time  is  incremented,  the  program  is  reoriented  at  the  initial 
starting  point,  and  computations  proceed  as  before  using  the  newly  computed 
temperatures.  This  is  a  "marching'*  procedure  where  the  program  advances  in  time 
computing  new  temperatures  based  upon  the  previous  values.  By  the  nature  of  the 
finite  difference  equations,  tin  program  requires  three  complete  radial  lines  and  the 
corresponding  nodal  data  in  memory  at  all  times.  Therefore,  in  order  to  have  the 
capability  of  considering  an  unlimited  number  of  radial  lines,  the  use  of  magnetic  tape 
for  the  storage  of  intermediate  temperatures  and  data  was  required.  This,  while 
necessary,  resulted  in  extremely  long  computation  times. 

A3  a  result,  the  program  was  revised,  whereby  all  the  nodal  temperature 
data  were  retained  in  memory  together  with  the  newly  computed  and  the  previous 
temperatures  for  each  node.  This  was  made  possible  by  limiting  the  number  of 
circumferential  points  (rad'al  lines)  to  27  while  retaining  the  original  number  of  30 
radial  points  (circumferential  lines).  The  effect  was  a  significant  reduction  in 
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running  time  and  input  cards  required.  The  latter  was  reduced  to  approximately  one- 
half  while  the  former  was  reduced  by  approximately  one-third. 

The  glossary,  flow  chart,  subroutine  identification  and  input  formats  for 
both  programs  are  contained  in  Appendix  I  to  this  report.  ' 

T3,  Test  Cases 

In  order  to  evaluate  the  performance  of  the  program  and  to  demonstrate 
the  utility  and  versatility  of  the  analysis,  numerous  preliminary  test  runs  were 
performed.  There  was  sufficient  variation  among  these  test  cases  so  that  each 
particular  routine  within  the  program  was  used.  Each  test  case  Is  described  below 
and  illustrated  in  Figures  6  through  9,  respectively. 

1)  A  center  perforated  grain  with  a  thin  steel  case  wall.  The  »•*. 
center  perforation  was  extremely  small  in  order  to  approxi¬ 
mate  an  end  burning  grain.  The  external  boundary  conditions 
were  as  follows: 

q  =  E  cos  0  =  400  cos  0  BTU/hr-ft2;  0  -  0*-90® 

q  “  ?0  <Trad4  ’  Tos4)  -  0. 3<r  (1004  -  Tos4)  BTU/hr-ft2;  0  =  90o-180* 
The  internal  boundary  condition  was  adiabatic. 

2)  Same  as  (1)  without  the  case  wall. 

3)  Same  as  (2)  with  the  following  external  boundary  conditions: 
q  -  E  cos  0  =  120  cos  0  BTU/hr-ft2;  0  =  0°-90° 

q  -  <7  <Trad4  -  Tos4)  0. 3(7  (1004  -  T0£,4)  BTU/hr-ft2;  0  =  90°-180° 

4)  A  30°  segment  of  a  star  grain  with  no  case  wall.  The  external 
boundary  condition  was: 

q  ;  h0  (TConv  -  Tos)  --  8.5  <580-TOS)  BTU/hr-ft2,  0  -  0u-30° 
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Propellant  Properties 

No  Case  Wall  Thermal  Conductivity  =  0. 14  BlTT/hr-ft-Tt 
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12  o  125  in.  Rad  _ J  hoaxer  Perforation  1.1251a. 


q-  0.26(580-Tia) 


The  internal  boundary  condition  was: 
q  =  ^  (Tconv  ”  Tos)  =  0. 25  <580-TOs)  BTU/hr-ft2 

Test  case  1  evaluated  the  external  surface  and  case-propellant  interface 
nodal  point  relationships.  A  significant  result  obtained  from  this  case  was  the  severe 
stability  requirements  imposed  on  the  solution  by  a  thin  case  wall  having  a  relatively 
higher  thermal  diffusivity  than  the  propellant.  These  stability  requirements  resulted 
in  small  iteration  times,  and,  consequently,  long  computation  times.  Since  the  validity 
of  the  case  wall  and  case-propellant  temperatures  was  demonstrated  within  the  program, 
the  remaining  test  cases  were  evuluated  without  a  case  material.  The  capability 
exists  to  consider  a  case  wall,  but  unless  it  contains  sufficient  thermal  resistance  or 
capacity  to  grossly  affect  the  gradients  within  the  grain  it  is  rot  analytically  necessary 
or  economically  feasible  to  do  so, 

The  results  of  test  cases  2  and  3  are  presented  in  Figures  10  and  11 
corresponding  to  12  hours  of  expoeure.  Intermediate  profiles  corresponding  to 
every  0. 5  hours  were  also  obtained  for  these  configurations. 

Test  case  4  illustrates  the  effect  of  a  small  thermally  symmetric  star 
segment  upon  the  temperature  gradient  through  a  grain  subjected  to  uniform  internal 
and  external  convection.  The  results  of  this  test  case,  presented  in  Figure  12  in  the 
form  of  an  isotherm  plot,  show  the  heat  sink  effect  obtained  within  the  grain  as  heat 
enters  from  both  boundaries.  At  the  center  of  this  sink  is  the  point  within  the  grain 
which  responds  last  and  slowest  to  the  thermal  environment.  Consequently,  this  is  a 
region  of  low  thermal  gradients. 
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To  demonstrate  more  extensively  the  effects  of  a  nonuniform  external 
environment  and  an  irregular  internal  boundary  upon  the  thermal  profile  within  a 
propellant  grain,  the  180*  star  configuration  illustrated  in  Figure  13  was  evaluated  for 
the  following  three  separate  sets  of  boundary  conditions.  The  initial  temperature  of 
the  grain  was  considered  to  be  54G8R. 

1)  External  Boundary  Condition 

0  -  0*  to  90* 

q  =  400  cos  0  BTU/hr-ft2 
0  =  90*  to  180* 

q  -  0. 3  a  (1004  -  Tqs4)  BTU/hr-ft2 
Internal  Boundary  Condition 
q  =  0  (adiabatic) 

2)  External  Boundary  Condition 

Same  as  (1)  above 
Internal  Boundary  Condition 

a  =  h  ^Tc  -  Tis)  =  5(650  -  Tjs)  (heated) 

3)  External  Boundary  Condition 

Same  as  (1)  above 
Internal  Boundary  Condition 

q  =  h(Tc  -  T«s)  =  5(300  -  Tis)  (cooled) 
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The  isotherms  existing  within  the  grains  for  the  above  three  environments 
are  illustrated  in  Figures  14  through  16.  These  demonstrato  the  effect  of  the  irregular 
internal  boundary  together  with  the  response  of  the  grain  to  either  heating  or  cooling 
from  the  inside.  The  resolution  of  the  solution  is  sufficient  to  outline  clearly  the 
"sink"  and  "source"  effects,  respectively,  for  both  internal  environments. 

While  these  test  cases  are  fictitious,  they  do  illustrate  the  capability  of 
the  program  to  predict  temperature  profiles  through  a  grain  when  subjected  to  adverse 
environments  for  subsequent  evaluation  of  the  effects  upon  ballistic  performance  and/or 
structural  integrity. 

B.  ADVANCED  GRAIN  DESIGN  ANALYSIS  -  PROGRAM  AGDA 

The  Thiokol-Elkton  Advanced  Grain  Design  Analysis  (AGDA)  computer 
program  represents  the  most  general  truly  three-dimensional  grain  design  computer 
program  available  today.  The  present  computer  program  provides  the  rocket  designer 
with  all  fundamental  geometrically  derivable  rocket  motor  performance  parameters. 
These  parameters  including  burning  surface  area,  chamber  volume,  propellant  weight, 
center  of  gravity  travel,  and  moments  and  products  of  inertia.  With  these  fundamental 
parameters,  interpretation  in  terms  of  ballistic  performance  tradeoffs  is  immediately 
and  easily  available. 

The  approach  utilized  in  program  AGDA  represents  an  entirely  new  concept 
in  the  analysis  of  solid  propellant  rocket  motor  grain  design.  This  new  concept,  which 
will  be  described  later,  yields  results  of  extremely  high  accuracy  through  use  of 
vector  algebra  and  calculus  techniques. 
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.  ISOTHERM  PLOT  FOR  180*  STAR  CONFIGURATION 


CONFIGURATION  (COOLED  INTERNAL  BOUNDARY) 


The  entire  system,  both  analysis  and  programming,  is  completely  modular 
in  design,  thus  readily  permitting  the  inclusion  of  additional  capabilities  and  techniques. 
Figures  17  and  18  are  representative  of  the  wide  class  of  solid  propellant  rocket  motors 
which  may  be  handled  today  by  Program  AGDA. 

It  Description  of  Program  AGDA 

Because  of  the  uniqueness  of  the  Elkton  Advanced  Grain  Design  Analysis, 
a  brief  description  of  the  theory  of  the  program  is  presented  herein. 

The  three  most  salient  characteristics  of  the  Advanced  Grain  Design 
Analysis  are  the  conceptual  simplicity  of  the  Elkton-Moore  method  (upon  which  it  is 
based),  its  generality,  and  the  accuracy  it  affords.  In  addition  to  providing  chamber 
volume  and  burning  surface  area  as  functions  of  burning  distance,  this  method  permits 
calculation  of  propellant  mass,  the  components  of  the  inertia  tensor,  and  propellant 
center  of  gravity. 

Briefly,  the  Elkton  Advanced  Grain  Design  Analysis  can  be  summarized 
in  the  following  eight  steps: 

1)  The  chamber  and  case  of  the  rocket  motor  are  represented 
as  mathematical  surfaces. 

2)  The  surface  representing  the  outside  of  the  chamber  is 
divided  into  a  large  number  of  pseudo-rectangular  elements. 

3)  Normals  to  this  surface  are  constructed  at  the  corners  of 
a  "rectangle"  such  that  they  travel  through  the  propellant 
and  intersect  the  case. 

4)  If  the  angle  between  any  pair  of  normals  emanating  from 
the  same  point  as  calculated  in  Step  3  exceeds  a  given 
small  value,  new  normals  are  constructed  between  those 
already  calculated  in  such  a  way  as  to  reduce  the  angular 
separation  between  successive  normals  to  a  tolerable  value. 
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5)  Each  of  the  normals  constructed  In  Steps  3  and  4  is  divided 
into  a  number  of  equal  increments  beginning  at  the  chamber 
surface, 

6)  Surface  area  between  each  set  of  four  points  of  corresponding 
burning  distance  i*  calculated.  Volume  and  momenta  and 
products  of  inertia  are  calculated  for  regions  of  propellant 
between  successive  sets  of  points, 

7)  Steps  5  and  6  are  repeated  for  all  sets  of  four  normals 
calculated  in  Step  4;  then,  Steps  3,  4,  and  5  are  repeated 
for  each  of  the  pseudo-rectangular  areas.  On  each  pass 
through  these  "loops, "  values  of  surface  area,  volume, 
etc.  corresponding  to  the  same  burning  distance  are  added 
together. 

8)  The  quantities  summed  in  Step  8  are  converted  to  the 
following  motor  parameters: 

a)  Burning  Surface  Area 

b)  Chamber  Volume 

c)  Unconsumed  Propellant  Mass 

d)  Moments  of  Inertia  of  Unconsumed  Propellant 

e)  Products  of  Inertia  of  Unconsumed  Propellant 

f)  Center  of  Gravity  of  Unconsumed  Propellant 

These  steps  are  shown  in  a  generalized  flow  chart  of  the  program  in 


Figure  19. 


Theoretical  aspects  and  derivation  of  the  mathematics  behind  the 
Advanced  Grain  Design  Analysis  are  presented  in  the  next  section  of  the  report.  Followir 
this  derivation,  the  method  of  application  to  a  generalized  solid  propellant  rocket  motor 
is  presented.  A  knowledge  of  the  eight  steps  outlined  above  i3  all  that  is  necessary 
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to  ensure  comprehension  of  the  majority  of  the  section  concerning  application,  A 
study  of  the  deviation,  however,  will  provide  more  insight  into  the  methou. 

2.  Theoretical  Aspects  of  the  Advanced  Grain  Design  Analysis  -  Uniform 
Temperature 

a.  Representation  of  Grain  and  Definition  of  Surface  Normals 

Given  a  general  grain  design  with  inner  and  outer  web  limiting  surfaces 
represented  by  2!  (x,  y,  z)  *-  0  and  y  (x,  y,  z)  M  0  respectively  (see  Figure  20),  an 
x-y  mesh  is  constructed  on  that  portion  of  the  plane  z  =  Zo  intercepted  by  the 
projection  of  0  (x;  y,  z)  3  0. 

The  regression  of  the  total  surface  is  obtained  by  considering  the 
regression  of  each  of  the  surface  elements  of  0  (x,  y,  z)  =  0  corresponding  to  a  block 
of  the  x-y  mesh.  Since  burning  proceeds  in  a  direction  normal  to  the  instantaneous 
burning  surface,  the  regression  of  any  surface  element  is  bounded  by  the  normals  to 
that  surface  evaluated  around  the  perimeter  of  the  element  or  approximately  by  the 
normals  at  the  corners  of  the  element.  To  Illustrate  this,  only  a  single  representative 
element,  such  as  the  one  shown  in  Figure  21,  will  be  considered.  With  respect  to 
this  figure,  the  unit  normals,  Nj,  are  defined  by: 

yv0(x,  y,  z) 
i  =  1,2, 3, 4 

where  Pj  (0)  represents  the  corners  of  any  rectangular  clement. 


Pi  (0) 


(1.  li 


A-3633 


b. 


Control  of  Accuracy 


•  * 


Since  the  accuracy  of  the  entire  calculation  is  directly  related  to  the 
angular  sepitaticn  of  the  hormals,  si  criterion  must  be  established  for  their  selection. 


if  8q  is  the  maximum  permissible  angular  separation  between  any  two  adjacent 


normals,  reducing  60  reduces  the  error  until,  as  0Q  approaches  aero,  the  error 


approaches  zero.  Thus  there  exists  a  value  of  90  which  will  ensure  any  predetermined 
degree  of  accuracy. 

.  *'  '  i  !  .  *..•*'  • 

The  angular  separation  of  consecutive  normals  is  considered  under 
two  separate  conditions,  each  of  which  is  handled  in  a  completely  independent  manner. 

The  first  of  these  to  be  considered  is  the  case  of  normals  emanating 
from  the  same  point. 

Initially  the  surface  curvature  is  assumed  to  be  minimal  and  normals 
are  constructed  for  the  coordinates  as  dictated  by  the  specified  mesh  interval  size. 

The  angle  between  these  normals  is  then  compared  with  the  maximum  permissible 
angle.  If  the  latter  has  been  exceeded,  a  new  coordinate  is  determined  using  the 
following  variable  mesh  relationship: 
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(1.17) 


where: 


V  (R,  2)  31  equation  of  the  core  head 

dfl  =»  maximum  allowable  angle  between  normals 

dR  =  maximum  allowable  mesh  increment,  because 

of  dd 

The  other  condition  Is  that  of  normals  emanating  from  the  same  point.  Stated 

mathematically,  the  problem  is  to  find  a  vector,  N,  in  the  plane  of  two  other  vectors 

^ 

(say  Nj  and  N2),  a  fixed  angle  9  away  from  Nj.  One  should  be  able  to  express  N  as 
follows: 

N  »  Ll^l+L2^2  (1.18) 

From  the  theory  of  linear  algebra,  this  Is  always  possible  when 
and  N2  are  co-planar  but  not  co-linear. 

■jb  *4* 

Proof  that  the  vector  N  is  in  the  plane  of  Nj  and  N2  may  be 
demonstrated  by  calculation  of  the  triple  vector  product  of  1^,  Njj,  and  N.  Since  this 
triple  product  represents  the  volume  enclosed  by  a  parallelopiped  of  sides  parallel 
to  the  given  vectors,  a  calculated  value  of  zero  (volume)  shows  that  all  three  vectors 
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Application  of  the  law  of  cosines  shows  that  L*  and  Lg  are  given  by: 


"  .  sin  9 

=  oos  9  -  Lg  cos  iY, 
where  9  =  angle  between  I?!  and  N 

.t  =»  angle  between  Ni  and  Ng. 

The  derivation  holds  for  0<Y  <  rr»  provided  that  9  f . 

Thus,  by  ues  of  the  angular  control  variable  90,  one  may  directly 
control  the  accuracy  of  the  resultant  motor  evaluation  for  both  situations.  By  choosing 
a  large  value  of  90  one  obtains  relatively  inexpensive  answers  suitable  for  preliminary 
design  purposes;  by  choosing  a  small  value  of  90,  answers  of  extremely  high  accuracy 
may  be  obtained. 


c.  Maximum  Burning  Distance 

Since  burning  along  the  i^  normal  proceeds  only  to  the  outer  web 
limiting  surface,  a  maximum  burning  distance,  d  maxj,  must  be  calculated.  This  is 
accomplished  by  the  simultaneous  solution  of  a  scalar  equation  of  the  normal  Ni  and 
the  equation  of  the  case  f  (x,y,  z)  =  0. 

The  coordinates  x*,  y^  and  Zj,  of  a  point  (d),  a  distance  d  along 
the  normal  from  the  inner  web  limiting  surface,  0  (x,y,  z)  =  0,  are: 

X1  =  xi*  +  d  cos  Qfj 

y*  “  yf  +  d  cos  (1.20) 

zi  =  zl°  +  d  COH  y  i 
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Note  thut  the  only  variable  In  tho  elemental  surface  area  equation 

(1.25)  Is  the  burning  diatanoe,  d;  tho  vectors  C1#  C2,  and  <$3  have  been  determined 

by  the  inner  web  limiting  surface. 

The  volume  of  propellant  consumed  as  a  result  of  the  surface 

k-1  k 

element  burning  out  tho  distance  from  "  2  AJ  to  d^  *  X  Aj  is  given  by  the 

j-1  J»1 

approximation  (Figure  22); 

Sk  +  Sk-i 

vk  “  - 2 -  Ak  <1,27> 

To  retain  over-all  program  accuracy  especially  with  respect  to  the 
volume  dependent  parameters,  *  an  extremely  accurate  method  of  evaluation  of  the 
propellant  volumes  adjacent  to  the  case  was  needed.  This  requirement  Is  attributable 
to  the  random  manner  In  which  each  of  the  four  normals  may  intersect  the  case  and 
the  significant  variations  in  the  respective  maximum  burning  distances  which  result. 

The  technique  employed  in  the  program  approximates  the  exact 
three-dimensional  solution  In  the  following  manner:  given  dmaXi,  dmaX2>  d^ax^ 
d*nax4  aro  tke  maxlmum  burning  distances  along  each  of  the  normals  N£, 

N3,  and  N4,  respectively,  the  two  larger  and  two  smaller  values  aro  independently 
averaged  to  obtain  DmaXL  (larger  maximum  burning  distance)  und  Dmaxg  (smaller 
maximum  burning  distance).  Using  these  parameters  the  various  shapc*3  of  remaining 
propellant  volumes  may  be  extremely  well  approximated  by  use  of  the  trapezoidal  rule. 
Pictorial  representation  of  these  shapes,  tho  applicable  logic  conditions  which  exist, 
ar.d  the  respective  equations  for  surface  and  volume  evaluation  aro  shown  in  Figure  23, 


♦Propellant  volume,  mass,  and  moments  and  products  of  Inertia 


^  3<jp  JOJ  (p*  s  =  (Vs 


e.  Elemental  Momenta  and  Products  of  Inortla 


Since  the  Increments  A,  Ax,  Ay,  and  Az  aa  well  as  the  maximum 

angle  between  normals,  0O,  may  be  arbitrarily  assigned  any  value,  sufficiently  small 

values  can  be  chosen  so  that  the  mass  of  any  block  can  be  treated  as  if  it  were 

concentrated  as  its  cantor  of  mass. 

If  the  eight  corners  of  the  volume  element  under  consideration 

have  coordinates  <xj,  ,  zj),  J  al,  2. . .  8,  the  coordinates  of  the  center  of  mass 

Xk,  are  approximated  by: 

8  8  8 

Xksi/8  2  xj,  Yk=:i/8  2  y,  Z^l/8  2  z,  <1.28) 

j-1  J-l  '  J»1 

The  moments  and  products  of  inertia  of  the  shell  are: 

<Wk  3  <\2  +  V)  Mk  <ixy>k  3  -W  Mk 

fyylk  3  $k2  +  ^k2)  wk  ^xy)k  3  "^k^k^k  (1.29) 

<!zz)k  3  (Xk2  +  ^k2)  Mk  (]yz)k  3  ~^k^k  Mk 

Because  of  the  cylindrical  symmetry  in  most  practical  rocket 

systems,  the  only  first  moment  to  be  considered  Is: 


Iz  =  Zk  Mk  (1  30) 

This  will  be  used  to  calculate  the  Z  coordinate  of  the  center  of  mass. 

f.  Total  Surface,  Moments  and  Products,  Mass,  and  the  Center  of  Mass 
To  review  briefly,  we  have  calculated  an  elemental  burning  surface 
area  Sk,  an  elemental  mass  Mk,  moments  and  products  of  inertia  (I^k,  (lyy)k, 


(izz)k,  (Ixy)k,  (lxy)k,  (Iyz)k»antl  flrsl  moment  (1^  at  each  value  of  burning 

k 

A 


distance,  d  ~  A^.  A^  f  A2  Aj  +  A2  <•  A^ . .  £  At  =  D( 


J  =  l 


max, 
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tho  avorugo  maximum  distance  botwoun  the  inner  and  outer  web  limiting  surfaces 
for  this  olomont. 


As  the  analysis  continues  from  block  to  blook  of  tho  x-y  mesh, 
these  quantities  are  summed  for  corresponding  values  of  k.  Thus,  after  all  the  blocks 
of  the  mesh  have  been  considered,  the  area,  momenta,  and  mass  for  the  "sheila" 
within  the  web,  as  shown  in  Figure  24,  are  derived. 

Finally,  since  only  the  moments  of  inertia,  products  of  inortia,  and 
mass  of  the  unconsumed  portion  of  propellant  are  of  interest,  tho  shells  must  be  added 
in  the  following  manner: 


K, 


max 


^xx/J  =  ^  Oxx)k 


k  3  J 


(1.31) 


where  “  xx  moment  of  inertia  of  tho  unconsumed  propellant  after  burning 

2  Ai  units; 

1=1 

^xx^k  *  **  moment  of  inertia  of  the  kth  shell 

kjnax  =  tota^  number  of  shell  required  to  evaluate  the  motor. 

Since  chamber  volume  (Vc)j  is  defined  as  the  volume  of  the  case 

minus  the  volume  of  propellant,  the  equation  becomes: 

M 

(Vc)j  -  v  - 

v'here  V  =  volume  of  ease 

p  3  density  of  propellant 

Mj  -  mass  of  remaining  propellant. 


{1.  32) 
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(1.33) 


Th®  Z  coordinate  <Z0)j  of  the. center  of  mass  after  burning  a  distanoe  dj  la  given  byi 

<!.>j 


where  (I  ),  and  It  are  defined  as  in  equation,  1. 31. 

Using  thir  method,  moment#*  and  prod  jots  of  inertia,  total  mass  of 
propellant,  center  of  mass,  and  burning  sur  face  area  of  a  solid  propellant  rocket  grain 
for  an  arbitrary  number  of  burning  distances  can  bo  calculated.  The  general  analysis' 
as  presented  bore  is  applicable  for  calculations  of  the  head-end,  aft-end,  and 

straight-through  portions  of  virtually  any  solid  propellant  grain  configuration.  The 

:i  .  <  •  .  ■  •  •  .  • 

application  of  this  analysis,  via  digital  techniques,  is  purely  restrioted  to  the 
ingonuity  cl  the  programmer,  the  sophistication  of  the  electronic  data  processing 
equipment  available,  and  the  order  of  priority  of  grain  designs  to  be  considered. 

3.  Theoretical  Aapeota  of  Nonuniforra  Tomparature  Qrain  Design  Anaiyels. 
a.  Preliminary  Considerations 

One  of  the  more  salient  characteristics  of  Program  AGDA  is  the 
fact  that  the  propellant  is  the  rocket  motor  under  consideration  is  divided  into 
literally  thousands  of  small  elements  during  the  analysis.  Since  the  coordinates  of 
each  small  element  are  known  in  addition  to  its  volume  and  surface  area,  a  "handle" 
is  readily  available  for  the  consideration  of  geometrically  dependent  propellant  properties. 
The  most  important  o?  these  properties  is  grain  temperature,  . 

The  effect  of  initial  grain  temperature  on  solid  propellant  rocket 
performance  is  well  understood,  provided  that  the  entire  "rain  is  at  tho  same 
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temperature,  In  general  a  warmer  motor  burns  for  a  shorter  time  and  at  a  higher 
chamber  pressure  than  an  identical  motor  at  a  lower  grain  temperature.  The  total 
impulso,  to  a  first  approximation,  is  independent  of  the  initial  grain  temperature,  at 
is  the  progrossivity  or  regressivlty  characteristics  of  the  pressure-time  or  thrust¬ 
time  trace. 


When  a  nonuniform  temperature  gradient  exists  within  the  grain  of 
a  solid  propellant  rocket  motor,  no  simple  method  exists  to  predict  motor  performance. 
The  slope  and  shape  of  the  pressure-time  and  thrust-time  ourves  may  be  altered 
appreciably  from  that  of  the  corresponding  curves  obtained  at  a  uniform  grain 
temperature.  The  principal  purpose  of  this  study  area  has  been  to  examine  in  detail 
the  effect  of  a  nonuniform  temperature  gradient  within  the  grain  of  a  solid  propellant 
rocket  motor  through  a  generalization  of  the  Thiokol  -Elkton  Advanced  Grain  Design 
Analysis  computer  program. 

b.  Derivation  of  Technique  of  Nonuniform  Temperature  Grain  Design 
Analysis 

The  burning  rate,  r,  of  a  solid  propellant  rocket  grain  at 
temperature  Tg  is 


rce 


1TkfV  VfO?c>  (1-34) 

In  this  equation,  r0  Is  the  burning  rate  at  grain  temperature  Tg0  and  chamber  pressure 
^CO1  andrj,  is  the  temperature  sensitivity  coefficient. The  function f  (Pc)  Is  usually 


of  the  form  (Pc/PCo)  • 


dx 


If  dx  is  an  increment  of  burn  distance  such  that  r  =  —  and 

at 


-58- 


rQ  *  — 2 ,  equation  1. 34  may  be  wrltton  as 
dt 


dx  *  dx0o 


'kP't  " 


<1. 35; 


The  value  of  "dx"  ia  required  to  perform  a  surface  regression 
analysis  which,  in  turn,  is  required  to  compute.’'P0. "  However,  P0  was  required 
Initially  to  compute  dx.  In  order  to  avoid  this  difficulty,  the  surface  regression  analysis 
may  be  performed  at  the  effective  chamber  pressure  Pco,  Since  f  (P0)  ■  1  when 


P0  *  PCQ,  the  equation  for  dx  becomes 
dx  ■  dxjje  ( 

After  the  entire  surface  regression  analysis  has  been  completed 

using  dx’s  from  equation  U  36,  the  equilibrium  chamber  pressure  may  be  defined  as 

P na  ^2  1-n 
At 

Since  r  =•  ^or  dt3  formulas  (L.35),  (1.36)  and  (1.37)  may  be  used  to  compute 


(1,36) 


(1.37) 


dx/k^g-Tgo1  f(pco)  =  ^ 

nk  “  Tgo)  -  '  rof 

r„e  K  6  *  f(Pe) 


r0f<pc) 


In  application  of  equation  (1. 38),  dxQ  is  chosen  to  be  a  fixed 
distance,  say  5%  of  the  web  thickness;  rQ  is  determined  by  experimentation  with  the 
propellant  under  consideration.  Finally,  if  f  (Pc)  has  its  usual  form,  the  equation 
for  dt  may  be  expressed  as 


(1.38) 


dx„  P, 


o  *co 


(1.39) 


Figure  25  is  intended  to  Illustrate  much  of  the  foregoing  theory. 
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Direction  of 
Burning 


Direction  of 
Burning 


Uniform  Grain  Temperature 

dx0  3  dx^  =  d.^2  3  dx3 
when  Tg0  »  Tgl  *  Tg2  ■  Tg3 


Nonunlform  Grain  Temperature 


cbtj^  <  tbc2  <  dx3 
when  T„,  <  Tf.o  <  Tg:* 


Lgl 


ann  uus  Ufcaj  < 


f 


Actuals  stance  between  burning 
surface  calculations 


=  dx  =  dxQe  ^  ^ 


Time  corresponding  to  above 
3urn  distance 


=  dt  = 


C'-XC 


r0  f (Pc) 


FIGURE  '*5.  COMPARISON  OF  ELEMENTAL  SURFACE  REGRESSION  BETWEEN 
UNIFORM  AND  NONUNIFCRM  TEMPERATURE  GRAINS 
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°,  Incorporation  of  Variable  Grain  Temperature  Effect  Within  the 

Framework  of  Program  AGDA 

The  actual  incorporation  of  the  foregoing  theory  within  the  frame¬ 
work  of  Program  AGDA  may  best  be  understood  through  a  brief  review  of  the  program 
as  previously  described,  A  simplified  Program  AGDA  for  the  uniform  temperature 
condition  is  shown  in  Figure  26, 

la  order  to  proceed  to  the  variable  temperature  case,  only  two 
modifications  to  this  logic  diagram  are  required.  The  first,  modification  simply  provides 
the  compute  ~  memory  with  a  table  of  grain  temperatures  for  various  coordinates.  Each 
temperature  is  immediately  translated  to  the  corresponding  value  of  A  (incremental 
burning  distance)  by  the  formula 


A  =  AQe 

rk  (Tg  -  Tg0) 

i 

where  Aq 

=  5  percent  of  the  web  distance 

=  temperature  sensitivity  coefficient 

Tg 

=  grain  temperature  at  this  point 

TgO 

=  reference  g^ain  temperature 

The  second  modification  involves  two  steps: 


1)  The  coordinates  of  the  center  of  the  block  of 
propellant  under  consideration  are  computed. 

2)  These  coordinates  are  used  to  find  the 
appropriate  value  of  A  to  use  for  this  particular 
block  by  the  formula 


(1.40) 


A  =  A 


(i,J)  + 


X  ••  X| 

*1+1” xi 


(Amj  -  AIf)  + 


£  -  yj 


yjM-yj 


(Ai  ju-i- Ai(  j)  (1.41) 
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where 

x,  y  are  coordinates  of  the  center  of  the  block 

xj,  y,  are  coordinates  for  which  vtluee  of  A  are  known 
i  exactly  from  the  table 

Aj  jta  value  of  A  associated  with  (xj,  yj). 

These  modifications  are  shown  in  proper  perspective  in  the  logic 
diagram  of  Frograpj  AGDA  with  variable  temperature*  which  is  shown  in  Figure  27,  An 
example  illustrating  the  importance  of  Inclusion  of  the  variable  grain  temperature 
capability  within  Program  AGDA  is  given  in  Section  4  (d)  of  this  report. 

The  approach  derived  is  not  only  unique  but  extremely  desirable  in 
that  it  manages  to  maintain  both  the  simpl  icity  of  the  Advanced  Grain  Design  Program 
and  the  separation  of  the  program  from  the  heat  transfer  analysis.  Since  the  latter 
is  an  extremely  complex  analysis,  this  separation  is  quite  desirable.  Further,  the 
appr  jaah  undertaken  can  be  extended  to  consider  effects,  similar  to  temperature,  which 
cause  a  positional  variation  of  turning  rate.  In  principle,  the  only  requirement 
necessary  is  knowledge  of  the  positional  variation  and  the  associated  effect  upon 
ballistic  characteristics. 

4.  Application  of  the  Theoretical  Aspects  of  Program  AGDA 
a.  Mathematical  Representation  of  the  Rocket  Motor 

Before  an  outline  of  the  application  of  Program  AGDA  to  a  real 
problem  is  given,  it  ia  desirable  to  outline  the  philosphy  by  which  the  present  program 
was  devised. 
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FIGURE  27,  riiOGIlAM  AGDA  WITH  VARIABLE  T  E  M  P  E II  AT  UR  3 
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The  first  step  in  applying  the  Advanced  Grain  Design  Analysis  used 
above  is  to  select  a  generalized  mathematical  model  of  a  rocket  motor,  II  this  model 
is  too  simple,  the  resulting  program  will  be  too  specific  and  have  limited  applicability. 

If  the  model  is  too  complex,  however,  program  utilization  will  suffer  because  of 
the  quantity  of  input  data  required.  The  model  selected  according  to  this  criterion 
is  described  below. 

The  case  is  divided  into  three  sections;  head  end,  aft  end,  and  strtight- 
through  portion.  Both  the  head-  and  aft-end  equations  must  be  either  concave  or  convex 
ellipses  of  revolution.  When  one  considers  the  totality  of  possible  ellipses,  thin 
requirement  is  reasonably  general. 

The  straight-through  portion  of  the  case  must  be  cylindrical; 
however,  both  the  length  and  radius  are  arbitrary.  Note  that  the  arbitrary  length 
feature  allows  the  length  of  this  section  to  equal  zero  when  spherical  ar.d  elliptical 
motors  ara  considered. 

The  surface  of  the  chamber  is  divided  into  a  head  end  and  a 
straight-through  portion.  As  shown  in  figure  2R,  the  head-end  inner  surface  may  be 
divided  into  as  many  as  five  general  conic  sections  of  revolution,  *  thus  allowing  for 
any  configuration  from  a  perfectly  flat  or  a  spherical  head  to  a  head  w’tb  extremely 
complex  curvature. 

,  . .  —  —i  ■  ■■  — ■■  ■  o  2 

•The  equation  of  a  general  conic  section  is:  Ax“  *  Bxy+Cy  +Dx+Ey+F=0.  Thus  the 
e^uaiioa_of  a  general  qjjniC-Section  of  revolution  about  the  Z  axis  is  A  (x^+y^'HS 
z*/x^  -\v2  +  cz^  *  D,Jx*+  y^  +  ez+f=0. 
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The  Inner  aur^  .ce  of  the  bond- 
end  web  la  raproaantod  by  up 
to  five  general  oonto  seotiona 
of  revolution, 


The  case  head  la  an 
ellipse  of  revolution. 


The  aft  end  la  an  ellipse 
of  revolution  (concave  or 
convex),  or  a  cone. 


fhs  b*tv  *ght- through 
portion  of  dv  caw.  hi 
of  arbitrary  length  and 
radius,  thus  allow 'ng, 
inclusion  of  elliptical 
and  spherical  motors. 


Up  to  straight  lines,  circles, 
parabolas,  ell'pses,  hyperbolas 
or  any  co.-nbmation  thereof  describe 
tne  "straight-through*  portion  of  the 
inner  web  limiting  surface. 


FIGURE  lid.  SCHEMATIC  DIAOItAM  OF  A  ROCKET  MOTOR 
ILLUSTRATING  COMPUTER  r-'ROQRAM  CAPABILITIES 
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The  straight-through  portion  of  the  chamber  is  considered  to  have 
any  number  of  symmetrical  sections  (N  star  points),  As  shown  in  Figure  28, the  surface 
representing  one -half  of  one  symmetrical  section  (star  point)  may  consist  of  up  to  ten 
general  conic  sections,  This  is  sufficient  to  fit  virtually  stnv  singly  perforated  grain. 

The  second  major  step  in  applying  the  AGDA  technique  is  to 
introduce  a  pseudo-rectangular  grid  on  the  surface  of  the  chamber  described  above, 

Two  techniques  are  used:  one  for  the  head-end,  and  the  other  for  the  straight-through 
portion.  The  head-end  grid  is  formed  by  projecting  an  x-y  mesh  upward  from  the  x-y 
plane  (which  is  perpendicular  to  the  rocket  axis)  to  the  chamber  head  surface.  The 
grid  for  the  straight-through  portion  is  obtained  from  the  points  defining  the  periphery 
of  this  portion  together  with  points  equally  spaced  along  the  motor  axis.  Both  grids 
are  shown  in  Figure  29. 

Once  the  surface  is  defined  and  the  grid  is  superimposed,  the 
remaining  steps  of  Program  AGDA  imolve  simple  vector  integration,  differentiation, 
multiplication,  and  addition.  The  mathematics  involved  are  identical  to  those  described 
in  Section  B-2  of  this  report. 

b.  Input  Generator 

One  of  the  most  significant  innovations  of  the  program  is  the  input 
data  generator.  An  extremely  simple  method  of  description  of  the  geometrical 
peculiarities  of  an  arbitrary  motor  was  achieved.  All  parameters,  other  than  several 
control  and  identification  words,  are  directly  obtainable  from  a  standard  engineering 
drawing  or 


sketch 


For  example,  trc-><-  sections  of  several  classes  of  rccket  motors 


1  be  pseudo-rectangular  grid 
the  grain.  The  X  -  Y  mesh  i 


is  superimposed  on  a  half 
s  vis/ble  through  the  cutaw 


star  point  of  the  Inner  surface  of 
ay  portion  of  the  star  point. 


FIuURE  20,  EXAMPLE  OF  TYPICAL 


PSEUDO-RECTANGULAR  GRID 


readily  specified  to  AQDA  through  the  input  generator  are  shown  in  Figure  30, 

In  rare  instances,  there  may  be  need  for  the  evaluation  cf  a  motor 
configuration  which  is  not  amenable  to  specifications  using  the  input  generator.  To 
accommodate  this  possibility,  the  original  input  option  (specification  of  all  mathematical 
functions  in  generalized  elliptical  form)  may  be  utilized. 

The  frequent  engineering  requirement  for  a  grain  of  uniform  head 
web  thickness  has  called  for  individual  consideration.  Since  the  case  head  is  usually 
an  ellipse  of  revolution  and  ellipses  do  not  bum  into  ellipses,  a  special  iterative 
technique  was  devised  to  permit  rapid  calculations  of  the  core  head  coordinates  a 
constant  distance  from  the  case  head.  Amplification  of  the  theory  behind  this  technique 
is  given  in  Appendix  II. 

A  sample  input  data  sheet  describing  the  parameters  required  by 
the  input  data  generator  section  of  AGDA  is  shown  in  Figure  32.  Reference  is  made 
to  the  sketch  of  a  generalized  rocket  motor  (see  Figure  31)  and  Table  I. 
c.  System  of  Checks  and  Balances 

In  extremely  complex,  and  sometimes  continuing,  programs  such  as 
the  Advanced  Grain  Design  Analysis  Computer  Program,  the  need  for  extensive 
checking  under  a  wide  range  of  conditions  is  of  paramount  importance.  To  this  end, 
the  best  programming  practices  have  been  exercised  in  the  writing  of  the  program, 
and  to  a  certain  extent,  several  novel  programming  innovations  have  been  Included  to 
reduce  significantly  the  'lost"  computer  time  following  an  indication  of  system 
difficulty. 
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1.  Simple  Star 


FIGURE  30.  THE  AGDA  INPUT  GENERATOR  CON8IDER8 
ANY  OF  THESE  THREE  HALF  8TAR  POINT  CROS3  SECTIONS 


-70- 


FIGURE  31,  SCHEMATIC  REPRESENTATION  OF  CROSS  SECTIONS  OF  GENERALIZED  80LID 
PROPELLANT  ROCKET  AMENABLE  TO  SOLUTION  WITH  AG  DA  SHOWING 
REQUIRED  INPUT  PARAMETERS 
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TABLE  I 


INPUT  DATA  quantities  required  by  the  agpa 
INPUT  GENERATOR* 


RCH 

Ellipse  ratio  of  case  head  ellipse 

ZCH 

Z  translation  of  center  of  case  head  ellipse 

ACH 

Major  axis  ox  case  head  ellipse 

RCA 

Ellipse  ratio  of  case  aft  ellipse 

ZCA 

Z  translation  of  center  of  case  aft  ellipse 

ACA 

Major  axis  of  case  aft  ellipse 

PORT 

Diameter  of  igniter  port 

AFTCUT 

Distance  from  Z  a  0  to  aft  cutback 

Rl,  R2,  R3 

Fillet  radii  in  grain  cross  section 

T1 

Least  port  diameter 

T3 

Greatest  port  diameter 

SN 

Number  of  star  points 

HWEB 

Head-end  web  thickness  along  transverse  motor  axis 

ZNPCT 

Percent  web  between  printouts 

ZEND 

Type  of  run 

ZRND 

Control  parameter  for  core  rounds 

ZNCORE 

Number  of  core  head  equations  used 

*  8ee  Figure  32 


The  least  complex,  and  probably  the  most  frequently  used,  ts  the 


input  data  check  option.  In  this  type  of  run,  the  previously  described  input  is 
processed  by  the  input  generator  portion  of  the  program;  and  the  first  two  pages  of 
program  output  are  produced  in  entirety.  Contained  within  the  input  generator 
portion  of  the  program  is  a  series  of  data  compatibility  checks.  If  the  specified  data 
fell  to  meet  the  requirements  on  one  or  more  of  these  controls,  the  problem  will  be 
terminated  as  an  Input  data  check  run  regardless  of  the  run  type  requested  in  the  input 
data. 

Probably  the  most  crucial  check  is  that  which  is  made  on  the 
continuity  of  segment  equations  describing  the  central  perforation.  Since  a  primary 
control  variable  is  the  terminal  (maximum)  value  of  x  for  which  a  specific  segment 
equation  is  valid,  one  of  the  coordinates  assumed  to  be  common  to  consecutive  segments 
is  immediately  available.  The  y"  and  y+  values  are  then  obtained  from  the  left  and 
right  segment  equations  and  the  following  requirements  imposed: 

I — ZJLL  s  o.  0001 

y 

Before  this  check  was  inserted,  occasional  difficulty  was  experienced  in  the  transition 
from  a  given  segment  to  an  adjacent  arc  segment.  Essentially,  the  incremented  value 
x  was  not  valid  for  the  particular  equation,  since  it  was  external  to  theeurvty.  .  In 
general  engineering  work,  a  certain  tolerance  is  unconditionally  associated  with  each 
set  of  calculated  values;  but, unfortunately,  this  flexibility  is  not  available  in  the  exact 
analytical  representations  from  the  computer. 
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Several  of  the  numerous  additional  checks  accomplished  in  this 
section  of  the  program  are; 

1)  Case  equation  check  requiring  zero  coefficients 
for  the  radial  translation  terms  to  assure  a  sur-  ' 
face  of  revolution  about  the  transverse  motor 
axis 

2)  Comparison  of  successive  x  coordinate  cutoff 
values  for  the  segment  equations  to  assure  a  half 
star  point  specification  whioh  is  a  singly  defined 
function  in  x 

•4  *  .  •  <  ‘  ■ 

3)  Evaluation  of  the  need  for  the  head-end  web 
thickness,  HWEB,  for  use  in  the  inhibited  igniter 
port  calculation.  If  the  latter  is  requested, 

PORT^  0,the  value  of  KWEB  on  the  transverse 
motor  axis  will  be  calculated  if  the  input  value 
was  zero.  A  note  to  this  effect  along  with  the 
computer  value  are  included  as  output 

4)  For  tbe  constant  head-end  web  option,  installed 
for  use  with  elliptical  case  equations,  the  ellipse 
ratio  is  compared  with  unity.  If  an  identical 
comparison  is  found,  the  calculational  option  is 
altered  and  the  appropriate  spherical  head-end 
case  equation  generated. 

The  next  option  in  order  of  increasing  complexity  is  the  "Point 
Periphery  Tape"  option.  This  run  will  include  all  the  check  features  mentioned  above, 
and  in, addition,  a  great  number  of  checks  to  assure  proper  flow,  minimum  degree  of 
calculation  sufficient  for  the  particular  option  and  the  analysis,  and  the  credibility  of 
certain  key  parametric  evaluations. 

Another  critical  check  occurs  in  the  routine  which  calculates  the 
maximum  burning  distance  along  a  specific  normal.’  If  this  value  exceeds  a  pre¬ 
determined  large  value  (actually  out  of  conceivable  range),  pertinent  parameters  are 
noted  as  output  and  the  problem  is  terminated. 
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Of  the  numerous  additional  oheeks  made  in  this  part  of  the  program 
there  arfl  two  which  are  probably  most  significant:  the  variable  me?1'  check  and  the 
point  periphery  generation  oheck. 

The  function  of  the  variable  mesh  check  is  to  assure  the  fineness  of 
grid  structure  required  to  meet  the  input  tolerances  in  the  most  expedient  raanrer, 

The  point  periphery  generation  cheok  assures  that  a  minimal  number 
of  calculations  will  be  done.  The  number  of  x  grid  evaluations  is  somewhat  less  than  in 
the  complete  analysis,  effecting  a  relatively  small  saving.  On  the  other  hand,  the 
number  of  y  grid  evaluations  psr  x  grid  is  reduced  to  one,  the  required  periphery  point, 
as  opposed  to  the  five  to  twenty  that  might  be  required  for  the  full  analysis  of  the  head 
end.  Most  significant,  though,  is  the  fact  that  this  option  involves  no  surface  regression 
calculations. 

The  program  output  associated  with  the  Point  Periphery  Tape 
Option,  in  addition  to  that  mentioned  above,  includes  the  following: 

1)  The  coordinates  of  each  grid  point  on  the  half  star 
point  periphery,  at  the  intersection  of  the  core 
head  and  the  straight-through  section.  Values  are 
tabulated  in  cartesian  and  polar  coordinates. 

2)  A  point-wise  tabulation  of  the  point  periphery  vectors 
in  terms  of: 

a)  The  angle  between  the  straight-through 
normal  and  the  y  axis,  (S,  Y) 

b)  The  angle  between  the  head-end  normal 
and  the  Z  axis,  (H,  Z) 

c)  The  angle  between  the  straight-through 
and  head-end  normals,  (S,  H). 
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3)  Component  and  total  check  results  for  initial  sur- 
face  and  initial  chamber  volume,  The  only  param¬ 
eters,  other  than  input,  supplied  by  the  main  pro¬ 
gram  for  these  analyses  are  the  points  as  described 
in  1)  above. 

4)  Initial  port  area 

5)  Surface  area  depletion  for  core  rounds  (if  requested). 

Thus,  throughout  Program  AQDA  a  series  of  checks  and  balances 

has  been  established  which  will  ensure  that  the  program: 

1)  Will  not  proceed  any  farther  than  necessary  with 
faulty  input  data. 

2)  May  be  "saved"  in  the  event  of  some  electro¬ 
mechanical  failure  in  the  computer  tape  units. 

3)  Will  provide  a  completely  independent  set  of  check 
results  to  help  substantiate  the  actual  program 
output. 

With  respect  to  the  latter,  the  rocket  engineer  is  assured  that  if  the 
"check  results"  and  the  actual  compuf  tr  results  agree,  his  answers  are  correct, 
d.  Uniform  Temperature  Example 

For  a  complete  appreciation  of  the  versatility  and  ease  of  use  of 
Program  AGDA,  reference  is  made  to  the  following  problem.  The  particular  motor 
under  consideration  is  illustrated  in  Figure  33. 

For  ease  in  handling  of  the  data  and  to  ensure  the  reliability  of  the 
input,  the  data  are  used  immediately  in  the  input  generator  program.  Any  inconsistencies 
in  the  input  data  are  immediately  noted  in  the  output. 


FIGURE  33.  SCHEMATIC  OF  SOLID  PROPELLANT  ROCKET  MOTOR 
FOR  WHICH  FULL  COMPUTER  OUTPUT  IS  PRESENTED  !N  TABLES  II  -  XI 


The  program  output  for  tho  complete  rocket  motor  evaluation  under 
consideration  is  shown  in  Tables  II  to  XI. 

Table  II  -  the  first  page  of  computer  output  simply  lists  the  input 
data  "exactly  as  It  appears  on  the  input  data  cards, "  but  with  variable  names  printed 
above  each  word. 

Tables  III,  IV,  and  V  are  the  "Point  Periphery  Tape"  output.  This 
tape  contains  the  coordinates  and  direction  cosines  which  describe  the  head-end  normals 
emanating  from  the  point  of  intersection  of  the  core  surface  of  revolution  and  the 
straight-through  star  design.  For  output  purposes,  the  direction  cosines  are 
expressed  as  angles  relative  to  the  coordinate  axes.  This  is  shown  pictorially  in 
Figures  34  and  35. 

At  the  conclusion  of  the  Point  Periphery  Tape  output,  the  supplementary 
check  results  are  printed  (see  Table  VI).  These  results  are  obtained  from  Point 
Periphery  Tape  data  only,  and  no  use  is  made  of  the  bulk  of  the  Advanced  Grain  Design 
Analysis  program.  Initial  surface  area  and  initial  chamber  volume  are  calculated 

The  next  three  pages  of  output  (see  Ta»  !<■>  v'i J.  VIII,  and  IX)  are  :u 
the  original  output  format.  Table  VI  shows  burning  surface,  chamber  volume,  and  the 
first  moment  of  the  propellant  along  the  Z  axis.  Table  VIII  shows  the  moments  of 
inertia,  while  Table  IX  shows  the  products  of  inertia. 

The  final  two  pages  oi  computer  output  (Tables  X  and  XI)  represent 
the  new  output  data  format.  These  pages  are  identical  except  that  Table  X  is 
independent  of  propellant  density,  while  Table  XI  utilises  this  quantity  for  determination 
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G6NERAIUS0  GRAIN  DESIGN  PROGRAM 


GRAIN  DESIGN  INPUT  GENERATOR 

COMPLETE  SOL  10  PROPELLANT  ROCKET  MOTOR  EVALUATION  TEST  CAib 


INPUT 


CONTROL  PARAMETERS 

*  ZNPCT 

>•0 

•  ZENO 

s  ZftNO 

•  ZNCORE 

e 

CASE  PARAMETERS 

AC 

•  AC 

*  zc 

•  l)I  AM 

•  ACA 

•  RCA 

• 

Z2.025 

1*4 

15.7 

44.03 

22.025 

1.4 

CROSS  SECTION 

PARAMETERS 

*  Rt 

• 

70 

N 

»  R3 

*  ri 

*  T3 

•  SN 

• 

•  A3 

0*1 

0.75 

13.A0C5 

30.85 

9.0 

ntf.na 

*  THETIM 

•  THET20 

»  THfcTZM 

»  THET3D 

•  ‘  THEM* 

• 

■  0*0 

iS.O 

ll.O 

37.0 

2*0 

2N.3 

^  CASE# 
/(•5*?6 

»  P&5ECG 
*•15*56 

*  PROPD 
0*0632 

•  PCRT  •  APTCUT  •  HWt*  • 
10. 0  9.735  6.6 

«  DEL  ft 

o  5 

8  DfcXY 

0*3 

•  DELZ 
0.5 

*  UELTHO  »  OELTHK  • 

5.C 

TABLE  II.  REPRODUCTION  OF  INPUT  DATA  TO  PROGRAM  AGDA 
FOR  ROCKET  SHOWN  IN  FIGURE  33 
(Blank  spaces  indicate  zero) 
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GENERAL IZSD  GRAIN  DESIGN  PROGRAM 
COMPLETE  SOLID  PROPELLANT  ROCKET  MOTOR  EVALUATION  TEST  CASE 
CHECK  OUTPUT  AND  CORE  ROUND  CORRECTION 


X 


0.000000 
0. 100000 
0. 600000 
0.724264 
0.900000 
1.200000 
1.443528 
t. 500000 
1.800000 
2.100000 
2.172792 

2.400000 

2.700000 
2.897056 
3.000000 
1.300000 
3.600000 
3.621320 
3.900000 
4.200000 
4.345585 
4.500000 
V8000CO 
5.069849 
5.100000 

5.400000 

9.700000 
8.794113 
6.000000 
6.190224 

6.200973 

6.213159 

6.226756 

6.241735 

6.258063 

6.275703 

6.294610 

6.306948 

6.313560 

6.333531 

6.354496 


V 

l 

R 

THt  t  A 

(StY) 

(H,n 

<  S»H) 

0.000000 

24.832141 

O.OOCCUO 

22*5000 

90.0000 

0.0000 

90.0'-'0n 

0.124264 

24.829965 

0.324718 

22.5000 

90.0000 

0, 7676 

90.0000 

0.248328 

24.821436 

C. 649435 

22.3000 

90.0009 

1.5363 

90*'W.0 

0.300000 

24.819454 

0.783938 

22.5000 

90.0000 

1.8548 

90.0000 

0.372792 

24.812543 

0.974153 

22.5000 

90.0000 

2.3059 

90.000U 

0.497056 

24.797279 

1.298371 

22.5000 

90.0000 

3.0776 

90.0000 

0.600000 

24.781306 

1.567876 

22.5000 

9C.OOCO 

3.7189 

90.0000 

0.621320 

24.777620 

1.623388 

22.5000 

9C.0000 

.>*8519 

90.0000 

0*745584 

24.753543 

1.948106 

22.5000 

90.0000 

4.6294 

90.0000 

0.869848 

24.725020 

2*273024 

22.5000 

90.0000 

5.4111) 

90.0000 

0.900000 

24.717425 

2*351813 

22.5000 

90.0000 

3.6011 

90. COCO 

0.994113 

24.692014 

2.597741 

22.5000 

90.0000 

6.1972 

90.U000 

1.118377 

24.654485 

2.922459 

22.5000 

90.0000 

6.9889 

90.0000 

1.200000 

24.627350 

1.135751 

22.5000 

90.0000 

7.5122 

90.0000 

1.242641 

24.612383 

3. 2471 76 

22.5000 

90.0000 

7.7867 

90.0000 

1.366905 

24.565657 

3.571H94 

22.5000 

90.00'i0 

8.3916 

90*0000 

1*491169 

24.514243 

3.896612 

22*5000 

90.0000 

9.404? 

90.0000 

1.500000 

24.310409 

3.919639 

22.5000 

90.0009 

9.4623 

90.0000 

1*615433 

24.458071 

4.221330 

22.5000 

90.0000 

1U.225& 

90.0000 

1.739697 

24.397063 

4.546047 

22.5000 

oo.opro 

11.0566 

90.0000 

1.800000 

24.165692 

4.703627 

22.5000 

99.0C00 

11.4637 

90.0000 

1.863961 

24.331139 

4.870765 

22.5000 

90.0000 

P.8932 

90.0000 

1.988225 

24.260193 

5.195482 

22.5000 

90.0000 

12.7516 

90.0000 

2.100000 

24.192010 

5.487565 

22o  5000 

90.0000 

11.5301 

90.000C 

2.112489 

24.184130 

5.520200 

22.5000 

90.0000 

11.6179 

90*0000 

2.236753 

24.102825 

5.544910 

22.5000 

90.0000 

14.4981 

90.0000 

2.361017 

24.016152 

6.169636 

22.5000 

90.0000 

15.39)7 

90. coo 0 

2.400000 

23.987831 

6.271502 

22.5000 

90.0000 

15.6781 

90.0000 

2.485281 

23.923967 

6.494353 

22.5000 

90.0000 

16.3061 

90.0000 

2.564C77 

23.862589 

6.70C251 

22.5000 

67.5003 

16*8941 

73.1059 

2.540124 

23.362342 

6.701067 

22.2756 

64.157C 

16.8965 

73.1373 

2.516738 

23.861594 

6.703530 

22.0512 

60.7989 

16.90)5 

73.2318 

2.494008 

23.860341 

6.707649 

21.8276 

57.4261 

16.9154 

‘  71.3899 

2.472022 

23.858581 

6.713430 

21.6059 

54*0383 

16.932C 

73.6118 

2.450867 

23.85631 2 

6.720871 

21.3870 

50.6354 

16.9535 

73.8974 

2.430629 

23.853541 

6.729963 

21.1718 

47.2166 

16.979o 

74.246) 

2.411390 

23.850259 

6.740690 

20.9612 

43.781b 

1 7 . 0 105 

74.6)80 

2.400000 

23.847974 

6.748154 

20.8335 

41.6435 

17.0)21 

74.9444 

2.394236 

23.846708 

6. 752288 

20.7678 

40.5266 

17.0440 

71.103  0 

2.373125 

23.842719 

6. 765286 

20.5802 

37.2500 

17.0315 

75*601 7 

2.363106 

23.838291 

6.779667 

20.3991 

33.9730 

17.1231 

76*1536 

TABLE  1!I.  REPRODUCTION  OF  AG  DA  "POINT  PERIPHERY 
TAPE" OUTPUT 
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GENERALISED  GRAIN  DESIGN  PROGRAM 
COMPLETE  SOLID  PROPELLANT  ROCKET  MOTOR  EVALUATION  TEST  CASE 


8 

V 

2 

R 

THETA 

I  S»Y) 

(Ht2) 

1S»H) 

6.3*6430 

2.149223 

23.833432 

6.795420 

20.2250 

30.6698 

17.1686 

76.7583 

6.399276 

§•336561 

23.623146 

6.812507 

20*0386 

27.3432 

17.2180 

77.4144 

6.422972 

2.325173 

23»  822445 

6.830885  ' 

19,9007 

23.9954 

17*2712 

78.1206 

6.447448 

2.313128 

23.816332 

6.850504 

19.7520 

20.6262 

17.3281 

78.8749 

6.472633 

2.306491 

23.809823 

6.871310 

19.6133 

17.3353 

17.3886 

79.6754 

6.498435 

2.299316 

23.802949 

6.893241 

19.4851 

13.8230 

17.4523 

80.5205 

6.324826 

2.293660 

23.795703 

6.916230 

19.3680 

10.3884 

17.5194 

81.4076 

6.331664 

2.289572 

23.788119 

6.940205 

19.2627 

6.9313 

17.5893 

82.3344 

6.570879 

2.287099 

23.780214 

6.965089 

19.1696 

3.4515 

17.6620 

83.2994 

6.606373 

2.286284 

23.772009 

6.990800 

19.0893 

-0.0512 

17.7373 

84.2968 

6.623244 

2.286614 

23.766868, 

7.006851 

19,0468 

-2.1998 

17.7844 

64.9216 

6.923244 

2.298138 

23.672172 

’  7.294707 

18.3633 

-2.1998 

18.6373 

84.5961 

7.223244 

2.309663 

23.572312 

7.583522 

17.7319 

-2.1998 

19.5112 

84.8689 

7.523244 

2.321188 

23.467124 

7.873190 

17.1469 

-2.1998 

20.4075 

84.8401 

7.823244 

2.332713 

23.356410 

8.163620 

16.6034 

-2.1998 

21.3278 

84.8093 

8.123244 

2.344238 

23.239977 

8.434735 

16.0973 

-2.1998 

22.2739 

84.7787 

8.423244 

2.353763 

23.117592 

8.746465 

15.6249 

-2.1998 

23.2479 

04.7420 

8.723244 

2.367287 

22.988992 

9.038752 

15.1831 

-2.1998 

24.2522 

84.7051 

8.790033 

2.369846 

22.959464 

9.103912 

15.0885 

-2.2002 

24.4802 

84.6968 

8.796120 

2.369896 

22.956009 

9.109782 

15.0789 

1.2763 

24.5008 

83.2938 

8.802218 

2.369575 

22.954161 

9.115587 

15.0670 

4.7762 

24.5212 

81.9012 

8.810253 

2.368580 

22.950736 

9.123088 

15,0478 

9.4166 

24.5476 

83.0930 

9.110253 

2.318826 

22.821029 

9.400727 

14.2802 

9.4166 

25.5395 

80.0219 

9.410253 

2.269072 

22.684520 

9.679956 

13.5568 

9.4166 

26.5690 

79.9461 

9.710253 

2.219317 

22.540878 

9.960642 

12.8741 

9.4166 

27.6394 

79*8653 

10.010253 

2.169563 

22.389695 

10.242663 

12.2288 

9.4166 

28.7544 

79.7791 

10.310253 

2.119809 

22.230522 

10.525915 

11.6182 

9.4166 

29.9135 

79,6870 

10.429692 

2.100000 

22.164814 

10.639007 

11.3842 

9.4166 

30.3967 

79.6485 

10.610253 

2.070054 

22.06283,9 

10.810299 

11.0397 

9.4166 

31.1369 

.79.5804 

10.910233 

2.020300 

21.886056 

11.095730 

10.4909 

9.4166 

32.4158 

79.4*23 

11.210253 

1.970546 

21.699459 

11.382127 

9.9696 

9.4166 

33.7628 

79.3695 

11.510253 

1.920792 

21.502233 

U.  669419 

9.4740 

9.4166 

35.  1866' 

79.2476 

11.810253 

1.871037 

21.293379 

11.957543 

9.0023 

9*4166 

36.6984 

79.1161 

12.110253 

1.321283 

21.071704 

12.246439 

8.5527 

9.4166 

38.3119 

78.9718 

12.238580 

1.800000 

20.972610 

12.370240 

8.3668 

9.4166 

39.0372 

78.9093 

12.410253 

1.771529 

20.835731 

12.536055 

8.1239 

9.4166 

40.0441 

78.8M4 

12.710253 

1.721774 

20.583648 

12.826341 

7.7145 

9.4166 

41.9170 

78.66 1 1 

13.010253 

1.672020 

20.313044 

13.117253 

7.3233 

9.4166 

43.9594 

73.4668 

13.310233 

1.622266 

20.020780 

13.408749 

6.9490 

9.4166 

46*2119 

78.2637 

13.610253 

1.572511 

19.702569 

13.700794 

6.5907 

9.4166 

4P.7283 

78.0332 

13.910253 • 

1.522757 

19.352159 

13.993352 

6.2473 

9.4166 

31.5909 

77.7830 

TABLE  IV,  HE  PRODUCTION  OF  PAGE  2  OF  AGDA  "POINT  PERIPHERY 

TAPE”  OUTPUT 
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GENERALIZED  GRAIN  DESIGN  PROGRAM 


COMPLETE  SOL  10  PROPELLANT  ROCKET  MOTOR  EVALUATION  mt  CASE 


X 

V 

z 

R 

THETA 

(S,Y) 

IH.23 

tS.M 

14.047463 

1.500000 

19.178724 

14.127326 

6*0930 

9.4166 

>3.0480 

77.6*97 

14.210233 

1.473003 

16.959932 

14.286392 

3.9180 

9.4166 

54.9280 

7  7.  6'JOO 

14.510233 

1.423248 

18.309560 

14.379883 

3*6020 

9.4166 

58.9621 

77. ITU 

14.763834 

1.377873 

18.021092 

14.347904 

3.  324? 

9.4166 

63.6179 

?*»*  ;234 

14.828248 

1.369126 

17.931823 

14.891320 

3.2793 

12.8753 

64.5032 

73. 6635 

14.872167 

1.337677 

17.840492 

14.934009 

3.2161 

16.3431 

63.420  7 

70.4782 

14.915423 

1.343354 

17.747043 

14.975814 

5.1472 

19.8200 

66.3712 

67.2500 

14.937839 

1.326786 

17.651610 

15.016387 

3.0690 

23.3062 

67.3550 

6  3.9)51 

14.999306 

1.307408 

17.334289 

15.036177 

4.9816 

26.8020 

68.3718 

63.6A38 

15.039606 

1.263478 

17.453133 

13.094442 

4.8834 

30*3076 

69.4218 

57.3465 

13.078603 

1.261048 

17.354268 

19.131242 

4*7806 

33.8236 

TO. 5046 

53.97.30 

19.116142 

1.234188 

17.231702 

13.166441 

4*6677 

37.3304 

71.620  7 

39. > jJ6 

15.152073 

1.204966 

17.147496 

13.199909 

4*3469 

40.3891 

72.7701 

47.1178 

13.137733 

1.200000 

17.130406 

19.203179 

4.3263 

41.4693 

72.9600 

46.5542 

15.166230 

1.173473 

17.041721 

1 5.231320 

4,4186 

44*4400 

73.9524 

43.6167 

15.218533 

1.139760 

16.934473 

15.261154 

4.2831 

48.0034 

75.1668 

40.1153 

15.247043 

1.106174 

16.832222 

13.287116 

4.1493 

31.3725 

76.339') 

36.7720 

15.273603 

1.070834 

16.726827 

15.311093 

4*0105 

54.7490 

77.53’9 

33.3734 

15.298086 

1.033945 

16.624339 

15.332986 

3.8665 

58.1331 

78.7600 

lO.O^U 

15.320389 

0.993584 

16.319368 

15.332703 

3.7181 

61.3234 

80.0054 

26.5813 

13.340196 

0.455933 

16.413635 

13.370130 

3.9658 

64.9168 

81.2688 

23.1432 

15.337996 

0.913242 

16.307463  15.383242 

3.4104 

68.3041 

82.5481 

17.6974 

15.363869 

0.900000 

16.268132 

13.390206 

3.3525 

67.9520 

83.024? 

13.4245 

15.373234 

0.873366 

16.200389 

15.398021 

3.2313 

71.7089 

83.8476 

16.2224 

15.366031 

0.830378 

16.09227 2 

13*408421 

7.0892 

73.1359 

33.1673 

.12.7148 

15.396291 

0.786414 

13.982936 

13.416361 

2.9240 

78*3853 

86.5094 

.  9.1753 

15.403937 

0.741617 

15.872899 

13*421778 

2.7564 

82.0373 

87.8650 

5.6072 

15.408868 

0.696139 

13.762443 

15.424584 

2.586? 

85*3327 

89.2295 

2.0140 

15.410416 

0.670597 

13.707111 

15.424999 

2.4917 

87*5083 

39.9121 

0.0869 

15.413326 

0.600000 

15.707111 

15.424999 

2.2292 

87*7707 

89.9121 

O'.  08u9 

15.422082 

0.300000 

15.707111 

13.424999 

1.1144 

18*8836 

89.9121 

0.0849 

15.423409 

0.221540 

15.707111 

15.424999 

0.8229 

69.1771 

89.9121 

0.0840 

15.425000 

0.000000 

13.707111 

15.424999 

0.00(30 

90* 0003 

89.9121 

’  C.C ^88 

TABLE  V.  REPRODUCTION  OP  PAGE  3  OP  AGDA  "POINT 
PERIPHERY  TAPE"  OUTPUT 


SURFACE 

VOLUME 


GENERALIZED  GRAIN  DESIGN  PROGRAM 
COMPLETE  SOLID  PROPELLANT  ROCKET  MOTOR  EVALUATION  TEST  CASE 

I 

SUPPLEMENTARY  CHECK  RESULTS 


TOTAL 

HEAD  END 

TOP  SECT* 
ST*  THRU 

CENTER  SECT* 

IT*  THRU 

APT  END 

4022.39 

509.36 

878*03 

26)4*79 

0*00 

9434.20 

2899*71 

63)4.57 

0*00 

PORT  AREA 

a  416*0262 

TABLE  VI.  CHECK  RESULTS 
(These  results  are  based  only  upon  "Point  Periphery  TapeH 
data  and  therefore  serve  to  validate  the  majority  of  AGDA. ) 
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CCMPUTEO  TRACS  OP  INTERNAL  GEOMETRY 


COMPLETE  SOUO  PROPELLANT  ROCKET  MOTOR  EVALUATION  TEST  CASE 


WEB 

Chamber 

CHAMBER 

FIRST 

l 

BURNOUT 

SURFACE 

VOLUME 

MOMENT 

(PERCENT) 

(SQ-IN) 

(CU-lN) 

(!N«*4) 

0 

5.57476 

03 

1.71516 

04 

3.8688b 

05 

S 

S.7019E 

03 

1.9012E 

04 

3.73026 

03 

io 

5.6624E 

03 

2.09016 

04 

3.58626 

05 

ts 

5.52676 

03 

2*27476 

04 

3.44206 

05 

20 

5.40466 

03 

2.45466 

04 

3.29766 

05 

25 

5.42026 

03 

2.63316 

04 

3.15096 

05 

30 

5.44076 

03 

2.8123E 

04 

J.OOOlr 

05 

35 

5.45336 

03 

2.99206 

04 

2.8449  •: 

05 

AO 

5.46206 

03 

3.17206 

04 

2.6350c 

05 

45 

5.46666 

03 

3.35226 

04 

2.52016 

05 

50 

5.46636 

03 

3.53256 

04 

2.34996 

05 

55 

5.4621E 

03 

3.71276 

04 

2.17416 

OS 

60 

5.4524E 

03 

3.89266 

04 

1.9921E 

05 

65 

5.43736 

03 

4.07226 

04 

1.80396 

03 

70 

5.41506 

03 

4.25116 

04 

1.60896 

05 

75 

5.38586 

03 

4.4292E 

04 

1.40666 

05 

80 

5.35186 

03 

4.60606 

04 

1.19716 

03 

85 

5.3232E 

03 

4.78186 

04 

9.79326 

04 

90 

3.30156 

03 

4.95686 

04 

7.53056 

04 

95 

5.26956 

03 

5.13096 

04 

5.17996 

04 

100- 

5.22736 

03 

1004 

4.44556 

03 

5.3027E 

04 

2. 75636 

04 

105 

2.45106 

03 

5.4002J 

04 

1.73046 

04 

no 

1.80146 

03 

5.4692E 

04 

1.05956 

04 

115 

1.29436 

03 

5.5193E 

04 

5.95466 

03 

120 

8.7868E 

02 

5.55466 

04 

2.83B&E 

03 

125 

5. 096.36 

G2 

5.57716 

04 

9.69546 

02 

130 

1.67186 

02 

5.58786 

04 

1.22556 

0? 

WEB 

THICKNfcSS  ■ 

6.600 

HEIGHT  « 

15.700 

DIAMETER  « 

44.050 

CASE 

i  VOLUME 

ta 

35894.089 

INHIBITED  PORT  DIAMETER  « 

10.000 

TABLE  VII.  PAGE  1  OF  SUMMARIZATION  OF  PROGRAM  AGDA  OUTPUT 
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C08PUT60  TRACE  OP  MOMENTS  OP  INERTIA 
COMPUTE  SOUQ  PROPELLANT  ROCKET  MOTOR  EVALUATION  TEST  CASE 


WEB 

BURNOUT  I-XX  I-YV  1-22 


(PERCENT) 

( IN«*5) 

t!N»»5) 

UN«*5) 

0 

1*34646 

07 

1.34646 

07 

i*0534b 

07 

5 

1*30546 

07 

1.30346 

07 

1*02806 

07 

10 

1*26196 

07 

1*26196 

07 

1*00006 

07 

19 

1*21716 

07 

1.21716 

07 

9*  70Q?h 

Ob 

20 

1*17106 

07 

1.17106 

07 

9*  3811h 

06 

25 

1*12306 

07 

1*12306 

07 

9*03946 

06 

30 

1*07266 

07 

1.07266 

07 

8.67346 

06 

39 

1*01996 

07 

1.01996 

07 

8.28336 

U6 

40 

9*64646 

06 

9.64846 

06 

7.869i> 

06 

43 

9*07326 

06 

9.07326 

06 

7.43136 

06 

30 

8.4736E 

06 

8.47366 

06 

6.97CU 

06 

35 

7*84916 

C6 

7.8491S 

06 

6.48596 

06 

bv 

7*19946 

06 

7*19946 

06 

5*9794*: 

06 

63 

6*52456 

06 

6.52436 

06 

5.4513b 

06 

70 

3*82406 

06 

5.82406 

06 

4*90266 

06 

73 

5.09756 

06 

5.09756 

06 

4*33456 

06 

80 

4*34656 

06 

4*3465E 

06 

3.74916 

06 

85 

3*56906 

06 

3.56906 

06 

3.14646 

06 

90 

2.7642E 

06 

2.76426 

06 

2.32666 

06 

93 

1*93186 

06 

1*93186 

06 

1.89056 

06 

too 

1.07896 

06 

1.07896 

06 

1.2447* 

06 

105 

6.8974E 

05> 

6.89746 

03 

8*38116 

05 

UO 

4.28276 

05 

4.28276 

09 

5.41976 

05 

115 

2.44066 

05 

2.44066 

05 

3.21246 

03 

120 

1.18196 

05 

1.18196 

05 

1.62096 

05 

125 

4.07066 

04 

4.0706E 

04 

5. 8056tr 

04 

130 

4.04896 

02 

6.30556 

01 

8.59046 

01 

TAfeLE  Vin.  PAGE  2  OF  SUMMARIZATION  OF  PROGRAM  AGDA 

OUTPUT 
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COMPUTED  TRACE  fiF  PRODUCTS  OF  INERTIA 


COMPLETE 

SOLID  PROPELLANT 

ROCKET  MOTOR 

evaluation  TEST 

WEB 

BURNOUT 

I-XY 

I-XZ 

I-Y2 

IPERCENTJ 

I IN*«5» 

UN«*5> 

1 1 N  *  •  5  ) 

0 

8.24416-02 

3.4223E-02 

1.23206-01 

$ 

8.03416-0? 

3.33526-02 

1.20076-01 

10 

7.80646-02 

3,24146-02 

1. 16696-01 

15 

7,56336-02 

3.14326-02 

1.13156-01 

20 

7.30516-02 

3.04086-02 

1.0947E-01 

25 

7.03016-02 

2.93306-02 

1.05596-01 

30 

6.73706-02 

2.81916-02 

1.01496-01 

35 

6.42596-02 

2. 69866-02 

9.7150t-02 

♦0 

6.0969E-02 

2.57116-02 

9.25596-02 

45 

5.75026-02 

2.43626-02 

8.77036-02 

50 

5.3861E-02 

2.29356-02 

8.2567c-Q2 

55 

5.00466-02 

2.14276-02 

7.71376-02 

60 

4.60626-02 

1.98326-02 

7.13966-02 

65 

4. 1912E-02 

1.81476-02 

6.5329u-02 

70 

3.76036-02 

1.63666-02 

5.8917c-02 

75 

3.31406-02 

1.44836-02 

4.21386-02 

80 

2.85386-02 

1.24946-02 

4.49796-02 

65 

2.3797E-02 

1.03906-02 

3.7403c-02 

90 

1.89186-02 

8. 16366-03 

2. 93896-02 

95 

1.3912E-C2 

5.80956-03 

2.C9146-0? 

100 

8.8328E-C3 

3.33696-03 

1.20136-02 

105 

5.8069E-03 

2.12416-03 

7.6468r- 03 

110 

3.6828E-C3 

1.31356-03 

4.7287C-03 

115 

2.1436E-03 

7.4426E-U4 

2.  .37936-03 

120 

1.0637E-C3 

3. 5736E-04 

1.28656-03 

125 

3.75226-04 

1.2273E-04 

4.4201.— 04 

130 

0. 00006  or. 

1.5602E  02 

O.OOOLc  00 

TABLE  IX.  PAGE  3  OF  SUMMARIZATION  OF  PROGRAM  AGDA 

OUTPUT 
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THlOKOL-ELRTnN 

AOVANCEC  CRAIN  OfiS IGN  COMPUTER  PROCRAM 
COMPUTE  SOLID  PROPELLANT  ROCKET  MOtOR  EVALUATION  TEST  CASE 


THEORETICAL  ANALYSIS  BY  F.E. MOORE 
PROGRAMMING  ANALYSIS  RY  O.H. FREDERICK 


BURNING 

CHAMBER 

PROP. 

PROP*  ■ 

PROPELLANT 

INERTIA 

.  URN 

SURFACE 

VOLUME 

VOLUME 

CO 

I-XX 

I-YY 

■'(.t 

#  t  ri*»f  u  \ 

i  nl  tu( 

1#*  1 1—  t  Aj  1 

Mu  1  * 

/ 

1 

ZJ 

o 

• 

<L 

I CU— IN  J 

l tU“  IN l 

\  INJ 

0 

5574.66 

17151.03 

38743.06 

9.986 

93500.6 

93500.6 

3 

5701.59 

19011.64 

36882.45 

10.114 

90655*6 

90655.6 

10 

5662.43 

20901.43 

34992c  66 

10*248 

87634*9 

87634. V 

11 

9526.67 

22747.39 

33146.70 

10*384 

64521.2 

64521.2 

20 

5404.55 

24546.06 

31348.03 

10*319 

61318.8 

81318.8 

29 

5420.16 

26330.97 

29563.12 

10*658 

77984.9 

77984.9 

10 

9440.68 

28122.73 

27771.36 

10.603 

74488*9 

74488.9 

19 

9453.31 

29919.56 

25974.33 

10.953 

70828.9 

70828.9 

40 

9461.95 

31719.94 

24174.13 

11.107 

67002.6 

67002.6 

49 

5466.61 

33322.22 

22371.87 

11.263 

63006.6 

63008.6 

70 

5466.34 

35324.93 

20369.16 

11.425 

58644.7. 

58844.7 

13 

5462.09 

37126.86 

18767.23 

11.584 

54308.0 

54306.0 

oO 

5452.42 

38926.37 

16967.72 

11.741 

49996.1 

49996.1 

63 

5437.26 

40721.74 

13172.35 

11.889 

43308.7 

45308.7 

10 

5415.02 

42510.87 

13383.22 

12.021 

40444.2 

40444.2 

/  9 

5353.84 

44291.97 

11602.12 

12.124 

35399.6 

35399.6 

:>0 

9351.81 

46039.71 

9834.38 

12.173 

30184.3 

30134.3 

.5 

5323.25 

47817.55 

8076.34 

12.126 

24784.6 

24784.6 

lO 

5301.46 

49567.61 

6326.48 

li.°03 

19195.9 

19195.9 

■i  5 

5269.46 

51309.21 

4384.88 

11.298 

13415.4 

13415.4 

1  JO 

9227.27 

53026.60 

2867.49 

9*612 

7492.2 

7492.2 

.  *X 

5701.89 

53026.60 

38743.06 

12*173 

93500.6 

93500.6 

MN 

5227.27 

17131.03 

2867.49 

9*612 

7492.2 

7492.2 

A  VO 

5438.34 

33269.35 

20624.73 

li*140 

55739.0 

55739.0 

1 J04 

4445.55 

1'5 

2450.99 

54001.86 

1892.23 

4789.8 

4789.8 

i  10 

1801.36 

54691.78 

1202.31 

2974.1 

2974.1 

1  18 

1294.28 

55192.89 

701.20 

1694.9 

1694.9 

120 

876.65 

55545.79 

348.30 

820.8 

820.8 

123 

909.63 

55771.24 

122.85 

282.7 

282.7 

« Id 

167.18 

55878.48 

15.61 

2.8 

0.4 

PER 

CENT  PROPELLANT  (  0- 

100  PCT  NEB ) 

BY  VOLUME 

*  64.18 

PER 

CENT  PROPELLANT  < 100- 

130  PCT  WEB) 

BY  VOLUME 

-  5. 13 

TABLE  X»  PAGE  1  OF  NEW  PHOGRAM  AGDA  OUJPUT 

(The  max, ,  min. ,  and  average  values  refer  only  to  O-AlOO  pet  web. 
The  100  +  surface  area  ts  the  surface  just  after  web  burnout.) 


TENSOR 

t-27. 


73184. 8 
7 1 3  6  fi  •  5 

69447.8 
67368*9 

65146.9 

62773.9 

60232.1 

57522.7 
SA&A6.6 
51606.0 
4J40J.? 
A50A9.M 
4152J.6 

37536.2 
34046.0 

30100.9 

26039.6 

21890.3 

17549.6 

13128.6 

8644.1 

73154.8 

8644.1 

45593.2 


5820.2 
376  3.7 

2230.9 

1125.6 

403.2 

0  #  fi 
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TH1QK0L-ELRT0N 

ADVANCED  GRAIN  DESIGN  COMPUTER  PROGRAM 

COMPLETE  SOLID  PROPELLANT  ROCKET  MOTOR  EVALUATION  TEST  CASE 

THEORETICAL  ANALYSIS  BY  P.E.HOQRE 
PROGRAMMING  ANALYSIS  BY  O.H. FREDERICK 


WtB 

BURNING 

PROP. 

SYSTEM 

SYSTEM 

PROPELLANT 

INERU4 

TENSOR 

i  URN 

SURFACE 

HEIGHT 

WEIGHT 

CG 

I-XX 

I-YY 

1-2/ 

POT 

1 SQ-IN) 

(LBS) 

(LBS) 

UN) 

- (SO- 

“P  i  »L8S  )* 

C 

5574.66 

2448.56 

3154.32 

7.6C4 

5909.2 

5909.2 

4623.4 

s 

5701.80 

2330.97 

3036.73 

7.610 

5729.4 

5729*4 

4511*6 

10 

5662.43 

2211.54 

2917.30 

7.609 

5538.5 

i  5538.5 

4389.1 

15 

5526.67 

2094.87 

2600.63 

7.601 

5341.7 

5341.7 

4257.5 

20 

5404*55 

1981.20 

2686.96 

7.583 

5139.3 

5139.3 

4117.3 

25 

5420.16 

1868.39 

2574.15 

7.555 

4928.6 

4928.6 

.  3967.3 

30 

5440.68 

1755.15 

2460.91 

7.515 

4707.7 

4707.7 

3806.7 

35 

5453.31 

1641.59 

2347.35 

7.461 

4476.4 

4476*4 

3635.4 

40 

5461.95 

1527.61 

2233.57 

7.389 

4234.6 

4234.6 

3453.7 

45 

5466.61 

1413.90 

2119.66 

7.294 

3982.1 

3982.1 

3261.5 

30 

5466.34 

1299.97 

2005.73 

7.172 

3719.0 

3719.0 

3099.1 

■  53 

5462.09 

1186.09 

1891.85 

7.017 

3444.9 

3444*9 

..  2346*6 

*0 

5452.42 

1072.36 

1778.12 

6.819 

3159.8 

3159.8 

2624.3 

65 

5437.26 

958.89 

1664,65 

6.569 

2863.5 

2863.5. 

2392.5 

70 

5413.02 

645.82 

1551.58' 

6.253 

2556.1 

2556.1 

2151,7 

75 

5365.84 

733.25 

1439.01 

5.854 

2237.3 

2237.3 

1902.4 

riO 

3351*61 

621.53 

1327.29 

5.349 

1907.6 

1907.6 

1645.4 

.35 

5323.25 

510.44 

1216.20 

4.706 

1566.4 

1566*4 

1380.9 

10 

5301.46 

399.83 

1105.59 

3.883 

1213.2 

1213.2 

1108.9 

)3 

5269.46 

289.76 

995.52 

2.820 

847.9 

847*9 

879.7 

too 

5227.27 

181.23 

886.99 

1.439 

473.5 

473.5 

546.3 

MX 

3701.89 

2448.56 

3154.32 

7.610 

5909.2 

5909.2 

4623.4 

MIN 

5227.27 

181.23 

886.99 

1.439 

473.5 

473*5 

546.3 

■WO 

5436.34 

1303.48 

2009.24 

6.338 

3522.7 

3522*7 

2881.5 

1004 

105 

4445.55 

2450.99 

119.59 

825.35 

0.761 

302.7 

302*7 

367.8 

110 

1801.36 

75.99 

781.75 

0.261 

183.0 

188.0 

237.9 

115 

1294.26 

44.32 

750.08 

-0.119 

107.1 

107.1 

141.0 

120 

678.68 

22.01 

727.77 

-0.394 

31.9 

51.  ) 

Tie  1 

125 

509.63 

7.76 

713.52 

-0.567 

17.9 

17.9 

25.5 

no 

167.18 

0.99 

706.75 

-0.648 

0.2 

0.0 

0.0 

TABLE  XI.  PAGE  2  OF  NEW  PROGRAM  AG  DA  FINAL  OUTPUT 

(This  page  includes  the  weight  of  Inert  parts  as  well  as  propellant 
density  to  obtain  actual  weight,  CG  travel,  etc;) 
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Port  Head 


Straight  Through  Section 


The  core  or 
central  perfora¬ 
tion  of  a  solid 
propellant  rocket 
motor. 


Star  Points 


[Angle 

fckS-H 


Magnified  view  of 
shaded  half  star  point 
of  1.  The  dots  indicate 
points  at  which  cartesian 
and  polar  coordinates  are 
computed  for  the  Point  ' 
Periphery  Tape  output. 


FIOTOE  84.  LOCATION  WITHIN  A  SOLID  PROPELLANT  ROCKET  MOTOR  OF  THE  ^PERIPHERY  POINTS. 


±__a _ 

Longitudinal  cross  section  of  rocket  motor  which  reveals  the  normal  H  and 
the  angle  Z- H 


Axial  cross  section  of  one  half  star  point  which  reveals  the  normal !?  and  the 
angle  S-Y 


The  angle  S-H  is  simply  the  angle  between  normals  H  and  S.  Reference  back  to 
Figure  34  will  show  a  three  dimensional  illustration  of  both  f  and  i?  and  the  angle 
S-H  ' 


FIGURE  36.  DEFINITION  OF  NORMALS  S*ANDll 
TOGETHER  WITH  THEIR  ASSOCIATED  ANGLES 
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of  propellant  weight,  etc.  Case  weight  and  payload  weight,  If  specified  to  the 
program,  are  used  in  conjunction  with  the  computed  propellant  center  of  gravity 
travel  to  compute  the  C.  G.  travel  of  the  entire  system. 

e.  Nonunlform  Temperature  Example 

The  importance  of  considering  the  effect  of  a  nonuniform  grain 
temperature  may,  perhaps,  be  shown  best  through  an  example  problem.  The  rooket 
motor  chosen  for  consideration  is  of  the  "broken  back"  star  classification.  A  cross 
section  of  one  half  star  point  of  this  rocket  motor  together  with  its  (relatively  severe) 
thermal  gradient  is  shown  in  Figure  36. 

This  rocket  motor  was  evaluated,  both  graphically  and  by  program 
AGDA,  for  two  temperature  conditions: 

1)  Uniform  grain  temperature 

2)  Variable  grain  temperature 

Computer  output  from  these  two  temperature  conditions  is  shown  in  Tables  XII  and 
XIII.  The  graph  shown  in  Figure  37  compares  the  computer  output  with  that  obtained 
graphically.  As  may  be  seen  from  the  figure,  this  comparison  is  quite  good.  The 
greater  percentage  of  difference  for  the  two  calculational  techniques  in  the  variable 
temperature  case  is  to  be  expected.  Any  small  errors  associated  with  the  graphical 
analysis  of  the  nonuniform  temperature  case  are  cumulative,  while  those  obtained 
from  the  uniform  temperature  case  are  not. 

Again  with  reference  to  Figure  37,  there  is  a  considerable  difference 
in  the  shape  of  the  surface  web  curve  for  the  uniform  versus  nonuniform  grain 
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TNIUKOL-ELKTON 

ADVANCED  GRAIN  DESIGN  COMPUTER  PROGRAM 

STRAIGHT  THRU  SECTION  BROKEN  BACK  STAR  VARIABLE  TEMP  OPTION  DISABLED 

THEORETICAL  ANALYSIS  BY  P.E.MOCRE 
PROGRAMMING  ANALYSIS  BY  D.H. FREDERICK 


WtB 

BURNING 

CHAMBER 

PRUP.  PROP. 

PROPELLANT 

INERTIA 

TENSOR 

BURN 

SURFACE 

VOLUME 

VOLUME 

CG 

i-XX 

l-YY 

1-22 

PCT 

1 SC-INl 

(CU-IN) 

(CU-IN) 

(IN) 

- t$Q- 

FT+CU-IN) 

0 

2634.05 

6556.54 

1737C.08 

7.850 

28072.6 

28072.6 

36323.3 

5 

2666.59 

7430.92 

16495.70 

7.850 

27138.3 

27138.3 

35452.4 

10 

2620.76 

8310.01 

15616.61 

7.850 

26163.4 

26163.4 

34503.8 

15 

2510.06 

9159.67 

14766.93 

7.850  ‘ 

25173«j5 

25173.3 

13495.6 

20 

2445.54 

9978.20 

13948.42 

7.850 

24171.6 

24171.6 

32423.8 

25 

2426.25 

10781.46 

13145.16 

7.850 

23145.0 

23145.1 

31269.4 

30 

2409.74 

11579.18 

12347.44 

7.850 

22085.2 

22085.2 

30080.0 

35 

2393.22 

12371.46 

11555.16 

7.S50 

20991.3 

2099KS 

28796.6 

40 

2376.71 

13158.30 

10768.33 

7.850 

19863.8 

19863.8 

27439.3 

45 

2360.20 

13939.68 

9986.94 

7.830 

18701.9 

18701.9 

26007.2 

50 

2343.68 

14715.61 

9211.01 

7.850 

17505.7 

17505.7 

24500.3 

55 

2327.17 

15486.10 

8440.52 

7.850 

16275.3 

16275.3 

22918.7 

60 

2310.65 

16251.14 

7673.48 

7.850 

15010,8. 

15010.8 

21262.7 

65 

2294.14 

17010.73 

6915.89 

7.850 

13712.4 

13712.4 

19532.7 

70 

2277.62 

17764.87 

6161.75 

7.850 

12360.5 

12380.3 

17729.5 

75 

2261.11 

18513.56 

5413.06 

7.830 

11015.4 

11016.4 

15833.7 

BO 

2247.56 

19256.12 

4670.31 

7.850 

9616.9 

9618.9 

13908.0 

85 

2Z43.14 

19995.82 

3930.80 

7.850 

8186.2 

8166.2 

11886*7 

TO 

2248.46 

20735.90 

3190.72 

7.850 

6712.1 

6712.1 

9783.1 

95 

2260.42 

21478.62 

2448.0C 

7.650 

5192.4 

5192.4 

7591.2 

100 

2279.00 

22225.26 

1701.36 

7.850 

3624.2 

3624.2 

5306.9 

MAX 

2666.59 

22225.26 

17370.08 

7.850 

28072.6 

28072.6 

36323.3 

MIN 

2243.14 

6556.54 

1701.36 

7.850 

3624.2 

3624.2 

5306.9 

AVG 

2378.86 

14604.72 

9321.90 

7.850 

16892.4 

16892.4 

23147.1 

100+ 

2276.10 

105 

1350.28 

22743.75 

1182.87 

2535.6 

2535.6 

no 

1052.03 

23136.73 

789.89 

1706.1 

1706.1 

2510.8 

115 

803.97 

23438.12 

468.51 

1063.3 

1363.3 

1569.1 

120 

334.86 

23666.32 

260.30 

571.3 

571.3 

845*6 

125 

375.92 

23821.76 

104.86 

232.1 

232.1 

344.5 

130 

137.53 

23908.10 

18.52 

3,3 

0.5 

0.7 

PER 

CENT  PROPELLANT  I  0- 

100  PCT  WEB ) 

BY  VOLUME 

*  65.49 

PER 

CENT  PROPELLANT  1100- 

130  PCT  NEB) 

BY  VOLUME 

»  7.11 

TABLE  XII.  FINAL  PAGE  OF  OUTPUT,  TEST  CASE  1 
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Straight  Thru  Section  Broken  Baok  Star  Variable  Temp  Option  Employed  With 

Grain  Temp  Shown  in  Figure  38 
THEORETICAL  ANALYSIS  BY  P.B.MOORE 
PROGRAMMING  ANALYSIS  BY  D.H.PAEOERICK 


MEB 

BURNING 

CHAMBER 

PRLP. 

PROP. 

PROPELLANT 

INERTIA  TENSOR 

BURN 

SURPACE 

VOLUME 

VOLUME 

CO 

1-XX 

l-YY 

i-li 

PCT 

ISQ-IN1 

ICC- IN) 

JCU-IN) 

t  IN) 

- ISO- 

Pt»CU-lN)-* 

0  ~ 

2634,03" 

6553,79 

17370.87 

7.850 

27254*5 

27254.5 

34486*2 

5 

2690,48 

7314.03 

16612.59 

7,850 

26418,0 

26418*0 

1)178.5 

10 

272B.I9 

8090.19 

158J6.43 

7.850 

29529.7 

25529.7 

32987.9 

15 

2737.88 

8864.80 

13061*82 

7*850 

24603*9 

24603.9 

32019.9 

20 

2694.69 

9646.74 

14279*88 

7*850 

23625.4 

23623*4 

30955.) 

.  25-  - 

_ 2111*2.4. 

—10637*38 

13489.24  . 

--1*150. 

..  .22190* 5_ 

2219CU5__ 

-29717.6.  . 

30 

2794.84 

11239.83 

12686.80 

7.850 

21500.9 

21500.9 

28524*2 

35 

2856.82 

12084.30 

11842.12 

7*  BOO 

20322.6 

20322.6 

27131*5 

40 

2904.70 

12976.34 

10950.28 

7.850 

19036.4 

19036.4 

2557ft* 8 

45 

2943.73 

13914.81 

10011*82 

7.830 

17641.5 

17641.5 

23837*9 

50 

2993.13 

14895.37 

9031*05 

7.850 

16138.2 

16138.2 

21970*5 

-  .33 _ 1034*33 

— 13910*11 . 

7996.23 

.  ,7*.850. 

.  14503*6 . 

.  14503*6 _ 198(2,1  . 

40 

3078*53 

17011.65 

6914*97 

7*850 

12740.9 

12740.9 

17590*8 

43 

3037.88 

18150.32 

5776*30 

7*850 

10833.7 

15833.7 

19075*8 

70 

2812.41 

19301.47 

4623.15 

7*850 

8832.5 

8832.9 

12387*1 

73 

2557.47 

20403.33 

3523.09 

7*850 

6822*3 

6822.3 

9624.3 

•0 

2298.65 

21461.11 

2465.51 

7,850 

4812.5 

4812.5 

8811.6 

_ as _ 

_ 1317.69 

— 22210*76- 

.  ..  -1635*86 

..  .7*850. 

- 3257*3.  - 

—  325743 — 

-.4621*0  . 

90 

1151.91 

22844.63 

1082.00 

7,850 

2147.5 

2147.9 

3060.2 

95 

834.19 

23273.67 

647*95 

7*850 

1298.0 

1298.0 

1856*6 

IOC 

624.34 

23595.26 

531*36 

7.850 

669.9 

669.9 

96147 

MAX 

3078.53 

23395*26 

17370,87 

7*850 

27254.5 

27234*!! 

3463^.2 

_Mia — 

_ 624*14 _ 6553*23 . 

-.311*36- 

.  .7*830 

-  .  669*9. .. 

-669*9 _ 

-  .((IaT.  . 

AVG 

2461.20 

13250.84 

8675.77 

7*850 

14789.5 

14789.9 

19678*6 

1004 

606.84 

103 

377.39 

23811.58 

1  .04 

234.5 

234.5 

337*8 

lie 

103.60 

23916.80 

9.82 

20.1 

20.1  * 

29*1 

1*9 

-J3*00  . 

24976.6? 

...  a.oo 

. 0*G  ... 

.  ..0*0.....: 

. JDcO 

120 

0.00 

23926.62 

0.00 

0.0 

0.0 

0.0 

123 

0.00 

23926.62 

0.00 

0.0 

0.0 

0.0 

130 

0.00 

23926.62 

0.00 

0.0 

0.0 

0.0 

PER  CENT  PROPELLANT  I  0-100  PCT  NE«I  BY  VOLUME  »  71.22 
PER- CENT  PROPILL'anT”  i- 100-1 30  PCTWEB)  BY  VOLUME  ■  1*38 
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temperature  cases.  It  is  conceivable  that  everything  from  a  "slightly  off  target"  to  a 
major  catastrophe  could  result  from  firing  such  a  rocket  motor  with  severe  thermal 
gradient.  These  effects  would  be  especially  pronounced  if  the  propellant  selected  had 
a  higher  than  average  temperature  sensitivity  coefficient  (TTj^. 


Iir  EFFECT  OF  STRAIN  ON  BURNING  RATE 

The  effect  of  strain  on  propellant  structural  Integrity  has  received  considerable 

i ' 

attention  in  the  past.  However,  only  limited  effort  has  been  directed  toward 
determining  whether  a  propellant  under  strain  will  burn  predictably  and  the  extent 
that  burning  rate  will  vary  as  a  function  of  the  induced  strain. 

Coy,  using  an  ammonium  perchlorate-polyurethane  propellant,  measured 
approximately  1. 5  percent  increase  in  burning  rate  for  each  percent  strain.  (Coy 
concluded  that  this  strain-burning  rate  relationship  was  produced  by  the  formation 
of  a  liquid  phase  at  the  propellant  burning  surface. ) 

Other  tests  on  plastisol  propellants  which  were  conducted  at  Atlantic  Research 
Corporation*  showed  a  definite  Increase  in  burning  rate  but  at  much  higher  strain 
levels  rv  .  hose  used  by  Coy  (on  the  order  of  100  percent).  In  contrast  to  the  work 
of  Coy  ai.  *  ARC,  the  Longhorn  Division*  of  the  Thiokol  Chemical  Corporation  static 
fired  a  Nike-Hercules  motor  that  had  experienced  a  severe  slump  problem.  The 
slump  (time-dependent  strain  due  to  gravity  force)  had  progressed  to  the  extent  that 
the  opposite  star  points  in  the  grain  almost  joined.  When  the  motor  was  static 
tested ,  the  measured  burning  rate  was  no  different  from  that  measured  in  laboratory 
control  tests.  Therefore,  it  was  concluded  that,  at  least  for  this  particular 
formulation,  burning  rate  was  unaffected  by  strain. 


♦Private  communication  with  Cortland  Robinson  of  Atlantic  Research  Corporation 
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In  these  studies,  eaoh  of  which  dealt  with  a  single  type  of  propellant,  no  general 
oriteria  for  all  propellants  could  be  established.  It  was  the  main  purpore  of  this  study 
to  determine  whether  such  criteria  existed  and,  if  possible,  to  determine  the  means  by 
which  they  could  be  established. 

TECHNICAL  APPROACH 
1.  Propellants 

The  first  propellant  investigated  was  an  HA/MAPO  hydrocarbon  propellant 
(designated  Hydrocarbon  "A"  in  this  report).  Depending  on  the  results  obtained  with 
this  propellant,  other  propellants  or  propellant  systems  were  to  be  considered  for 
evaluation.  Before  the  end  of  the  program,  three  other  propellants  were  studied:  HB/ 
epoxy  (Hydrocarbon  "B"),  a  polyurethane,  and  a  plastisol. 

The  standard  JANAF  physical  properties  of  those  propellants  tested  are 

given  below; 

JANAF  Maximum  Stress,  Strain  at  Maximum 

Propellant  Type  Modulus.  e3.2.psl»* **  psi _  Stress,  c3.2.tn/ln*» 


Hydrocarbon"  A" 

590 

100 

0.24 

Polyurethane 

545 

77 

0.25 

Hydrocarbon  "B"  (Mod.)* 

1565 

113 

0.17 

Plastisol 

400 

90 

0.20 

*This  material  was  purposely  formulated  at  a  high  epoxy -to-carboxy  ratio  in  order 
to  achieve  a  higher  modulus  material.  The  normal  values  for  this  formulation  are 
600  psl,  80  psi,  and  0-25  in. /in. ,  respectively. 

**The  value  3.2  refers  to  the  effective  gauge  length  used  in  the  determination  of  the 
physical  properties  shown. 
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2.  Tests  and  Procedures 


T^o  distinct  areas  were  to  be  Investigated  at  the  beginning  of  this  program . 
One  area  was  concerned  with  measuring  the  change  in  burning  rate  of  specimens  with 
induced  uniform  strain  fields  ranging  from  0  to  12  1/2  percent  at  pressures  of  500,800, 
and  1200  psi.  However,  the  strain  levels  induced  across  the  face  of  a  grain  by  firing  or 
thermal  sources  are  nonuniform  because  of  the  various  geometries  encountered  in  solid 
propellant  grain  designs .  Therefore,  a  second  type  of  test  in  which  a  nonuniform  strain 
distribution  was  produced  in  a  propellant  sample  was  planned.  The  burning  rate  of  this 
type  of  sample  was  to  be  determined  at  various  positions  within  the  sample  under 
different  strain  concentrations  and  compared  with  the  uniform  strain  data  for  possible 
correlation  with  actual  motor  firings .  The  uniform  strain  field  teats  were  given  the 
greater  initial  attention  since  it  was  felt  that  the  results  obtained  would  elucidate  some 
of  the  experimental  problems  that  might  be  encountered  in  the  more  complex  nonuniform 
strain  field  tests . 

When  the  uniform  strain  field  tests  on  hydrocarbon  "A"  were  completed,  it 
was  found  that  virtually  no  difference  existed  between  the  linear  burning  rate  of  the  un¬ 
strained  strands  and  those  strained  up  to  12  1/2  percent.  In  the  light  of  the  conflict  that 
existed  between  these  results  and  those  of  Coy  and  ARC,  it  was  conjectured  that  strain 
does  not  affect  burning  rate  directly  but  instead  affects  a  material  property  which,  in 
turn, can  be  related  to  changes  in  burning  rate.  A  quantitative  parameter,  which  might 
bo  directly  related  to  burning  rate,  is  the  change  In  /olume  with  strain.  This  property, 
referred  to  here  as  Poisson's  ratio,  is  defined  as  the  ratio  of  the  lateral  strain,  <33,  to 

the  longitudinal  strain,  e^,  in  a  material  subjected  to  a  uniaxial  tensile  stress. 
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When  a  material  is  stretched,  the  cross-sectional  area  changes  as  well  as 

length.  Poisson's  ratio,  v,  is  the  parameter  relating  tnese  changes  in  dimension 

C2  _  Percent  Change  in  Width  &W/Wp 
«T  Percent  Change  in  Length  AL/Lo 

If  a  material  is  incompressible,  an  elongation  would  result  in  changing. the 
cross-sectional  area  but  not  the  volume.  Since  linear  burning  rate  is  independent  of 
cross-sectional  area,  no  change  in  burning  rate  would  be  expected.  On  the  other  hand, 
if  the  material  is  compressible  and  incurred  a  change  in  volume  as  a  result  of  elongation, 
the  resulting  change  in  density  should  cause  a  change  in  the  burning  rate. 

It  can  be  shown  that  if  the  volume  of  a  material  remains  constant  when 
subjected  to  very  small  strains,  v  is  a  constant  and  equals  0.50.  Generally,  materials 
increase  in  volume  when  subjected  to  a  tensile  strain,  resulting  in  a  v  less  than  0-50. 

For  most  materials  Jlies  between  0.2  and  0-5  and  approaches  0.50  for  rubbers  or 
liquids . 

When  a  material  is  subjected  to  a  strain  sufficiently  large  to  prohibit  the 
use  of  infinitesimal  strain  theory,  1  is  no  longer  a  constant  but  a  function  of  the  strain: 


v  .  1 
a  -  l 


l 


rl/2  -1/2 


a 


where  v 


v  o/ 
Poisson's  ratio 


a  -  Principal  extension  ratio  (1  +  e^) 

V  -  Volume 


VQ  -  Original  volume 


(2.2) 
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If  the  material  is  inoompressible,  t.e.,  V/V0  a 


1 


i  -  a  + «,) 


-1/2 


1,  equation  (2.2)  reduoos  to: 


where  v/v^  ?"l 

It  can  be  seen  from  equation  (2.2)  that  the  lower  the  value  of  v,  the  greater  the  volume 
inorease  and  hence  the  greater  the  density  ohange  due  to  the  strain,  Therefore,  It 
can  be  expected  that  the  burning  rate  of  the  material  which  displays  the  greatest  depar¬ 
ture  from  Incompressibility  would  be  most  affected  by  strain.  The  dimensionless  ratio 
v/V(,  henceforth  referred  to  as  the  "volumetric  change  index,"  is  representative  of  the 
departure  from  incompressibility.  To  teat  this  hypothesis  the  linear  burning  rates  of 
the  four  candidate  materials  were  determined  over  the  pressure  ranges  mentioned  above. 
In  addition,  the  volumetric  change  Index  was  determined  as  a  function  of  strain  for  the 
same  materials.  To  accomplish  this  an  apparatus  for  the  experimental  determination 
of  vwas  developed,  based  on  fundamental  concepts  previously  devised  under  corporate 
funding. 

A  discussion  of  the  methods  and  equipment  used  for  the  three  procedures 
and  a  discussion  of  the  results  arc  given  below. 

a.  Uniform  Strain  Field  Tests 

The  burning  of  the  propellant  was  accomplished  In  a  standard  Auto¬ 
clave  Engineering  Company  strand  burner  (see  Figure  38).  The  strand  burner  consisted 
of  a  3-lnch-dlameter  by  12-lnch-deep  cylindrical  burning  chamber.  The  sample  holder 
was  lowered  into  this  chamber.  In  addition  to  supporting  the  strand,  the  sample  holder 
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provided  the  means  for  igniting  the  propellant.  The  burner  was  also  equipped  with 
three  clocks:  one  clock  measured  the  time  to  burn  through  the  first  half  of  the  strand: 
another  measured  the  time  to  burn  through  the  last  half;  and  the  third  clock  measured 
total  burning  time  and  provided  a  check  on  the  other  two.  The  bomb  was  pressurized 
with  nitrogen  by  a  reciprocating  booster  pump.  The  maximum  allowable  pressure  for 
the  bomb  was  10, 000  psia.  Any  pressure  buildup  in  the  bomb  during  burning  was  ab¬ 
sorbed  by  a  1750-milliliter  surge  tank.  Burning  rates  can  be  measured  to  an  accuracy 
of  +  .005  in. /min. 

The  specimen  consisted  of  a  1/4-inoh  by  1/4-inch  by  6-inch  strand  of 
propellant.  The  specimen  ends  were  inserted  through  two  Teflon  collars  and  bonded 
with  Armstrong  A-l  cement.  In  this  form  it  was  similar  to  the  end  bonded  tensile 
specimen®  developed  at  the  Thiokol-Redstone  Division.  Four  constant  strain  devices 
were  revised  to  provide  a  means  for  achieving  the  desired  strain  levels  in  the  strand. 
Each  device  was  capable  of  independently  straining  two  strands.  Once  the  desired 
strain  level  was  achieved,  the  sample  was  coated  on  all  sides  with  Armstrong  A-l 
cement.  This  procedure  is  shown  in  Figure  39.  The  cured  adhesive  bonded  to  the 
strandjand  since  the  coating  was  stiff,  the  strain  level  in  the  grain  was  maintained.  In 
addition,  the  cement  acted  as  an  Inhibitor  to  prevent  side  burning  during  firing. 

Three  wires  were  placed  through  the  strands  spaced  2-1/2  inches 

apart  in  the  direction  of  burning.  These  wires  were  connected  to  the  three  clocks  of  the 

strand  burner.  The  circuitry  was  designed  so  that  as  the  propellant  burned  and  the 

wires  were  severed,  they  started  and  stopped  the  timers,  causing  the  burning  time 
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between  each  wire  to  be  recorded.  One  additional  wire  was  Inserted  through  the  propel¬ 
lant  near  the  top  of  the  strand,  which  acted  as  an  igniter . 
b.  Nonuniform  Strain  Field  Test 

The  apparatus  necessary  to  induce  a  nonuniform  strain  field  into  solid 
propellant  sample  is  shown  in  Figure  40.  A  triangular  sample  of  propellant  2  inches  by 
2  inches  by  2  inches  by  1  inch  thick  was  cast  into  a  vise  made  of  1-inch  thiok  Plexiglas 
to  make  a  direct  bond  between  vise  wall  and  the  sample.  When  the  sample  cured,  its 
face  was  covered  with  a  birefringent  coating  1/8  inch  thick  with  a  thin  film  of  reflective 
paint  at  the  interface  between  the  coating  and  the  propellant.  In  this  manner,  the  strain 
levels  in  the  sample  could  be  determined  by  PhotoStress  (reflective  photoelasticity 
technique;.  When  the  strain  geography  in  the  sample  for  a  particular  amount  of  vise 
wall  displacement  was  known,  a  new  sample  could  be  similarly  prepared  without  the 
birefringent  coating  and  paint.  The  sample  might  then  be  compressed  by  the  vise  to 
achieve  the  previously  determined  strain  pattern.  Ar.  igniter  wire  can  be  inserted  in  the 
1/4-inch-diameter  groove  at  the  lower  portion  of  the  vise. 

The  bomb  assembly  (which  has  been  utilized  previously  to  study  the 
effect  of  defects  on  burning)  consisted  of  two  basic  parts:  a  windowed  container  and  an 
Internal  assembly  to  hold  the  sample  (see  Figures  41  and  42).  The  bomb  has  a  window 
and  a  vented  stainless  steel  cylinder.  Two  1-inch-thick  by  1/2-inch-wide  by  5-inches- 
long  Herculite  windows  held  in  place  with  steel  plate  and  screws  were  located  on  opposite 
sides  of  the  bomb.  An  outlet  through  the  side  of  the  bomb  near  the  bottom  provided  a 
means  for  pressurizing  the  bomb.  Two  other  connections  provided  a  means  for  pressure 

gauge  connection  and  a  safety  blow-out  disk. 
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FIGURE  40.  TEST  JIG  FOR  NON-UNIFORM  STRAIN  VERSUS  BURNING  RATE  TEST 
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FIGURE  41.  WINDOW  BOMB  INTERNAL  ASSEMBLY 


The  speoimen  can  be  mounted  on  the  Internal  assembly  sc  that  it  can 
be  viewed  through  the  bomb  window .  The  system  can  then  be  sealed  and  pressurized 
with  nitrogen  to  the  desired  pressure  level  and  ignited. 

A  Wollensak  Fastax  16mm  high  speed  camera  can  be  used  to  photograpl 
the  surface  of  the  propellant  as  it  burns.  The  propellant  should  provide  sufficient 
illumination  to  permit  good  resolution;  however,  a  standard  photo  flash  unit  can  be  U3ed 
to  obtain  a  picture  of  the  specimen  before  ignition.  Timing  reference  marks  super¬ 
imposed  on  the  film  can  be  used  to  provide  for  burning  rate  measurements, 
o.  Poisson's  Ratio  Test 

An  apparatus  for  idle  determination  of  the  volumetric  changes  of  pro¬ 
pellant  has  been  developed  and  consists  of  a  point  light  source,  an  end-bonded  cylindrical 

g 

tensile  specimen, .  a  10X  lens,  an  Instron  tensile  tester,  and  a  sheet  if  graph  paper. 

This  is  shown  schematically  in  Figure  43.  A  specimen  of  known  diameter  was  placed 
in  the  path  of  the  point  light  source.  The  image  was  magnified  ten  times  by  means  of 
the  lens  and  projected  on  the  graph  paper.  The  highly  magnified  diameter  was  measurec 
and  the  volumetric  changes  are  calculated  from  the  measured  lateral  and  longitudinal 
strams,  using  equation  2.2.  It  should  be  mentioned  that  since  the  optics  of  this  device 
are  still  in  the  preliminary  stage  of  development,  it  was  necessary  that  the  average  of 
at  least  five  tests  were  reported.  The  preparation  of  samples  must  be  performed  with 
extreme  care  to  assure  uniaxial  stress  conditions  during  deformation. 
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d.  Discussion  of  Results 


The  results  of  the  linear  burning  rate  versus  strain  tests  are  given  in 
Tables  XIV  through  XXIV  and  Figure  44  for  the  four  propellants  tested.  It  can  be  seen 
from  the  graph  that  the  influence  of  strain  on  linear  burning  rate  is  different  for  the  four 
formulations.  The  data  for  hydrocarbon  "A"  show  virtually  no  difference  between 
burning  rates  of  the  unstrained  strands  and  those  strained  up  to  12  1/2  percent.  This 
independence  of  strain  on  burning  rate  for  this  particular  propellant  supports  the  Long¬ 
horn  Division*  results  on  the  slumped  Nike -Hercules  motor. 

It  was  originally  intended  to  obtain  data  at  0,  4,  8  and  16  percent 
strain  and  500,  900  and  1200  psi  burning  test  pressure.  It  became  increasingly  difficult 
to  achieve  16-percent  strain  in  the  type  of  specimen  used  without  sample  rupture.  For 
hydrocarbon  "A"  the  maximum  strain  achieved  was  12  1/2  percent.  For  the  other 
materials  even  this  was  not  possible,  and  10-percent  strain  was  used.  In  addition, 
since  virtually  all  propellants  experience  a  linear,  log-log,  pressure-burning  rate  trace, 
it  was  decided  to  use  500  and  1500  psi  rather  than  500,  900  and  1200  psi.  This  was  also 
done  to  permit  a  wider  spreau  in  the  data.  Therefore,  some  of  the  results  are  reported 
at  different  strains  and  pressures. 

The  data  for  the  polyurethane  propellant  show  an  increase  of  4.0 
percent  in  burning  rate  at  1000  psi  and  0.10  in.  /in.  strain.  The  burning  rate  of  the 
modified  hydrocarbon  "B"  increased  9.1  percent  at  0.10  in. /in.  strain  and  1500  psi, 
and  the  plastisol  exhibited  an  increase  in  burning  rate  of  27  percent  at  0.10  in. /in. 


and  1500  psi. 
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TABLE  XV 

BURNING  RATE  OF  HYDROCARBON  "A11  (4%  STRAIN) 


Average  Pressure, 


Test  Number  oai 

,  ^ — p-r —  .  i  K — 

1 

2 

3 

4 
'6 
8 

Aver. 

1 

'  '  2 
.  3 
'  '4 

5 

6 

;  7 

Aver. 

.  1 

2 

3 

.  4 
,  5 
6 
7 

Aver. 


Burning  Rate,  in/sec 


588 

0.416 

593 

0.404 

563 

0.414 

550 

0.391 

565 

0.404 

565 

0.  389 

~570 

Aver. 

0.403 

1008 

0.508 

1015 

0.594 

1010 

0.541 

1008 

0.511 

1001 

0. 508 

1000 

0. 589 

1000 

0.  595 

1006 

Aver. 

0.  549 

1313 

0.576 

1318 

0.569 

1313 

0.  556 

1325 

0.  587 

1325 

0.588 

1325 

0.587 

1325 

0.582 

1320 

Aver. 
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TABLE  XVI 


«* 


BURNING  RATE  OF  HYDROCARBON  "A"  (8%  STRAII 


Test  Number 


Average  Pressure, 
psi 


Burning  Rate,  ln/sec 


1 

2 

3 

4 

5 
8 

7 

8 
9 


575 

0.  388 

563 

0.407 

563 

0.  394 

563 

0.415 

563 

0.404 

563 

0. 394 

550 

0.  396 

563 

0.410 

563 

0.411 

Aver.  863 


Aver.  0.402 


1 

1120 

0.568 

2 

998 

0.  585 

3 

1000 

0. 529 

4 

1120 

0.542 

5 

1120 

0.575 

6 

1000 

0.494 

7 

1000 

0.483 

8 

1000 

0.497 

Aver. 

1045 

Aver.  0. 532 

1 

1313 

0.  563 

2 

1313 

0.553 

3 

1338 

0.587 

4 

1313 

0,578 

5 

1318 

0^  558 

6 

1313 

0.550 

Aver. 


1318 


Aver. 


0.565 


TABLE  XVII 


BURNING  RATE  OF  HYDROCARBON  "A”  <12-1/2%  STRAIN! 


Te*t  Number 


Average  Pressure, 
_ gsl _ 


Burning  Rate,  tn/seo 


1  . 

575 

■> 

0.404, 

2 

575 

0.442 

3 

575 

0.401 

4 

563 

0.401 

5 

563 

0.404 

6 

587 

0. 407  * 

.7 

575 

0.390 

8 

562 

0.393 

Aver. 

572.> 

Aver. 

0.405 

1 

975 

0.488 

2 

963 

0.478 

3 

950 

0.475 

4 

950 

0.472 

5 

950 

0.  521 

.6 

950 

0.514 

7 

963 

0. 529 

8 

975 

Q.  497 

Aver. 

960 

Aver. 

0. 484 

* 

X 

! 

t 

1300 

0.  586 

2 

1300 

0.545 

3 

1300 

0.  535 

4 

1300 

0.534 

6 

1300 

0.  555 

6 

1300 

0.  547 

7 

1300 

0.539 

3 

1300 

0.562 

Aver. 

1300 

Aver. 

0.  550 
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TABLE  XVIII 


1  575  0.160 

550  0.157 
538  0.157 
538  0157 
500  0.158 


Aver.  0.1578 


1075  0. 169 

1075  0.131 

1088  0  183 

1075  0.*8.i 

10  1063  0.  if  £ 

Aver.  (TUTe 
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TABLE  XDC 


BURNING  RATE  OF  POLYURETHANE  PROPELLANT  (10%  STRAIN) 


t 

Test  Number 

Average  Pressure* 
pal 

s 

Burning  Rate,  in./seo 

.,1  . 

538 

*  1 

0  460 

2 

530 

0.159 

3  • 

530 

0.161 

Aver.  0.160 

5  .  •  ... 

1068 

0.166 

6 

1050 

0.170 

7  ■■  .  . 

1075 

0.168 

8 

1063 

0.167 

Aver.  0.168 
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TABLE  XX 


LINEAR  BURNING  RATE  OF  HYDROCARBON  ”B”  (MOD)  (0%  STRAIN) 


Test  Number 


X 

2 


3 

4 

5 


Average  Pressure* 
pal 

S6S 

SSO 

1565 

1560 

1565 


Burning  Rate,  in./seo 

0.246 

0.247 

Aver.  0.247 

0.396 

0.396 

0.396 

Aver.  0.396 
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TABLE  XXI 


T INE  ab  MTBtnwn  RATE  OF  HYDROCARBON  »B”  (MOD)  (6%  STRAIN) 


Test  Number 
1 


Average  Pressure , 
psl  _ 


563 

550 

550 

550 


1550 

1800 

1800 

1575 


Burning  Rate  >  tn./aeo 

0.264 
0.256 
0.265  . 
0.256 


Aver.  0.258 

0.409 

0.411 

0.412 

0.415 


Aver.  0.412  . 
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TABLE  XXn 


LINEAR  BURNING  RATE  OF  HYDROCARBON  "B11  (MOD)  (10%  STRAIN) 


Test  Number 


Average  Pressure » 
_ psl _ 


Burning  Rate,  tn./seo 


1  650 

675 
675 
650 


0.269 

0.287 

0.285 

0.274 


Aver.  0.281 


5 

6 
7 


1600 

1590 

1595 


0.418 

0.437 

0.443 


Aver.  0.432 
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TABLE  XXIII 


BURNING  RATE  OF  PLASTISOL 


Test  Number 


Average  Pressure , 
si  '  ' 


Rate,  in. 


-122- 


Aver. 


0.32.1 
0 !  32^ 
0.317 
0.323 
0.317 
0.320 


Aver. 


0.320 
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TABLE  XXIV 


1  '*®**rt**XH«  i  Wes**!  wnsfps 


It  oan  be  seen  that  the  ourves  for  the  strained  oondltion  are  not 
neoessarily  parallel  to  their  respective  unstrained  curves,  thus  indicating  a  different 
pressure  dependence  for  each  case.  It  appears  that  strain  may  increase  the  pressure 
dependence,  although  the  data  are  inconclusive. 

Consideration  was  given  to  the  "mass  burning  rate"  in  addition  to  the 
linear  burning  rate.  When  a  strain  is  imposed  on  a  strand,  it  will  permit  an  increase 
in  time  necessary  to  burn  the  propellant  completely  (l.e. ,  a  10-percent  strain  will 
produce  a  10-percent  longer  burning  time  for  an  incompressible  material).  The  mass 
burning  rate  is  thus  reduced  by  the  time  increase  required  to  consume  the  sample. 
Therefore,  mass  burning  rate  is  a  direct  function  of  the  reciprocal  of  the  strain  for 
the  incompressible  case.  The  mass  burning  rate  calculations  were  based  on  the 
following  relationship: 

w  =  rAp 

where  w  ~  Mass  burning  rate,  lb/sec 

r  =*  Linear  burning  rate  (taken  from  experimental  data), 
in .  /sec 

2 

A  =  Burning  surface  area,  in. 

3 

P  =  Material  density,  lb/in. 

The  value  of  A  changes  with  strain  and  was  calculated  from  the 
Poisson's  ratio  data.  From  this  value  a  new  density  of  the  material  was  computed.  The 
combination  of  the  linear  burning  rate  (r),  area  (A),  and  density  (p)  of  the  material  In 
the  strained  condition  permitted  the  calculations  of  the  mass  burning  rate,  w.. 
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The  computed  mass  burning  rate  shown  in  Tables  XXV  through  XXVIII 
and  Figure  45  experienced  a  deorease  with  strain  for  all  propellants  except  the  plastiacl. 
This  is  not  necessarily  inconsistent  since  it  is  conceivable  that  strain  can  inorease  the 
burning  rate  to  a  degree  which  will  completely  offset  the  time  to  burn  the  additional 
length  of  the  specimen. 

It  should  be  mentioned  here  that  the  burning  rate  calculations  were 
based  on  three  assumptions: 

1)  Perfect  homogen°ity  exists  throughout  the  strained 
sample.  This  would  mdicate  that  the  dewetting  action 
due  to  strain  was  also  uniform  along  the  length  of  the 
sample.  This  is  probably  not  true,  and  localized 
areas  of  dewetting  probably  exist.  This  assumption 
affects  the  calculation  of  the  density  in  the  deformed 
state,  and  the  linear  burning  rate  is  not  really  constant 
throughout  the  length  of  the  strand.  However,  the 
linear  rates  measured  include  the  gross  summation  of 
such  areas,  and  the  results  should  approximate  the 
actual  behavior . 

2)  The  assumption  is  also  made  that  the  dimensions  of  the 
strand  are  1/4  inch  by  1/4  inch  along  the  entire  length 
and  that  after  straining  the  new  cross  section  is  uniform 
also.  Where  this  nonuniformity  does  not  enter  into  the 
determinations  of  linear  burning  rate,  it  does  affect  the 
mass  burning  rate  calculations. 

3)  The  mass  of  the  inhibitor  does  not  influence  the  mass 
burning  rate.  This  assumption  w as  verified  in  the 
linear  burning  rate  studies  by  comparison  with  control 
data  on  uninhibited  strands.  The  difference  was  felt 
to  be  negligible. 

These  factors  in  conjunction  with  the  nature  of  the  plastisol  propel¬ 
lant,  limit  the  reliability  of  mass  burning  rate  results  obtained  In  this  indirect  manner. 
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MASS  BURNING  RATE  VERSUS  STRAIN.  HYDROCARBON  "A" 
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MASS  BURNING  RATE  VERSUS  STRAIN,  POLYURETHANE  PROPELLANT 


£ 


to 

a 


CD 

a 


s  *s 

«  QJ 
& 


a  »\ 


Im 

H  , 


cti 

gjN 


*? 


O 

p-t 

H 


*? 


to 


CO 

a 


to 

a 

« 


in 

co 

to 

o 


*? 


o 

H 

<8 

CO 


? 

o 

f-H 

tg 

OS 


w 

I 

o 


3  M 

£  <§ 

U)  CO 


o 

o 

to 


*? 


o 

© 

to 


to 

a 

« 


o 

£> 

to 

o 


■? 


tg 

<?> 


a  s  a  a 


c? 

o 


© 
pH 

M  * 

OS  CO 

9  to 


© 

o 

a 


H 

o 

00 

to 

CO 

CO 

rH 

ft 

» 

• 

• 

o 

o 

o 

o. 

i 

to 

N 

co 

t3» 

os 

to 

to 

© 

S  W 

W  rW 

d 

• 

© 

table  xxvtt 


?  ?  <?,  <J  7  n 


a  s 
S  s 


•El!  c. 

lj  - 


“5  5 

c3  .  © 

II 

^  « 


I 

in  2 

O  -  « 

>  a  • 
•£  o 


5  s 


s  s 


•? 

CO 

1 

*? 

? 

o 

rf 

*? 

o 

fH 

1 

o 

fH 

i 

© 

iH 

1 

O 

3 

09 

1 

O 

tH 

S 

c 

H 

d 

t 

fH 

• 

Q 

5  ° 

o  © 

s  w 


§  § 

UP  UP 


-129- 


2 


n 


CO 


*? 


*? 


2  « 
a 


M 

in 


><1 

in 


a  a 


M 

CO 

• 

to 


H 

CO 

% 

to 


M  «  «  P5 

I  I  I  I 

o  o  o  © 


X  H 

CO  <0 


,8 

eo 


K 

cc 

oo 


Q> 

^  *CQ 

o 

o 

© 

© 

o 

o 

o 

o 

§  °J 

in 

in 

in 

in 

rH 

£ 


8 


3^ 

"a 

o 

CM 

CO 

m 

in 

o 

o 

to 

(n  « 

o 

o 

d 

d 

a  1 

h3  « 

v 

f 

3« 

>  ^ 


in 

CM 

F— I 

CO 

• 

© 


o 

CM 

w 

CO 


00 

in 

o 

• 

o 


kOO 

&  d 

CO 

a 


£ 


<N 

05 

lO 

o 


OQ 

<N 

in 

o 

• 

o 


§  a 

t  a 

o  £ 

&<  CO 


-130- 


A  6533 


M 


a  a  &  a  a 


M 


A  d  U  6  d 


gOt  *  00B/qt  ‘0}By  flujuing  sfltJW  q{  x  oo»/qt  'o^wg  gujuang  sut)^ 

b 


-131- 


FIGURE  45. 


The  different  degrees  that  the  various  materials  were  affected  by 


strain  Initiated  the  effort  to  determine  whether  a  criterion  for  these  changes  in  burning 
rate  could  be  established  based  on  the  change  in  volume  due  to  strain  (Poisson’s  ratio) 
rather  than  strain  "per  ae."  The  value  of  the  expression  for  an  incompressible 
material  (equation  2.3)  was  calculated  for  strain  values  for  the  four  propellants  deter¬ 
mined  experimentally.  The  results  are  shown  in  Figure  46  and  Table  XXIX.  It  can  be 
seen  that  the  departure  of  the  Poisson's  ratio  for  the  individual  propellants  from  the 
incompressible  material  curve  is  in  direct  sequence  with  the  degree  of  influence  of 
strain  on  burning  rate.  The  data  for  hydrocarbon  "A,"  which  showed  virtually  no  effect 
of  strain  on  burning  rate,  follow  quite  closely  the  curve  for  an  incompressible  material, 
Indicating  that  this  particular  formulation  exhibits  a  negligible  change  in  density  due  to 
strain.  On  the  other  hand,  the  plastisol,  which  displayed  the  departure  from  the  incom¬ 
pressible  curve,  also  had  the  most  strain-sensitive  burning  rate.  In  order  to  examine 
this  behavior  more  closely,  the  percent  increase  in  linear  burning  rate,  A'r/r0,  is 
compared  with  the  volumetric  change  index,  v/vi,  at  the  various  pressures  and 
strains  tested.  The  results  are  shown  on  page  135. 
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LINEAR  BURNING  RATE  AND  VOLUMETRIC  CHANGE  DATA 


Percent  Increase  tn 


Strain 

Propellant  (6,)  % 

Pressure,  psi 

Linear  Burning  Rate 
due  to  Strain, 
Lt/to 

Volume  trio  Change 
Index  (vA>j) 
x  100 

Ar/r  0  +  fli’/p 

Hydro¬ 

carbon  "A' 

'  4 

500 

-1.7 

100.0 

101.7 

1320 

0.35 

100.35 

8 

500 

1.95 

99.36 

101.31 

1320 

1.20 

100.56 

121/2 

500 

1.219 

100.7 

101.92 

1320 

2.777 

103.48 

Poly¬ 

urethane 

10 

500 

1.394 

97.42 

98.81 

1000 

3.96 

101.38 

Hydro¬ 

carbon  "B 

"  5 

500 

4.453 

90.25 

94.70 

(Mod) 

1500 

4.040 

94.29 

10 

500 

13.77 

87.96 

101.7 

1500 

9.09 

97.05 

Plastisol 

10 

500 

39.06 

69.89 

109.0 

1500 

26.32 

96.2 

An  interesting  relationship  appears  to  exist  between  the  two  parameters 
at  the  10-percent  strain  level,  i.e-: 

Ar/r0  +  (yA^i)  ^constant  « 100  (3.5) 

Since  only  the  0  and  10-percent  strain  effects  were  obtained  on  all  but 
hydrocarbon  "A, M  it  is  not  possible  at  this  time  to  determine  whether  this  relationship  is 
consistent  at  intermediate  strain  values.  This,  perhaps,  could  be  a  subject  for  future 
study. 
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Dr.  Landel  of  JPL  was  contacted  to  determine  the  value  of  Poisson's 
ratio  at  10-percent  strain  for  the  material  tested  by  Coy.  *  The  value  of  Poisson's  ratio 
for  that  particular  material  was  0.415.  The  rise  in  burning  rate,  Ar/r0,  was  14  percent 
and v/v^ equaled  88  percent;  he  correlation,  at  least  for  linear  burning  rate,  is  thus 
further  substantiated.  From  the  above  table  it  rppears  that  equation  2.4  is  valid  with  an 
accuracy  of +  0.05. 

A  similar  table  for  mass  burning  rate  is  shown  below. 

•  *  MASS  BURNING  RATE  AND  VOLUMETRIC  DATA 


Percent  Change  in 
Mass  Burning  Rate 


Propellant 

Strain 

(«!),% 

Pressure, 

pal 

Due  to  Strain, 

&  w /f  o 

Volumetric  Change 
Index,  (V)Vj)Xl00 

Aw/w0  +^M) 

Hydro¬ 

carbon  "A" 

4 

500 

1320 

-8.259 

-7.480 

100.00 

90.74 

92.51 

8 

500 

1320 

-18.049 

-15.86 

99.36 

83.311 

83.50 

12  1/2 

500 

1320 

-22.22 

-24.67 

100.7 

78.48 

76.03 

Poly- 

urethane 

10 

500 

1000 

-16.214 

-13.816 

97.42 

81.206 

83.604 

Hydro- 

carbon  "B" 

5 

500 

1500 

-5.906 

-6.266 

90.25 

84.34 

83.98 

10 

500 

1500 

-7.564 

-11.434 

87.96 

80.40 

76.53 

Plastisol 

10 

500 

+13.94 

69.89 

83.83 

1500 

+3.53 

73.42 

*Dr ,  Landel  was  contacted  for  this  information  since  Coy  had  since  left  the  employ  of 
JPL.  Dr.  Landel  had  been  consulted  before  to  thl3  as  to  the  advisability  of  pursuing 
the  Poisson's  ratio  approach. 
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Unlike  the  previous  case,  the  expression  A  w'/v  Q  +  £vAh)  hs  &  ooa- 
stant  but  chafes  with  strain.  For  the  case  of  hydrocarbon  "A,*’  which  approximates  an 
incompressible  material,  the  relationship  decreases  with  strain.  It  can  be  expected  that 
for  an  incompressible  material  such  as  hydrocarbon  "A, "  decrease  in  value  of  this  expression 
with  strain  will  be  maximum.'  Also,  as  the  material  deviates  from  incompressibility  the 
values  of  Aw/w0  will  approach  zero  and  then  become  positive.  This  can  be  seen  from  the 
values  for  the  plastisol.  The  data  would  Indicate  that  the  strain  dependence  of  A'W/w  is 
slightly  more  pressure  sensitive  than  that  for  Ar/i*c-  The  correlation  is  thus  more  com¬ 
plex  than  that  for  and  additional  work  at  more  pressures  and  strains  is  necessary 

before  any  definite  criteria  for  mass  burning  can  be  established. 

The  attempt  to  develop  a  test  specimen  suitable  for  measuring  burning 
rats  in  a  nonuniform  strain  field  was  unsuccessful.  The  contributing  factors  were: 

1)  The  mass  of  the  propellant  sample  had  to  be  relatively 
small  due  to  the  limitations  set  up  by  the  dimensions  of 
the  window  bomb.  This  created  problems  in  the  Photo- 
Stress  analysis,  since  the  fringe  patterns  were  so  close 
together  that  an  accurate  correlation  of  the  burning  rate 
with  position  was  difficult,  especially  at  higher  strains. 

2)  The  size  of  the  sample  also  requires  that  extreme  care  be 
taken  in  the  selection  of  the  PhotoStress  plastic.  The  end 
effects  at  the  outer  perimeter  of  the  specimen,  due  to  re¬ 
lative  reinforcement  of  plastic  and  material,  at  times 
extended  across  the  entire  specimen. 

3)  The  PhotoStress  plastics  used  all  broke  away  from  the 
vise  grip  walls  at  very  low  strains.  The  most  promising 
material  used  was  a  Hysoi  #2085  PhotoStress  plastic  with 
a  modulus  of  approximately  660  psi. 
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4}  Th8  technique  in  its  present  form  could  not  assure  that 
the  strain  patterns  induced  in  the  sample  would  not  be 
altered  as  the  sample  material  was  consumed.  The 
most  promising  specimen  evaluated  was  the  wedge. 

For  this  reason,  along  with  the  more  promising  results  being  generated 
in  the  uniform  strain  field  studies,  it  seemed  advisable  that  the  latter  effort  in  this  area 
be  discontinued  in  January  1963  until  a  better  definition  of  uniform  strain  effects  could 


?oe  achieved. 


IV,  PROPELLANT  SLUMP  ANALYSIS 

Propehant  structural  integrity  must  be  of  prime  concern  to  engineers  associated 
with  the  design  and  development  of  solid  propellant  rocket  motors  since  failure  of  a 
motor  system  may,  in  many  instances,  be  traced  to  deficioncies  in  the  design  and/or 
physical  properties  of  the  propellant  charge.  Effort  in  this  area  during  the  past  year 
has  been  directed  toward  the  development  of  a  means  by  which  the  effects  of  internal 
pressurization,  thermal  shrinkage,  and  axial  acceleration  loads  upon  propellant  physical 
integrity  could  be  predetermined.  Since  any  propellant  structural  investigation  must 
consider  the  complex  time-dependent  (or  viscoelastic)  properties  of  the  material  being 
examined,  a  limited  amount  of  time  and  effort  was  expended  in  the  area  of  propellant 
characterization.  The  second  phase  of  the  program  concerned  the  development  of  a 
method  for  analyzing  structural  problems  in  axisymmetric  grains  having  arbitrary  end 
geometries  and  straight-through  ports. 

Subsection  A  specifies  a  means  of  analytically  characterizing  the  viscoelastic 
properties  of  a  solid  propellant  by  linear  differential  operators  and  the  time-dependent 
elastic  modulus,  E(t).  The  method  and  its  accuracy  are  illustrated  by  examining  the 
time-dependent  properties  of  a  typical  hydrocarbon  propellant  (TP-H-1011). 

Subsection  B  presents  a  technique  for  resolving  structural  problems  in  axisymmetric 
grain  designs  subject  to  loadings  of  pressure,  temperature,  and  axial  acceleration.  The 
analysis  is  formulated  in  terms  of  two  stress  functions,  §  and  Y,  and  a  coupled  pair  of 
governing  elliptic  partial  differential  equations.  The  solution  is  sought  through  a  finite 
difference  approximation  to  the  governing  syscem.  The  pertinent  boundary  conditions, 
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together  with  their  corresponding  finite  difference  forms,  are  specified.  A  digital 
computer  program  based  upon  the  above  formulation,  which  is  being  completed,  is 
discussed. 

A.  CHARACTERIZATION  OF  PROPELLANT  TIME- DEPENDENT  MATERIAL 

PROPERTIES 

Problems  associated  with  predetermining  the  structural  integrity  of  Bolid 
propellant  rocket  motors  defy  resolution  for  many  reasons.  Among  these  are 
difficulties  associated  with  tho  inherently  complex  mechanisms  that  govern  the  response 
of  composite  solid  propellants  to  broad  spectrum  environmental  and  operational  loads. 
That  is,  present  composite  3olid  propellants  are  members  of  a  wide  class  of  polymeric 
materials  whose  structural  characteristics  are  typified  by  both  elastic  and  viscous 
response  when  subjected  to  given  tensile,  shear,  or  compressive  forces.  If  propellant 
structural  designs  arc  to  be  reliable,  they  must  necessarily  take  Into  account  the  complex 
physical  properties  that  these  highly  loaded  polymers  exhibit. 

1.  Determination  of  Linear  Viscoelastic  Model  Constants  from  Uniaxial 
Creep  and  Stress  Relaxation  Tests 

In  classical  elasticity  where  linearity,  isotropy,  and  homogeneity  are 
usually  assumed,  the  relationship  between  stress  and  strain  in  a  Hookean  body  is  given 
by 

a  -  Ee 

For  a  linear  viscoelastic  material  (l.e. ,  a  material  that  satisfies  Boltzmann's 
superposition  principle)  it  can  be  inferred  that  the  Laplace  transformed  stresses 


(3.1) 
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0 

are  related  to  the  Laplace  transformed  strains  by  an  analogous  equation  of  the  form 

cr  (a)  »  s  E(s)  «  (s)  (3, 2) 

The  transformed  elastio  modulus  E(s)  is  usually  associated  with  the  transfer  function 
of  a  generalized  Maxwell  mechanical  model,  as  illustrated  in  Figure  47.  This 
association  arises  naturally  from  the  Isothermal  viscoelastic  stress-strain  law  that 
governs  the  response  of  an  incompressible  linear  viscoelastic  solid  to  a  uniaxial  force. 

This  relationship  may  be  expressed  by  either  hereditary  integrals  or  the  following  linear 
differential  operator  form0 

P<r(t)  a  Q*  (t)  (3.3) 

where 


P  = 


and 


Q  = 


n 

E 

i  =»  o 


n 

I 

i  =  o 


Sti 


Q 


i 


Jr 

0 


ti 


The  employmentof  (3. 3)in  solutions  to  propellant  structural  problems  has  been  limited 

because  of  the  inherent  difficulties  that  the  determination  of  the  so-called  "model 

constants" presented.  ®  These  constants  are  not  available  explicitly  but  may  be 

obtained  implicitly  from  experimental  uniaxial  stress  relaxation  or  constant  stress 

8  11 

(creep)  data  wherein  both  stress  and  Btraln  are  known  functions  of  time,  ’ 

An  accurate  knowledge  of  propellant  relaxation  times  in  the  millisecond 
range  is  essential  to  the  meaningful  use  of  equation  (3,  3).  In  the  past, stress  relaxation 
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data  were  unavailable  from  a  reliable  experimental  apparatus  capable  of  measuring 
stress  decay  at  abort  f‘mes  while  eliminating  transient  dynamic  effects  in  the 
experimental  system  /or  these  times,  Difficulties  that  resulted  have  been  resolved 
at  Thickol  during  the  past  year,*®  Many  experimental  problems,  such  as  those 
mentioned  above,  were  overcome;  and  the  measurement  of  stress  response  for  finite 
input  strain  histories  has  progressed  so  that  now  error  in  empirical  data  does  not 
contribute'  significantly  td  the  over-all  error  in  the  analytical  formulation  of  stress 
analysis  solutions  based  upon  differential  operators, 

In  a  tensile  relaxation  test  having  a  finite  loading  time,  the  input  strain 
may  be  characterized  by  {see  Figure  48) 

«  t«>  3  V  -  Ao  (3 

In  this  case,  the  time-dependent  relaxation  modulus,  as  determined  from 
the  transf&r  function  of  the  mechanical  system  illustrated  in  Figure  47,  may  be  rep¬ 
resented  as  a  function  of  time  :s 

Eft)  *  Cn+2  ♦  8  Cje  ’  ^  (3 

.  ..  i  =  1 

where  n  is  some  conveniently  chosen  integer. 

.'he  time-dependent  stress  as  measured  during  the  relaxation  test  may  be 
characterized  as* 


♦See  Appendix  17 
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FIGURE  44.  INPUT  STRAIN  VERSUS  TIME 


(3.6 


c<t)  -  Ac 


An+2l  +  An+1  +  2  Ate 

1  »  1 


-  otjt 


^n+2(t-ti)  +  An+l  +  ^  Aje 

i  3  1 


-  q(t-tj) 


u  (t-tx) 


where  the  Aj's  in  (3.6)  are  related  to  the  in  (3. 5)  by 
°n+2  a  An +2 


n  Cj 

E  oT  a  Vn 
i  =  i  1 


(3.1 


-£t 


For  t>  ti 


Ai 


°  ll)  ~  A0  (  An+2tl  +  -  Aj 

l  =  1 


i-**  iiyait) 


(3. 


If  the  stress  response  in  relaxation  is  known  for  t  >  t<  and  (3. 8)  is  fit  to 
these  data,  it  can  be  deduced*  that  the  relaxation  modulus,  (3. 5), and  the  material 
constants,  Pj,  Q,,  are  specified  in  that  Cj,  Pj,  Qj  can  be  determined. 

Difficulties  inherent  in  the  fitting  of  Prony  type  series  <3. 8) to  experimental 
data  have  been  emphasized  elsewhere®  and  will  not  be  considered  here.  However,  the 
fit?.ingof(3.8)todata  where  the  decay  factors  dj  have  not  been  prescribed  previously 
gives  rise  to  a  system  of  transcendental  equations  that  must  be  solved  by  unwieldly 
algebraic  techniques,^  This  problem  may  be  averted  by  employing  a  collocation  method 
due  to  Schapery, 12  which,  in  essence,  states  that  in  fitting  viscoelastic  data  it  is 
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♦See  Appendix  IV. 


sufficient,  in  a  least  squares  sense,  to  specify  the  decay  factors  as 


(3.9} 


2.  Time  Dependent  Characterization  of  TP-H-1011 

The  analysis  given  above  has  been  used  to  investigate  the  time-dependent 
mechanical  properties  of  TP-H-1011,  a  hydrocarbon  solid  propellant  used  in  the 
Minuteman  propulsion  system.  Typical  data  from  uniaxial  stress  relaxation  tests 
conducted  on  this  propellant  are  illustrated  in  Figures  49,  60,  and  51. 

The  80°F  Isothermal  temperature  tests  disclosed  that  the  response  of 
TP-H-1011  could  be  adequately  represented  if  the  constants  Cj,  Pj,  and  Qj  were 
specified  as 


Jl. 


p6 

1 

Q6 

= 

9,794.4 

c6 

a 

7.38 

P5 

= 

505.55 

q5 

= 

522,234.5 

C6 

- 

77.8 

P4 

- 

2780.5 

q4 

= 

1, 474, 683. u 

C4 

sc 

119.8 

P3 

= 

1402. 9 

q3 

= 

501,832.026 

C3 

a 

153.9 

P2 

= 

70.132 

q2 

= 

16,563.32 

C2 

= 

434.6 

Pi 

= 

0.347 

Qi 

= 

55.522 

Cl 

= 

8,720.7 

Po 

- 

1.5  x  10"4 

Qo 

= 

0.  023 

Cs 

= 

145.0 

This  specification  is  equivalent  to  the  use  of  a  13-  element  Wiechert  model  (n  =  6).  The 
accuracy  of  the  method  is  best  illustrated  by  considering  Figure  52  where  data 
(stress)  from  a  relaxation  test  are  compared  with  the  results  obtained  from  the 
analytical  characterization. 
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FIGURE  4S-  RELAXATION  MOMJLVB  VERSUS  LOG  (TIME),  TP-H-1011 


FIGURE:  §0.  RELAXATION  MODULUS  VERSUS  LOG  (TIME),  TP-B-101I 
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FIGURE  a.  RELAXATION  MODULUS  VERSUS  LOG  (TIME),  TP-H-10I1 
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Log  (Time),  sec 

FIGURE  52.  STRESS  VERSUS  LOG  (TIME),TR-H-1011 


Although  these  results  are  fragmentary,  they  provide  basic  information 
concerning  the  variation  of  modulus  with  time.  This  information  is  necessary  to 
extend  the  elastic  results  of  the  finite  difference  elastic  analysis  described  in  sub¬ 
section  B. 

B.  STRESS  ANALYSIS  OF  AXISYMMETRIC  GRAIN  DESIGNS  HAVING  ARBITRARY 

END  GEOMETRIES  AND  STRAIGHT-THROUGH  PORTS 

Solid  propellant  rocket  motors  are  exposed  to  a  wide  range  of  loadings 
during  operation  and  storage.  Basic  among  these  are  loadings  due  to  internal  pressure 
and  acceleration  during  operation  and  propellant  changes  in  volume  due  to  thermal 
shrinkage  (expansion)  during  cure  and  subsequent  storage. 

Principal  effort  In  this  area  has  been  directed  toward  the  development  of  an 
elastic  analysis  for  axisyir.metric  grain  designs  having  arbitrary  end  geometries  and 
straight-through  ports.  The  propellant  charge  is  assumed  to  be  linear  elastic, 
i80tropic;  and  homogeneous  and  bonded  to  a  rigid  case  on  its  outer  periphery  head 
and  aft  ends.  The  stresses,  strains,  and  displacements  to  be  determined  are  those 
induced  by  either  internal  pressurization,  uniform  thermal  shrinkage,  or  a  constant 
acceleration  load  applied  in  the  axial  direction. 

The  initial  step  in  any  structural  analysis  of  a  solid  propellant  rocket  grain 
is  to  obtain  a  solution  of  the  elastic  field  equations  subject  to  given  boundary  conditions 
and  a  prescribed  external  loading.  In  particular,  if  the  grain  is  axisymmetric  and 
the  prescribed  loads  are  rotatlonally  symmetric,  the  governing  system  is  considerably 
simplified,  Jn  this  case,  the  analyst  has  several  methods1^  at  his  disposal  f^om  which 


hu  may  soak  the  solution.  However,  such  methods  present  Inherent  difficulties  that 
may  be  traced  to  the  basic  need  for  satisfying  the  mixed  boundary  conditions  typical  of 
this  class  of  problems.  An  approximate  numerical  technique  originally  developed 
by  Southwell1**  and  more  recently  extended  to  the  analysis  of  solid  propellant  rocket 
grains  by  Parr15  18>  ^  was  found  to  be  especially  applicable  to  the  solution  of  t’ 
problem. 

X.  Stress  Functions 

Elastic  bodies  having  forms  that  are  solids  of  revolution  and  loaded  so  that 
axial  symmetry  la  maintained  in  the  resulting  deformation  have  been  shown  to  be  amenable 
to  numerical  solution  using  stress  functions  and  relaxation  techniques.  14» 18  In 
particular,  for  the  conditions  described  above,  if  two  functions,  ♦  and  Y  ,  of  the 
spacial  coordinates  r  and  z  (but  not  #)  are  defined  in  such  a  way  that  they  satisfy  the 


the  equilibrium  and  compatibility  conditions  are  satisfied  if  the  stresses  are  defined  by 


(4.2) 
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(4.4) 


T  «  iO 

rz  r  8z 

I2  is  the  effective  body  force  due  to  acceleration  and  is  equal  to  the  load  in  g's 
multiplied  by  the  mass  density  of  the  propellant. 

The  associated  strains  and  displacements,  expressed  as  functions  of  $ , 
Y,  and  their  partial  derivatives  may  be  obtained  from  the  axisymmetric  stress-strain 
and  strain-displacement  relations  of  elasticity.  For  brevity  they  are  not  delineated 


here. 


The  relations l)-(4. 5)and  their  corresponding  strains  and  displacements 


may  be  generalized  to  include  the  effects  of  an  axial  acceleration  and  propellant 
expansion  due  to  thermal  shrinkage  during  cure  (see  references  16  and  18). 

2.  Formulation  of  Problem  in  Terms  of  Dimensionless  Quantities 
UBing  standard  notation*®*  we  define  (see  Figures  53  and  54). 


V  ~  *r 


11  "  b 


*9  = 


z  -  bIz 

Z  "  TT 


6  ^  —  6  7 

T1  -  z 


$  =  b^E  0 


V  Y 

TPtl  =  *z 


Y  =  b2E$ 


u  s  w 

n  K 
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OKECTION  OF  A 


154* 


FIGURE  53.  AXIAL  CROSS  SECTION  OF  CYLINDRH 
WITH  ARBITRARY  HEAD  AND  AFT  END 
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FIGURE  54.  NONDIMENSIONAL  COORDINATE  SYSTEM 


In  ordor  that  the  governing  system  £  ,  1),  stresses  (*i.2)-{4.5),  strains  and  displacements 
may  be  expressed  In  a  more  compact  nondimonsional  form.  These  may  now  be  written 
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On  a  Sat  head-end  boundary,  *\  was  specified  a?  zero.  This  boundary 


Side  D  fp  **  fg  {in)j 
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The  boundary  conditions  are  the  same  as  those  for  Side  B;  namely, 
Acceleration,  Thermal  and  Pressure  Loading 
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was  either  fixed  or  free.  For  a  fixed  end  condition,  the  displacements  are  zero  for 
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Acceleration,  Thermal  and  Pressure  Loadings 


If  X)*  ie  free,  the  boundary  conditions  take  the  form  of  specified 


stresses,  which  are; 
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the  boundary  conditions  for  each  pair  of  adjoining 
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The  stress  functions  i  and  ijt  must  satisfy  not  only  the  governing  differential 
equations  bat  also  the  boundary  conditions  specified  by  (6. 2)-(6, 21).  In  most  instances, 
these  general  boundary  conditions,  when  coupled  with  their  respective  representations 
(5.2)-$.  12),  will  be  used  as  stated.  However,  in  certain  instances  (e.g. ,  on  curved 
boundaries)  the  use  of  these  conditions  is  not  desirable  from  a  numerical  analysis 
standpoint.  To  obtain  numerical  representations  that  are  "nice, "  recourse  must  be 
made  to  physical  intuition  and  algebraic  manipulation  of  the  various  relations  of 
classical  elasticity. 

4.  Finite  Difference  Formulation  of  Solution 
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The  five-point  nodal  molecule  in  Figure  55  illustrates 
central  finite  difference  representation  of  (5, 2)  in  relation  to  node  0  is 
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If  one  or  more  points  of  a  nodal  molecule  lie  outside  the  boundary  of  the 
alterations  must  be  made  in  (7.  i).  Such  points,  aptly  termed  factitious, 


for  both  boundary  and  interior  nodes  when  these  nodes  are  £ 
nodes  by  a  boundary. 
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the  mic£3ection,  as  it  does  in  Figure  54,  three  possibilities  are  ap 
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i)  Ncotjs  2  and  3  are  fictitious; 


2)  Node  3  is  fictitious,  or 


3}  Node  2  is  fictitious. 


Similarly,  if  the  curved  segment  (aft  end)  of  the  boundary  lies  to  the  right 


of  the  mid-section,  as  it  dees  in  Figure  34,  three  possibilities  are  again  applicable  to 


the  nodes  centered  in  the  interior.  However,  if  node  2  is  fictitious,  possibility  (3) 


above  is  applicable,  and  only  the  distinct  cases  need  consideration; 


1)  Node  1  is  fictitious,  or 


2)  Nodes  1  and  2  are  fictitious. 


The  five  distinct  possibilities  are  illustrated  in  Figure  56, 


The  first  case  is  identified  by  (1)  in  Figure  56,  Without  lose  of  generality 


it  may  be  assumed  that  0  and  lji  can  be  expanded  about  node  0  in  a  two-dimensional 


Taylor  series  of  the  form 


X  =  X0  +  <n-  ib) 


( 8X\  ,  /SX\ 
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The  fictitious  nodes  2  and  3  must  be  replaced  by  their  corresponding 


boundary  nodes,  A  and  B,  The  points  1,  A,  4,  and  B  have  coordinates  (ri0  ;h,p0), 


(  rjj-  $Ah,  pc),  (rj0,  c0-h)  and  (r^,  pQ  +  fi3h),  respectively.  When  these  are  substituted 


into  (7.2)  (after  terms  having  factors  of  h  greater  than  three  are  truncated),  four 
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linear  equations  result  in  terms  of  the  first  and  second  derivatives  of  X  with  respeot 
to  p  itnd  f\  at  node  9  c 


X  \  *  y  k  (mS  „|£  fd2x) 

x‘  °'alHr^/0 

SB  -  X0  +  *£h  (tt\  ♦  ItaK?^ 

\a^/o  2  \?p“/0 

When  these  systems.  (7.  3)-(7. 4),  are  solved  for  the  first  and  second  derivatives  of  X 
and  these,  in  turn,  are  substituted  for  the  corresponding  derivatives  of  0  and  i)r 
in  (f,2}4  the  finite  difference  representation  of  the  governing  system  is 
Xj  sx  *  Kj  08+  K,  0A*  K4  04  +  X,  0O.  -  0 

%H+K2  *b'k3  *AtK4  *4"  %  *0 

-K,  0i-K2  0a*K6  «*„-<> 

where  the  K^’s  ate  constants  that  must  be  computed  for  each  interior  node  of  type  (1). 
In  this  instates  the  Kj*s  were 
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In  a  similar  manner,  finite  difference  expressions  for  the  governing 
system  for  Irregular  Interior  nodes  of  types  (2),  (3),  (4),  and  (5)  may  be  derived.  These 


{7.10) 

(7-11/ 


were: 

Type  (2) 

^  *  t1  '  2^]  *2  +  Ka  0A  +  [l  +  2^"]  04  '  [2  +  Ka]  «o  “  0 

Kl  »1*  *2*«3  *A+[l+2^-l*4  <7. 12) 

-  [2  *  Kel  *0  -  K,  -  K3  0a*K60o-O 

Type  {3) 

01  +  K2  $b  +  03  +  K4  04  +  K7  0o  a  0 
ijfl  +  K2  +  $ 3  +  ^4  ^4  +  K7 

-  0!  -  03  *  ?0O  =  0  (7.13) 

Typ_e_(4j 
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K2  ♦'  +  f1  -  itr]  *2  +  Rl‘  *3  +  I1  +  2?Tj  *4  ,7' 141 
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(7.1 


K3  0O  +  K2  0b  +  K1  03  +  ^4  04  +  0o  ■  0 
K3  to  +  K2  ^  K1  ♦a  +  **4  ^4  +  KS  t0 
-  K*  03  -  K3  0O  +  Kj  0o  -  0 


where 


JL  .  2  +  iil-J] 

SB  P0  SB 


Kj*  indicates  that  the  string  0-1  has  been  cut  by  a  boundary  at  C;  therefore,  A  must 
be  replaced  by  C  and  by  ®c  •  For  a  fully  loaded  head  end  and/or  a  partially  loaded 
conical  aft  end  (see  Figures  57-60),  similar  expressions  could  be  derived  for  use  in 
computations  at  irregular  Interior  nodes  created  by  these  additional  curved  boundaries, 
b.  Stresses.  8trains,  and  Displacements 

The  finite  difference  forms  of  the  equations  used  to  calculate  the  stresses, 
strains,  and  displacements,  (5.3)-(5.12),  are  (see  Figure  56): 

Regular  Interior  Nodes 

(^p)osiL^L)  [to +  d  '  v)0oJ  "  vP0  tiqZ  +  P0i  (7 

<u^)o = {[f  (T10,] s  -  Oo2|.  a  +  v)  (7. 

(c  P>°  =  f^'p  V2  {.o  ^2  ~  U)  +  P0  (!  “  V)  (02  -  04>  (7. 

-  2h  +  (1  -  V)  0oj|  -V110  Z  +  6 

(€e)o  =  %^[^oMi-v)0o]-v^  Z  +  6  (7. 
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Direction  of  Acceleration 


FIGURE  59.  NONDXMENSXONAL  COORDINATE  SYSTEM 
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FIGURE  60.  IRBERfTTJkR  nrmirat  u/v>w» 
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c.  Boundary  Conditions 

In  addition  to  satisfying  the  governing  differential  system  at  each  point  in 
the  interior,  the  stress  functions  must  also  be  compatible  with  the  boundary  conditions 
(see  IV-B-2)  imposed  by  the  loading  and  particular  geometry  being  considered.  The 
finite  difference  form  of  these  constraints  may  be  expressed  as: 

Side  A  (p  =  a) 

The  constraints  that  all  considered  loadings  place  on  the  stress 
function  on  this  boundary  are  summarized  by  two  relations: 

h  +  T  l1  +  h]}0o 
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(7.27) 


♦o  *  0 


(7.28) 


-171- 


Sides  B* 
of  the 

flat-end  cylinders  that  he  studied. 

Side  B  [p  ■  fj  (T|)] 


This  finite  difference  form  of  the  boundary  conditions  for  Side  A  and  those  for 
»"id  D+,  below,  are  the  same  as  those  given  by  Parr*8  for  the  free  boundaries 


On  thla  fixed  boundary  the  displacements  are  identically  zero;  and, 
using  (6. 6)  and  (7. 17),  one  condition  that  the  stress  functions  must  satisfy  is 


[to  +  U-v)*o]  'VPo^z  +p04  a  0 


<7.29 


When  applied  to  <7. 18),  the  second  boundary  condition,  (6  7),  does  not 
present  a  form  that  is  very  tractable  from  a  numerical  analysis  standpoint.  However, 
since  u^  is  invariant  with  respect  to  arc  length  on  this  boundary, 


3  0 


(7.30 


but 


=  i“2L  5a  +  lEn  & 


3u 


9u- 


3  s  3  T)  ds  dp 
Hence,  using  (7.30)  and  (7.31) 
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On  this  boundary 
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By  applying  (7.33)-(7.37)  to  (7,34),  &  more  tractable  boundary  condition  Is  obtained} 
namely, 

a  *  va.(6j,.(i.v)Ml.  x fa 

dp  dp  dri  dn 

*  J 

-  >  V(Y-H)].  -fii-  =  0 

K£  +V)  *■  J  1  +V 

Using  backward  differences,  in  finite  difference  form  (7.38)  may  be  expressed  as 
«Ch  [*o  "  t4  +  V(0O  -  04)] 
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The  equations  (7 , 29)  and  (7 . 39)  are  the  most  useful  form  of  finite  difference  relations 
that  can  be  formulated  for  Side  B. 

Side  B»  (H=  B+h) 
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2  P  o 
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Side  0(0  =  1) 


Since  up  *  0  on  this  fixed  boundary,  the  first  relation  that  the  stress 
functions  must  satisfy  (see  7.29)  la 


♦o  +  (1  “  v)  0O  - 


Vh0  2 
1  +  V 
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1  +  v 


is  also  identically  zero  on  this  boundary;  and  the  fact  that  axial .  ..placement  is 
invariant  with  respect  to  ri  here  implies  that  e^is  also  zero.  Thus,  after  (7. 1)  is 
used  to  eliminate  the  fictitious  values  62,  V  2  in  (7.21),  the  second  condition  that  the 
stress  functions  must  satisfy  is 

Side  D  paf2  (n) 

The  Side  D  displacements  are  zero,  as  on  Side  B;  therefore,  after 
the  geometric  terms  are  replaced  by  those  of  Side  D  and  subscript  1  is  substituted  for 
the  subscript  2,  the  boundary  conditions  (see  7.29  and  7.39)  are: 
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Side  D*  n  3  8  =  0 


Fixed  End  Conditions; 


Equation  (5,28),  when  coupled  with  {?,  17),  gives  the  first  condition 

♦oM*-v)0o  +  p0  5  *0 

po 

Sines  is  invariant  with  respect  to  p  on  this  boundary 


Now  the  strain-displacement  equation  that  relates  the  shearing  strain 


Y,  and  the  displacements  u„  and  u«ie  given  by 
Pvj  P  '< 
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When  (7.47),  (5.3),  and  (5. 12)  are  employed  and  the  fictitious  nodes  03,  t3  are 

eliminated  by  using  (7,1),  the  second  boundary  condition  takes  the  form 

01-20o-2  ^  +  (1  -  2V)  [ 0X  +02+04 

2 

-4<*o  +  Tp;  <«4-<2]-  tH  [5tvkZ]  =  0  (7.49) 

Free  End  Conditions: 

When  (5. 10)  is  integrated  along  this  boundary,  >)ls  known  at  each  nodal 
point  and  for  all  loadings  on  this  boundary  is  given  by 

.  <°l  -  °2>_  tr  (7.50) 
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From  equations  (6.21)  and  (7.25): 

^  -^3=0  (7.51) 
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The  result  of  eliminating  $g,  using  (7. 1),  is 

~~  (04-02>+*2+*4‘2*o  +  2  *1-2  "  0  <7.51 

«Pq 

The  relations  (?,  50}  and  (7,52)  are  the  most  useful  forms  of  the  boundary  conditions 
for  D*  when  this  side  is  unconstrained. 


Corner  Nodsa 


At  adjoining  boundaries  the  satisfaction  of  the  boundary  conditions 
defined  abo”e  consists  mainly  of  algebraic  manipulations  of  these  equations.  A  complete 
list  is  given  below. 
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6.  Calculation  of  Stresses,  Strains,  and  Displacements 

(7. 17)  -  (7.26)  are  the  equations  used  to  calculate  the  stresses,  strains, 
and  displacements  at  regular  interior  points.  The  equations  used  to  calculate  these 
quantities  at  other  than  regular  interior  points  are  listed  below. 

Irregular  Interior  Nodes 

Since  the  finite  difference  relations  used  to  calculate  up  and 
remain  the  same  for  all  nodes,  they  are  not  specified  under  the  different  types  below. 
Also,  since  several  of  the  relations  (7. 17)-  (7.26)  involve  derivatives  with  respect  to 
only  one  spacial  coordinate,  the  absence  of  a  representation  for  a  stress  or  strain  for 
any  one  type  of  interior  point  ii3ted  below  will  serve  to  Indicate  that  the  corresponding 
relation  (7.17)-(7.26)  is  to  be  used, 
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(7.7a) 

(7.73) 
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Type  (2L 
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1 

-  Same  as  (8.6) 
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Type  (5) 

The  equations  used  to  calculate  the  normal  stresses  an^ 

normal  strains  (  «e,  «n)  are  the  sarae  as  those  of  Type  with  6C  replacing  6  A* 

The  equations  used  to  calculate  the  shearing  stress,  ip^,  and  the  shearing  strain, 

y  ,  are  equations  (8,15)  and  (8.16),  respectively. 

Pr\ 

6,  Method  of  Numerical  Solution, 

If  the  finite  difference  equations  that  govern  the  elastic  state  ox  the 

propellant  grain  are  written  for  each  node,  the  values  of  the  stress  functions  at  each 

nodal  point  and,  hence,  the  dimensionless  stresses,  strains,  and  displacements  at 
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each  nodal  point  may  be  obtained  from  the  simultaneous  solution  of  the  resulting  set  of 
linear  equations  and  the  subsequent  employment  of  equations  (7,  l?)-(7.26)  and  (8, 1)~ 

(8, 18),  For  realistic  propellant  grains,  this  system  may  consist  of  as  many  as  2500 
equations.  Although  the  solution  of  such  a  large  linear  system  presents  certain 
difficulties,  there  are  a  number  of  methods  by  which  the  solution  may  be  obtained 
accurately  and  economically.  A  two-variable  relaxation  routine,  which  has  been  used 
extensively  in  the  past  by  Southwell15  and  Allen,20  was  used  in  the  digital  program 
based  upon  the  analysis  described  above.  The  method  is  fundamentally  point  iterative 
and  fairly  versatile  in  handling  the  large  linear  systems  common  to  axisymmetric  stress 
problems. 

Basically,  the  values  of  the  stress  functions  at  each  nodal  point  muse  be 
obtained  from  a  matrix  equation  of  the  form 

[ci,  j]  •  [xi]  =  [a,]  (t  =  1, n)  (j  =  1, n)  (9.1) 

where  j  is  the  coefficient  matrix,  [x^j  is  the  solution  vector  with  entries  0jj, 
j  and  j  is  the  constant  or  residual  vector. 

Initially,  all  entries  in  the  solution  vector  are  set  equal  to  zero.  New 
trial  values  for  the  stress  functions  0jj,  ^  j,  at  a  nodal  point  Pj  j  for  the  (n  +  1) 
iteration  are  obtained  from  the  formulas 
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where  (9. 2)  and  (9. 3)  are  used  for  nodal  points  in  the  interior  and  on  the  boundary, 
(n)  (n) 

respectively.  Git  j  are  the  residuals  obtained  by  evaluating  the  governing 
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finite  difference  relations  at  Pj  j  using  values  of  the  stress  functions  j,  Y|  j 
from  the  iteration,  The  Pj  are  relaxation  factors,  which,  ir.  general,  range  over 
[i,  2).  For  a  particular  geometry,  the  convergence  of  the  iterative  process  is  a 
strong  function  of  v  and  the  values  of  Pj  that  are  chosen.  As  v  decreases  in  value  from 
0. 5  {incompressible  propellant  material  as  measured  by  Poisson's  ratio),  the  number 
of  iterations  needed  for  convergence  decreases.  Values  of  w  1.7  are  near  optimum 
for  relaxing  the  interior,  while  0g  and  usually  vary  with  the  geometry. 

The  solution,  for  a  given  mesh  size  h,  is  assumed  to  be  "good"  when  the 
absolute  value  of  all  residuals  is  less  than  a  specified  number  5 : 

fi.jI<5  ’WjI  < 5 

The  value^of  5  indicates  the  degree  of  accuracy  of  the  solution  for  a 
given  meoh  size.  As  is  common  to  mo;  iterative  methods,  the  time  involved  in  solving 
any  given  problem  will  be  proportions  to  the  accuracy  desired. 
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A  description  of  the  complete  program  and  comments  on  its 
versatility  are  presented  in  subsection  7  below, 

7.  Description  of  Computer  Program  Stress  Analysis  of  Axisymmotrio 

Grain  Designs  Having  Arbitrary  End  Geometries  and  Straight-Through 

Ports 

The  block  diagram  in  Figure  61  illustrates  th6  digital  computer  program 
developed  to  analyze  structural  problems  in  axisymmetric  grain  designs  having 
arbitrary  end  geometries  and  straight-through  ports.  Although  the  existing  program 
is  not  capable  of  accounting  for  the  addition  of  a  head-end  web,  with  minor  modifications 
it  can  be  readily  adapted  to  handle  9uch  configurations  for  loadings  by  intei.  1  pressure, 
thermal  shrinkage,  and  axial  acceleration.  A  description  of  the  essential  features 
of  the  computer  program  (input,  calculations,  output)  is  given  below. 

Basically,  three  separate  routines  comprise  the  main  program: 

1)  Finite  difference  input  grid  generator 

2)  Relaxation  routine 

3)  Stress-strain-displacement  routine 

The  input  grid  generator  (see  Figure  61)  specifies  the  spacial  coordinates 
of  all  nodal  points  in  the  interior  and  on  the  boundary.  The  generator  also  classifies, 
by  type,  both  boundary  and  interioi  nodes.  Since  the  method  of  solution  is  fundamentally 
point  iterative,  the  generator  has  also  been  written  so  that  it  specifies  the  iteration 
indices  upon  which  the  relaxation  routine  is  based. 
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Tho  relaxation  routine  makes  use  of  systematic  over-relaxation  to  solve 
for  the  values  of  0  and  ijiat  all  nodal  points,  The  values  of  0  in  the  solution  vector 
[Xj]  are  improved  to  a  certain  degree  by  reduction  of  the  corresponding  entries  in 
the  residual  vector  [rJ  ,  These  values  of  0  are,  in  turn,  used  to  Improve  the  values 
of  i|f  in  the  solution  vector.  This  process  is  continued  until  ail  residuals  ere  less 
than  some  specified  number  dL.  When  all  residuals  are  less  than  dL,  the  values  of  0 
and  ^  which  are  in  storage,  are  used  to  calculate  the  stresses,  strains,  and 
displacements  at  each  node. 

Since  the  axial  displacement  at  any  nodal  point  involves  evaluating  integrals 
of  derivatives  of  the  stress  functions  0  and  t  with  respect  to  the  spacial  coordinate  P, 
numerical  integration  along  a  radial  line  P=  PQ  must  be  employed  to  calculate 
The  remaining  stresses,  strains,  and  radial  displacement  at  a  nodal  point  are  readily 
obtained  by  using  the  appropriate  finite  difference  equations,  (7. 17)-{7. 26)  and  (8. 1)- 
(8. 16),  which  specify  the  values  of  these  quantities  at  each  point  of  the  bounded  region. 
Since  these  expressions  involve  only  first  derivatives  (see  equations  5.3-5. 12),  these 
values  can  be  calculated  to  a  considerable  degree  of  accuracy. 

The  value  of  the  dimensionless  stresses,  strains,  and  displacements 
(u^  u0,  %  ,  eg,  ap,  crQ,  o^,  t rvj  ,YPti  )  at  each  nodal  point  of  the  bounded  region 
(see  Figure  61)  are  printed  as  output. 

Although  the  existing  analysis  and  computer  program  does  not  account  for 
the  viscoelastic  (time-dependent)  properties  that  most  propellants  exhibit,  the  results 
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(true  stresses,  strains,  and  displacements),  as  obtained  from  the  analysis  above,  may 
be  used  directly  in  analyzing  problems  involving  short  and  long  times  by  employing  the 
glassy  and  equilibrium  properties  of  the  propellant.  The  extension  to  time  space  of 
the  elastic  results  obtained  from  this  analysis  is  being  investigated.  Apparently,  a 
knowledge  of  the  variance  of  the  elastic  solution  for  a  range  of  the  elastic  modulus, 

E(t),  is  necessary  for  incorporating  the  viscoelastic  propellant  properties  into  the 
elastic  solution. 

This  program  will  be  used  to  investigate  the  effects  of  variation  in  the 
geometric  parameters  X  (1/d  ratio),  a  (a/b  ratio),  and  V  (ellipse  ratio)  upon  the 
stress-strain  fields  of  these  grain  configurations  when  they  are  subjected  to  loadings 
by  internal  pressure,  temperature  variation,  and  axial  acceleration.  These  results, 
in  the  form  of  design  charts,  will  provide  an  additional  structural  tool  to  design 
engineers. 
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V.  CONCLUSIONS  AND  RECOMMENDATIONS 


The  technical  discussions  presented  in  Sections  II,  III,  and  IV  outline  the  results 
obtained  from  the  studies  of  heat  transfer,  influence  of  strain  on  burning  rate,  and 
propellant  slump.  Each  of  these  Interdependent  studies  generated  specific  conclusions, 
which  are  given  in  this  section.  Recommendations  are  also  presented  based  on  the  work 
accomplished  on  this  contract. 

The  two-dimensional  heat  transfer  program  developed  under  this  contract  has 
extended  the  capability  to  predict  temperature  profiles  within  propellant  grains  to 
situations  more  closely  representing  the  actual  environments  experienced  by  today's 
propulsion  systems.  Subject  to  the  characterisation  by  experimental  means  of  the 
materials  being  considered  by  the  analysis,  the  program  will  handle  a  wide  variation 
in  the  temperature  dependency  of  thermal  properties.  The  program  has  the  demonstra¬ 
ted  ability  to  consider  either  a  single  thermally  symmetric  section  or  an  entire  motor 
cross  section  with  the  same  comprehensive  consideration  being  given  to  time- 
dependent  and  circumferentially  variant  external  environments.  The  above  capabilities 
and  versatility  are  not  known  to  be  available  in  any  other  single  program  throughout 
the  industry. 

Apart  from  the  above  capabilities,  the  resolution  and  apparent  accuracy  of  the 
results  were  extremely  rewarding.  This  is  not  always  the  situation  with  numerical 
solutions  of  this  type.  In  all  instances,  expected  effects  were  produced  and 
apparently  realistic  reproductions  of  the  anticipated  thermal  profiles  resulted. 


-192- 


In  these  respects,  the  two-dimensional  heat  transfer  program  provides  a  capability 
directly  applicable  to  subsequent  analysis  and  programs  for  ballistic  and  structural 
evaluation  of  grains  subjected  to  adverse  environments. 

Since  satisfactory  results  were  obtained  from  the  two-dimensional  heat  transfer 
program  and  an  extended  capability  was  shown,  the  following  recommendations  regarding 
the  use  and  possible  extension  of  the  analysis  arc  submitted. 

1)  Because  of  the  demonstrated  success  of  the  numerical  approach, 
the  basic  analysis  might  be  extended  to  consider  three- 
dimensional  configurations. 

2)  An  experimental  program  should  be  undertaken  to  ascertain  the 
necessary  temperature  dependency  of  thermal  properties  of 
propellants.  A  similar  effort  is  required  regarding  the 
coefficients  for  determining  the  heat  of  polymerization  during 
cure.  Only  upon  completion  of  these  efforts  can  the  full 
capability  of  the  program  be  realized. 

3)  Since  the  resolution  of  program  results  has  been  demonstrated, 
application  of  these  results  to  areas  other  than  the  ballistic 
analysis  should  be  investigated.  Specifically,  consideration 
should  be  given  to  the  development  of  a  thermal-structural 
analysis  to  evaluate  stress-strain  profile  within  propellant 
grains  subjected  to  adverse  thermal  environments.  This 
analysis,  of  course,  should  not  be  limited  to  thermal 
consideration  but  should  also  allow  for  pressure  and  dynamic 
loading,  since  the  moat  realistic  situation  involves  the 
superposition  of  all  these. 

These  recommendations  are  directed  toward  the  development  and  combination  of 
modular  program  systems  in  order  to  provide  an  integrated  system  for  complete  motor 
design  evaluation. 
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The  following  conclusions  resulted  from  the  study  program  to  advance  the  state 
of  the  art  of  solid  propellant  rocket  design  by  high-speed  electronic  digital  computation 
techniques: 

1)  The  Advanced  Grain  Design  Analysis  has  proved  to  be  an  ex¬ 
tremely  effective  tool  for  the  rocket  design  engineer.  The 
method  combines  conceptual  simplicity  with  extreme 
computational  accuracy  and  is  applicable  to  a  wide  class  of 
rocket  motor  designs. 

2)  The  Advanced  Grain  Design  Analysis  (Program  AGDA)  is  the 
only  computer  program  in  existence  capable  of  computing  the 
propellant  Inertia  tensor  components  as  they  vary  with  time. 

3)  Ability  of  the  Elkton-Moore  Method  to  consider  nonuniform 
thermal  gradients  throughout  the  solid  propellant  grain  has 
been  clearly  established.  This  is  the  only  known  technique 
available  that  is  amenable  to  such  considerations. 

4)  The  introduction  of  engineering-oriented  software  through  the 
Input  Generator  has  proved  to  be  of  tremendous  value  to 
computer  operations. 

Certain  areas  concerned  with  computerized  grain  design  by  AGDA  deserve  further 
consideration.  Thiokol  recommends  further  effort  as  follows: 

1)  While  Program  AGDA  computes  virtually  all  the  geometrically 
derivable  performance  parameters,  it  does  not  consider  such 
fundamental  parameters  as  chamber  pressure,  thrust, 
delivered  ISp,  or  total  impulse.  Therefore,  the  most  logical 
extension  to  Program  AGDA  would  provide  a  rigorous  Interior 
ballistic  analysis  capable  of  considering  the  following  effects: 
erosive  burning,  pressure  variations,  nozzle  erosion,  and 
ablation. 

2)  The  wide  class  of  solid  propellant  rocket  motors  that  may  be 
treated  by  Program  AGDA  represents  only  a  small  fraction  of 
total  imaginable  rocket  motors.  Additional  classes  of  rocket 
motors  that  might  be  considered  are  the  dendrite  or  forked 
wagon  wii-.el  and  the  conicel. 
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3)  Since  Program  AQDA  is  the  only  solid  propellant  rocket  grain 
design  computer  program  In  the  industry  today  capable  of 
considering  geometrically  dependent  propellant  properties, 
this  capability  could  be  extended  to  include: 

a)  Variation  in  burning  rate  with  three  -dimensional  grain 
temperature  variations 

b)  Variation  of  burning  rate  caused  by  three-dimensional 
strain  distributions 

4)  The  engineering-oriented  software  created  and  used  with  Program 
AGDA  should  be  applied  to  other  programs  requiring  large 
amounts  of  input  evolving  from  complicated  geometric 
configurations.  Specifically,  multi -dimensional  heat  transfer 

and  stress-strain  analyses  are  envisioned.  Further,  the  output 
of  such  programs  should  be  adapted  so  that  it  is  usable  directly 
by  Program  AGDA.  This  would  establish  a  relatively  complete 
solid  propellant  rocket  simulation  computer  program  system 
capable  in  the  most  general  sense  of  forecasting  the  failure  or 
predicting  the  performance  of  a  wide  variety  of  solid  propellant 
rocket  motors. 

5)  Future  effort  should  be  directed  toward  establishing  methods 
of  semi-automatic  o;  automatic  solid  propellant  grain  design. 

Fruition  of  such  concepts  as  the  solid  propellant  grain  design 
library,  a  catalog  of  solid  propellant  fuel  characteristics, 
automatic  combination  and  scaling  of  parts  of  rocket  motors 
from  the  solid  propellant  grain  design  library,  and  automatic 
generation  of  solid  propellant  fuel  characteristics  by 
consideration  of  the  past  history  of  every  solid  propellant 
rocket  motor  firing  using  this  propellant  should  be  realized. 

The  results  of  the  study  of  the  influence  of  strain  on  burning  rate  showed  that 
burning  rate  is  affected  by  strain  only  when  the  material  ceases  to  behave  as  an 
incompressible  material.  This  condition  occurs  in  filled  polymers  (such  as  composite 
solid  propellants)  when  the  tensile  strains  have  reached  a  level  that  causes  a  severing 
of  the  oxidizer-binder  bond.  This  phenomenon,  also  referred  to  as  "dewetting, "  can 
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occur  at  very  low  strains,  as  in  the  case  of  brittle  materials,  or  at  relatively  high 
strains  for  more  viscous  material.  Therefor  all  propellants  are  not  affected  to 
the  same  degree  by  strain. 

The  basis  for  a  definite  criterion  on  the  effect  of  strain  on  burning  rate  has  been 
shown  through  the  investigation  of  volumetric  change  index,  and  a  correlation  between 
percent  change  in  burning  rate  and  a  volumetric  change  index  is  given.  The  latter  term 
is  a  representation  of  the  departure  from  incompressibility  that  the  material  experiences 
due  to  strain.  For  the  linear  case,  an  empirical  relationship  has  been  shown  to  exist 
for  the  range  of  strains  and  pressures  and  for  the  propellants  tested.  The  fact  that 
the  results  obtained  here  correlated  quite  well  with  the  data  of  a  study  conducted  at 
another  agency  adds  further  support  to  this  conclusion.  For  mass  burning,  a  similar 
but  more  complex  criterion  has  been  shown  to  exist.  However,  both  cases  require 
more  study  before  a  general  correlation  for  all  classes  of  solid  propellants  can  be 
established. 

From  the  results  achieved  in  this  program,  certain  areas  ar^  recommended  for 
future  work. 

1)  Examine  a  greater  variety  of  propellants  (including  double-base 
types)  over  wider  ranges  of  strain  and  pressure  to  establish  a 
firm  correlation  among  burning  rate,  strain,  and  chemical 
structure. 

2)  Investigate  the  effect  of  temperature  upon  the  relationship 
between  strain  and  burning  rate, 

3)  Determine  the  extent  to  which  compressive  strains  influence 
burning  rate.  If  a  similar  correlation  could  be  extended  to 
Include  compressive  as  well  as  tensile  strain,  a  convenient 
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input  for  use  in  grain  design  computer  programs  for  ballistic 
performance  predictions  could  be  realized.  For  example,  if 
the  thermal  history  of  a  motor  grain  and  the  thermal  expansion 
properties  of  the  propellant  are  known,  the  strain  geography 
can  be  determined.  From  this  information,  along  with  the 
knowledge  of  the  strain  dependence  on  the  burning  rate,  the  net 
burning  behavior  for  the  motor  can  be  computed.  This  becomes 
increasingly  important  as  the  burning  rates  of  propellants  are 
increased. 

The  results  achieved  here  are  not  limited  in  their  use  to  the  solid  propellant 
industry.  The  experimental  logic  employed  should  also  apply  to  the  pyrotechnic  industry 
or  to  any  Industry  when  the  control  of  burning  rate  is  important. 

The  computer  program  conceived  under  this  program,  'Stress  Analysis  of 
Axisymmetric  Grain  Designs  Having  Arbitrary  End  Geometries  and  Straight-Through 
Ports,"  provides  a  suitable  method  for  predetermining  the  response  of  these  grain 
designs  to  pressure,  thermal,  and  axial  acceleration  loads. 

It  is  recommended  that  this  analysis  be  extended  to  encompass: 

1)  A  thorough  investigation  of  the  effects  of  the  parameters  v  ,  X 
(L/D  ratio),  a  (a/b  ratio),  and  -tf  (ellipse  ratio)  upon  the  stress - 
strain  field  in  axisymmetric  grain  designs  having  arbitrary  end 
geometries  and  straight-through  ports 

2)  A  modification  of  the  computer  program  in  order  to  incorporate 
the  effects  of  head-end  webs  and  arbitrary  aft  cutbacks 

3)  A  parameter  study  to  examine  the  effects  of  the  variables  in 

(1)  upon  the  strain  field  in  motors  having  fully  loaded  head  ends 

4)  Verification  of  the  studies  in  (3)  by  means  of  a  correlation  with 
actual  stresses  and  strains  in  experimental  models  having  full 
head-end  webs 

5)  Incorporation  of  the  above  recommendations  Into  engineering 
design  charts 
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TWO-DIMENSIONAL  HEAT  TRANSFER  PROGRAM 
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Symbol 

A.  GLOSSARY 

identification 

Dimensions 

A 

Experimental  Constant  for 

None 

B 

Heat  of  Polymerization 

Experimental  Constant  for 

None 

cp 

Heat  of  Pol,  mertzatlon 

Specific  Heat 

BTU/lb-°R 

E 

Heat  Flux 

BTU/hr-ft2 

Fo 

Shape  Factor 

None 

Internal  Film  Coefficient 

BTU/hr-ft2-*  R 

ft 

ho 

External  Film  Coefficient 

BTU/hr-ft  -R 

Kl.i 

Thermal  Conductivity- 

BTU -ft/hr  -ft2  -°R 

n 

Experimental  Constant  for 

None 

Pj,k 

Heat  of  Polymerization 

Polymer  Fraction 

None 

Q 

Heat  of  Polymerization 

BTU/lb 

R 

Gas  Constant 

BTU/lb-°R 

rU 

Radius 

In. 

T 
x  c 

Internal  Convective  Temp. 

°R 

Tu. 

Initial  Temp. 

0  R 

Iq 

External  Convective  Temp. 

°R 

^rad 

Radiant  Temp. 

0  R 

FORTRAN 

EN1 

EN2 

CP 

E 

GEOM 

HGC 

HGCS 

1HRCON 

EN3 

ALF(4) 

HOP 

R 

ALF(5) 

TEIN 

ALF{1) 

TECON 

TERAD 
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A',  GLOSSARY  (Continued) 


Symbol 

Identification 

Dimensions 

FORTRAN 

A  r 

Forward  Difference  In 

Radial  Direction 

In 

ALF(3) 

A  9 

Time  Increment 

Min 

DELTAT 

A  0 

Forward  Difference  in 
Circumferential  Direction 

Rad 

ALF{2) 

8 

Time 

Min 

TIME 

0 

Density 

lb/ft3 

RKO 
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B.  LjSl ^FUNCTIONS  AND^UBRQUTINE5..ITilA.NSEER..V.1SLCTORS) 
USED  BY  TWO-DIMENSIONAL  BEAT  TRANSFER  (LISTED 
SEPARATELY  FOR  PARTS  I  AND  ID 

Main  Program  -  Part  I 

TAPES1' 

EXIT** 

LOAD*** 

Main  Program  -  Part  II 

SQRTF 
COSF 
SINF 
TAPES* 

EXIT** 


1 

*  Subroutines  and  Functions  for  Part  II 

j 


Subroutines  for  Part  I  (no  iUnctions) 


*TAPES  subroutine  is  unique  to  the  Elkton  IBM-7070  System  but  can  fee  easily  ^avisad 
to  any  other  system.  The  tape  designation  at  the  Elkton  division  is  as  follows: 


FORTRAN  Tape  Designation 

Symbolic  Representation 

Use 

1 

None 

FORTRAN  Package 

O 

« 

IN 

Loader 

3 

LET 

Output 

4 

INTER  C 

Scratch 

5 

INTER  B 

Scratch 

6 

INTER  A 

Scratch 

This  may  be  revised  by  altering  the  equivalence  statement  to  the  specifications  of  another 
system.  As  an  example,  if  the  input  for  a  system  was  FORTRAN  tape  5  and  output  was 
FORTRAN  tape  6,  the  equivalence  statement  would  be  as  follows: 

EQUIVALENCE  (NTAPES(2),  INTERC),  (NTAPES  (3),  INTER  B),  (NTAPES  (4), 
INTER  A),  (NTAPES  (5),  IN),  (NTAPES  (6),  LET). 

**EXIT:  The  EXIl  subroutine  is  unique  to  the  Elkton  system.  It  is  not  required  for 
application  of  the  program  to  the  IBM;7090,for  instance. 

***LOAD;  Subroutine  LOAD  is  aroutlne  used  to  provide  for  multiple  core  loads.  This 
will  not  be  needed  In  the  IBM-7090  but  the  following  cards  will  have  to  be  deleted: 
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Part! 


185  CALL  LOAD 
END 

Part  II 

All  the  cards  that  contain  DIMENSION,  EQUIVALENCE,  and 
COMMON  statements. 


The  program  will  be  compiled  as  one  deck,  and  not  as  a  two-part  program. 


SYMBOL  DELTAT  TOT  I  RHO I 


1-7 


SHEET 


♦Phi  line  separating  two  external  boundary  zones. 
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I-iO 


1-12 


To  be  used  ior  each  pom*  considered  beginning  with  0  - 


definition  of  code  word 


WORD  1: 

COLUMNS  1-5: 

00000 

-  None  of  the  conditions  listed  below 

00001 

-0-0 

. 

00002 

-  0  a  0  max 

WORD  2: 

COLUMNS  6-10: 

00000 

-  None  of  the  conditions  listed  below 

00001 

-  Internal  point 

00002 

-  Internal  surface  point 

00003 

-  Propellant-case  interface 

00004 

~  External  surface  point 

WORD  3: 

COLUMNS  11-15: 

00000 

-  None  of  the  conditions  listed  below 

00001 

-  Entering  propellant  on  constant  0 

00002 

“  Leaving  propellant  on  constant  0 

WORD  4: 

COLUMNS  16-20: 

00000 

-  None  of  the  conditions  listed  below 

00001 

-  Entering  propellant  on  constant  r 

00002 

-  Leaving  propellant  on  constant  r 

00003 

-  Entering  and  leaving  propellant 

on  constant  r 

I— .13 


A-C589 


Page  1  of  4 


D.  FLOW  CHARTS 


I  Prelim, 
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Two-Dimensional  Heat  Transfer  Input  Format  (Revised) 

CPE-70129-A 


E.  TWO-DIMENSIONAL  HEAT  TRANSFER  INPUT  FORMAT  (REVISED) 


1  08-1 

1-60 

>—* 

Each  word  contains  the  beginning  of  each  constant  radial  line. 

NPHIB  (30) 

WORD 

COLUMNS 

FORMAT  1 

CARD 

3 

-j 

0 

eo 

i 

1-18 


WORD  1  2  3  4  5  6  7 _ 8 

7777771  TTYo  11-20  21-30  31-40  41  50  51-60  61-70  71-80 

■.QLVJMN5  . . . . — - - - - - - - - 

FORMAT  F10.0  F10-0  F10.0  F10-0  F10.0  F10.0  F10.0  F10. 0 

A0's,  forward,  for  each  radial  line  beginning  with  constant  phi  line  eq^ai  to  zero. 


1-19 


1-23 


1-25 


COLUMNS 

FORMAT 


WORD  1 : 


WORD  2: 


WORD  3: 


WORD  4: 


DEFINITION  OF  CODE  WORD 


COLUMNS  1-5:  00000  > 

00001  - 
00002  - 

COLUMNS  6-10:  00000  - 

00001  - 
00002  - 
00003  - 
00004  - 

COLUMNS  11-15:  00000  - 

00001  - 
00002  - 

COLUMNS  16-20:  00000  - 

00001  - 
00002  - 
00003  - 


None  of  the  conditions  listed  below 

0-0 

0  *  0  max 

None  of  the  conditions  listed  below 
Internal  point 
Internal  surface  point 
Propellant-case  Interface 
External  surface  point 

None  of  the  conditions  listed  below 
Entering  propellant  on  constant  0 
Leaving  propellant  on  constant  0 

None  of  the  conditions  listed  below 
Entering  propellant  on  constant  r 
Leaving  propellant  on  constant  r 
Entering  and  leaving  propellant  on 
constant  r 
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ELLIPTICAL  CASE.  CONSTANT  WEB 

Quita  often  a  specific  mission  may  require  the  design  of  a  grain  having  a  neutral 
trace,  a  high  impulse,  and  a  fairly  short  burn  time.  If,  in  addition,  there  is  a  low 
L/D  requirement  imposed,  the  problem  becomes  one  of  significantly  Increased  complex¬ 
ity  relative  to  the  straightforward  task  of  designing  a  straight-through  star  having  a 
neutral  trace.  Once  we  become  aware  of  the  considerable  influence  of  end  effects  on 
the  motor  performance  characteristics,  we  realize  that  we  must  have  a  progressive 
head  end  grain  design  to  counterbalance  the  inherent  regressivity  of  the  aft  end.  While 
loading  density  may  also  be  a  problem,  we  immediately  realize  that  the  elliptical  case 
head  design  will  be  more  effective  than  that  of  the  spherical  design,  since  the  portion 
of  our  over  -all  length  relegated  to  the  head  and  aft  ends  will  be  considerably  less  when 
using  the  former.  Unfortunately,  we  have  just  eliminated  the  end  case  design  directly 
amenable  to  mathematical  solution  in  closed  form.  Fully  aware  of  the  worth  of  a  uni¬ 
form  thickness  head-end  web,  we  endeavor  to  find  a  generalized  mathematical  relation¬ 
ship  between  the  surface  which  is  normal  to, and  a  fixed  distance  from, our  elliptical 
case  for  purposes  of  head -end  evaluation. 

An  iterative  scheme  was  derived  and  implemented  for  the  calculation  of  points  on 
this  resulting  nonelliptical  surface.  The  following  is  a  presentation  of  the  derivation 
by  Mr.  V/.  G.  Andrews.  For  this  purpose  it  is  sufficient  to  consider  an  social  cross 
section,  as  shown  in  the  figure  that  follows. 
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CROSS  SECTION  OF  SOLID  PROPELLANT  ROCKET  HEAD  END 
TO  ILLUSTRATE  MATHEMATICAL  CONSTRUCTION  OF  A  CORE  HEAD  SURFACE 
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If 


ia  the  elliptical  esse  equation,  then  the  equations  of  the  ellipse  in  parametric  form  are 
given  by, 

R0  =  a  cos  0 

0  s  9  s  it 


ZQ  =  b  sin  0, 

where  9  is  the  eccentric  angle  of  the  point,  P(RC,  Z^),  on  the  ellipse. 
Differentiating  these  parametric  equations,  it  is  seen  that 
dRc 


d0 


-a  sin  8 


- —  =  b  cos  0. 

d9 

Hence  the  slope  of  the  tangent  lire  to  the  ellipse  at  any  point  (Rc,  Zc)  =  (a  cos  9,  b  sin  9) 
is  given  by, 

dZc  _  b  cos  9 

n  - - - , 

dRc  a  sin  9 

It  follows  thac  the  slope  of  the  normal  is 


a  sin  9 
b  cos  8 


n-3 


Applying  the  two-point  slope  for  multi 

a  sin  9  _  Z  -  b  sin  9 
b  cos  9  R  -  a  cos  9 


Let 

R  -  a  cos  9  Z  -  b  sin  9  _  ^ 

b  cos-  9  a  sin  9 


Then 

R  -  a  cos  9  3  b  cos  0  *  t  -  R  -  Rc 
Z  -  b  sin  9  =  a  sin  9  •  t  =  Z  -  Zc 

and 

D2  -  (Z  -  Zc)2  +  (R  -  Rc)2 

9  2 

D2  =  (a  sin  9  •  t)Z  +  (b  cos  9  •  t)  , 

where  D  is  the  constant  web  thickness. 

This  implies  that 

+  _  D  _ 

*  (a2  sin2  9  +  b2  cos2  9)1'2 

Since  c2  =  a2  -  b2,  where  c  is  the  distance  from  the  center  of  the  ellipse  to  the  focus, 
the  expression  inside  the  radical  simplifies  to 

(a2  -  c2  cos2  9)  ^ 

An  examination  of  the  equations, 

R  =  (a  +  bt)  cos  9 

0  *  0  stt 

Z  -  <b  +  at)  sin  0 

11-4 


Sr. 


indicates  that  the  minus  sign  in  the  above  expression  for  t  is  applicable  for  the  given 
range  of  0,  and  in  fact  for  0  *  0  *  2  rr  , 

It  is  clear  that  the  coordinates  of  any  point  cf  the  burning  surface  are  given  in 
terms  of  the  auxiliary  variable  8  by  the  equations, 


E  =  a  - 


(a2  -  c2  cos2  0)1'2 


,  2  2  2  1/2 
fa*  -  c  cos  0) 


cos  8 


sin  8  . 


0  s:  0  s 


It  follows  that  any  attempt  to  solve  R  =  h  (0)  or  Z  =  g  (0)  for  9  will  ultimately 
lead  to  the  evaluation  of  a  fourth  degree  equation  in  cos  0. 

The  preceding  analysis  includes  the  special  case  of  the  sphere  (a  =  b)  and  it 
follows  that  analytical  expressions  for  surface  area,  volume,  etc. ,  can  be  obtained 
from  the  parametric  equations,  although  these  expressions  will  in  general  be  complex. 

Implementation  of  the  above  functional  relationships  in  the  Advanced  Grain  Design 
Analysis  Program  was  accomplished  by  the  determination  of  bounding  values  of  0  on  the 
eccentric  circle, 

R2  +  Z2  =  a2  . 

Following  the  standard  program  technique  for  the  analysis  of  the  head  end  web 
surface  regression  pursuant  to  construction  of  a  core  head  mesh,  an  initial  value  of 
R  on  the  core  head  is  obtained  (zero  initially).  The  angle  associated  with  thi3  coordi¬ 
nate  on  the  eccentric  circle  is  then  determined  from 


9,  =  cos-1  (  —  . 

la 


11—5 


If  the  R  coordinate  value  is  thon  incremented  by  a  web  thickness  and  the  respective 
angle  is  evaluated  from 

<J2  ■  cos”1  „(R  t_P).  t 
a 

it  is  Immediately  evident  that  the  desired  angle,  9,  is  bounded  by  the  relationship 
(see  the  figure  on  pago  n-2), 

0|  >01  02  « 

Since  the  actual  bounds  on  0  have  been  determined,  it  is  an  extremely  simple  matter 
to  calculate  the  precise  value  of  the  angle,  to  any  desired  degree  of  accuracy,  by  use 
of  the  "Method  of  False  Position"  (iteration  technique). 

Worthy  of  note  is  the  fact  that  in  this  specification  of  the  bounds  on  9  in  terms  of 
the  known  parametric  values,  for  complete  generality  the  range  of  evaluation  is  an 
absolute  minimum. 
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PICTORIAL  DEMONSTRATION  OF  THE  INSIGNIFICANCE 
OF  SURFACE  NORMAL  BENDING  RESULTING  FROM  A 
NONUNIFORM  GRAIN  TEMPERATURE 
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PICTORIAL  DEMONSTRATION  OF  THE  INSIGNIFICANCE  OF 
SURFACE  NORMAL  PENDING  RESULTING  FROM  A 
NONUNIFORM  GRAIN  TEMPERATURE 
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determination  of  linear  viscoelastic  model  constants 
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DETERMINATION  OF  LINEAR  VISCOELASTIC  MODEL  CONSTANTS 


FROM  UNIAXIAL  CREEP  AND  STRESS  RELAXATION  TESTS 


The  mathematical  model  constants  in  the  linear  differential  operator  specification  of 
time-dependent  tensile  stress,  <r  (t),  to  the  corresponding  time-dependent  tensile  strain, 


n  j, 

I  Pi-i 
i  =  o  St 


P^i  ]^(t)»[E  Qjli  ]  *(t) 

at1/  V  i  =  o  btl  I 


(IV-l) 


are  not  available  explicitly  but  may  be  obtained  implicitly  from  experimental  data  in 
which  both  stress  and  strain  are  known  functions  of  time.  It  is  clear  that  (TV-1),  after 


normalizing  with  respect  to  Pn,  can  be  written  as 

- +  Pji-1  — — —  + . *-  Pi  -S—  +  o  I 

L  at  atn_i  at  J 


->n-i 

a  .  .  n  a 


rr(t)  = 


(IV-2) 


Q  +■  Q  .  2 _  h  .  - - +  Q.  __2_  +  Q_  e(t) 

^n  n  m-1  n-l  1  ° 

.  at  at  at 


In  a  stress  relaxation  test,  the  input  strain  can  be  characterized  by  a  ramp  type 
loading  function  (see  the  figure  that  follows)  of  the  form 


e(t)  =  AQt-A0  (t-tj)  u(t-tx) 


(IV-3) 


where 


1  ts-ti 


Therefore,  the  stress-strain  relationship  (IV-2)  may  be  expressed  as 


h 


FIGURE  48,  INPUT  STRAIN  VERSUS  TIME 


IV-2 


inlL  + . +  P1A+Pn  c(t)a 

U"  "'St”-1  '«  °J 

%  nr  *  Q"-l  ~ZT  * . .  «1  ^  *  %  [Aot-Aoft-tO.utt-til] 


(IV -4) 


Taking  the  Laplace  transform  of  (IV -4),  it  is  clear  that  in  s-space 


Qn^  tQn-l  sn-U  ...  Qt3  +  Q0  [  Aq  :  A0e' 
sn  +  Pjj.j  Sn-1  +  .  .  .  PjS  +  P0  J  L  S2  S2 


(IV -5) 


Without  loss  of  generality,  the  factors  of  the  denominator  in  (IV-5)  may  be  assumed  to 
be  linear  and  distinct.  Therefore,  since  the  rational  function  of  s  is  proper,  it  may  be 
expanded  in  partial  fractions  to  give 


,  4  A  An+2  Ap+i  «  At  p,  -tisi 

a(s)  -  A0  -* — >+ - +  E  - «.  •  I-e  1 

a  i=i(s  +  ai)J  1  J 

The  inversion  of  (IV-6)  via  contour  integration  yields 

a(t)  -  A0  An+2t  +  An+1  +  Aieait  +  . . .  Ane"®^ 

u  ..." 


(IV-6) 


_Ao  An+2<t_tl)  +  An+1  +  Ale 


an  expression  for  a  in  time  space. 


-Vt-ti) 


+  . . .  +  Ane 


•  [uft-ti)] 


(IV -7) 


If  (IV -7)  were  fit  to  data  obtained  from  a  stress  relaxation  test,  the  constants  a^, 

Aj  would  be  known  and  the  model  constants  Pj,  Qj  could  be  directly  determined  from 

n 

the  relations  -^n-1  =  ^  (a) 

Pn-2  =  Va> 


(IV-8) 


Pn.u  =  K  (a) 


IV-3 


and 


„  n  n  n 

Qn  =  AtKx  (ax  =  0)  +  A2K1  (a2  =  0)  +  . . .  +  AnKx  (an  =  0) 

* 

*  Vik?<“>+V2  (iv'9) 

I 

Vk  *  AlKkn  (al  “  0)  +  A2Kk+l  <a2  =  °)  *  '  ■  •  ■  *  AnVl  (“n  ’  °> 

+  An+1  ^k+1  ^ ^  +  An+2  ^k^ 
n 

where  K-^  is  defined  to  be  a  function  of  (a)  in  that  it  is  the  sum  of  all  products  formed 
by  taking  the  a's  i  at  a  time  and  multiplying  them  together.  The  above  relations 
follow  immediately  when  (IV-6)  is  cleared  of  fractions  and  the  results  are  equated  to 
the  corresponding  polynomials  of  s  in  (TV~5).  For  example,  if  n=3,  then 

p2  ~  <al  +  a2  +  a3) 

pl  “<<*1*2  +  ala3  +  a2a3>  (TV-10) 

Po-  (aia2^3) 
and 

Qg  —  Aj  (Ct2  +  Otg)  +  Ag  +  CCg)  +  Ag  (Ot.^  +  Q-g)  +  A4  (Ctj  +  &2  +  &g)  +  Ag 

Q2  “  Ai  (a2a3)  +  A2  (a1a3)  +  Ag  (o^ag)  +  A4  (c^ag  +  axa3  +  a2a3) 

+  Ag  (ax  +  a2  +  a3)  (TV-U) 

=  A4  (CC|CC2CCg)  +  Ag  (CX^OCg  +  CC-^  Ctg  +  G.g&3) 

Q0=  A5  (CtiClgag) 


IV-4 


